Problem:
天体観測

ある晴れた夜の帰り道、空を見上げるとそこには無数の星。強く輝く星々、かすかに見える星々、いろ
いろな明るさの星々が夜空を彩っています。

あなたはふと思いました。この星空に自分の星座があったらどんなに素敵だろうと。そこであなたはひとつの基準を決め、部屋の窓から見える星々を適当につなげて自分の星座を作ることにしました。その基準とは、「どの２つの星を選んでも、それらの明るさの差がある定数 d 以下になるような星のグループを作り、その中で一番大きいものを自分の星座にしよう！」というものです。例えば、図のような窓から見える夜空を考えてみましょう（外側の長方形は窓枠です）。

この夜空には、明るさがそれぞれ 1，12, 2，4，1，8, 3，5，4 の９つの星がありますが、d = 2 とすると例えば以下のような３つの星座 A, B, C ができます。

大きさが 12 の星座 A

大きさが 10 の星座 B

大きさが 16 の星座 C

星座の大きさを次のように決めることにしました。ある星座の星をすべて含むような、窓枠に平行な辺からなる長方形のうち、面積が最も小さいものを選びます。この長方形の面積をその星座の大きさとします。例えば、上の夜空では星座 A, B, Cの大きさはそれぞれ12, 10, 16になるので、星座Cが最も大きい星座となります。

N 個の星の位置と明るさ、および整数 d が与えられたとき、一番大きい星座の面積を求めるプログラムを作成してください。星の位置は窓枠の左下隅を原点とした座標で与えられ、軸は図のような向きとします。星座を構成する星が１つの場合や、星々が軸に平行な直線上にある場合は、その星座の面積は 0 となることに注意してください。

入力

入力は以下の形式で与えられる。

N d
x1 y1 b1
x2 y2 b2
:
xN yN bN

1 行目に星の数 N (1 ≤ N ≤ 200000) と整数 d (0 ≤ d ≤ 109) が与えられる。続く N 行に、i 番目の星の座標を表す整数 xi (0 ≤ xi ≤ 2000000) と yi (0 ≤ yi ≤ 2000000)、明るさを表す整数 bi (0 ≤ bi ≤ 109) が与えられる。入力される星の座標はすべて異なる。

出力

一番大きい星座の面積を１行に出力する。

入出力例

入力例

9 2
1 1 1
1 5 12
2 3 2
3 2 4
4 4 1
5 1 3
5 3 8
6 5 5
7 2 4

出力例

16

Solution in C++11 (judge verdict: Time Limit Exceeded)

#include <iostream>
#include <algorithm>
#include <map>
using namespace std;

class Star
{
public:
	long long int x, y;
	int b;
	int boundary;
};
bool operator > (Star s1, Star s2)
{
	return s1.b > s2.b;
}
bool operator < (Star s1, Star s2)
{
	return s1.b < s2.b;
}

int n, d;
Star star[200000];
long long int maxX = 0, maxY = 0, minX = 0, minY = 0;
long long int ans = 0;

void calcBoundary()
{
	int boundary = n - 1;
	star[n - 1].boundary = boundary;
	for (int i = n - 2; i >= 0; i--)
	{
		if (star[i + 1].b - star[i].b > d) boundary = i;
		star[i].boundary = boundary;
	}
	
	for (int i = 0; i < n; i++)
	{
		Star s;
		s.b = star[i].b + d;
		boundary = upper_bound(star + i, star + n, s) - star - 1;
		star[i].boundary = min(star[i].boundary, boundary);
	}
}

void solve()
{
	int boundary = -1;
	for (int i = n - 1; i >= 0; i--)
	{
		if (star[i].boundary != boundary)
		{
			boundary = star[i].boundary;
			ans = max((minX - maxX)*(minY - maxY), ans);
			maxX = star[i].x;
			minX = star[i].x;
			maxY = star[i].y;
			minY = star[i].y;

			for (int j = i + 1; j <= star[i].boundary; j++)
			{
				maxX = max(star[j].x, maxX);
				maxY = max(star[j].y, maxY);
				minX = min(star[j].x, minX);
				minY = min(star[j].y, minY);
			}
		}
		else
		{
			maxX = max(star[i].x, maxX);
			maxY = max(star[i].y, maxY);
			minX = min(star[i].x, minX);
			minY = min(star[i].y, minY);
		}
	}

	ans = max((maxX - minX)*(maxY - minY), ans);
}

int main()
{
	cin >> n >> d;
	for (int i = 0; i < n; i++) cin >> star[i].x >> star[i].y >> star[i].b;
	
	sort(star, star + n);
	calcBoundary();
	solve();

	cout << ans << endl;

	return 0;
}