CodeForces 528 B. Clique Problem

Calendar Clock iconCalendar Clock icon

Code Forces

目次

# 問題

https://codeforces.com/contest/528/problem/B

一直線上に並んだNN個の点の座標xix_iが与えられます.
お互いにある条件を満たす2つの点同士をすべてつなげる時、お互いに直接繋がった点の最大の点の数(いわゆる最大クリーク)を答える問題です.
条件とは、点xix_iにはそれぞれ重みwiw_iが与えれていて、2点xi,xjx_i, x_jxixjwi+wj|x_i - x_j| \geq w_i + w_jを満たすことです.

# 入力

N
x_1 w_1
x_2 w_2
...
x_N w_N
  • N - 点の数(1N2000001 \leq N \leq 200000).
  • x_i - 点の座標(0xi1090 \leq x_i 10^9).
  • w_i - 点の重み(0wi1090 \leq w_i 10^9).

# 出力

C

# 入出力例

4
2 3
3 1
6 1
0 2
3

# 解説

NN2×1052 \times 10^5と大きいので計算量的にO(NlogN)O(N \log{N})以下でないと間に合いません.
これはソート1回と線形探索1回が許される程度の計算量です.

問題文中でも言及されている通り、一般的な最大クリーク問題を多項式時間で解くのは難しそうです.
全探索で2N2^Nで答えが求まりますがここではそれは使えません.

そこでこの問題における"重み"の意味を考えます.
よくあるグラフ問題の重みとは違って実はもっとシンプルなものであることが分かります.
例えば座標1010で重み55の点であれば10±510 \pm 5の位置の点とは繋がらないということです.
つまりここでの"重み"とは長さのことです.
105=510 - 5 = 5を始点、10+5=1510 + 5 = 15を終点とする線分として捉え直します.
"以上"なので始点と終端どちらかは含めませんので、ここでは終点を含まず[5:15)[5:15)とします.
(整数でいうと[5:14][5:14]です.)

それを踏まえて上の入力例を見てみましょう.

まず1つ目の点です.

2 3

座標がxで重みが3です.
線分として捉えると始点-1、終点4の線分([1:5)[-1:5))です.
範囲0x1090 \leq x \leq 10^9を超える場合はカットします.
始点0、終点4となります.

座標をo、それ以外の線分の部分を-で表してプロットすると以下のようになります.
0から7の番号は座標です.

   0 1 2 3 4 5 6 7
1: - - o - -

次に2つ目の点です.

3 1

線分として捉え直すと、始点2、終点3です.
プロットすると、

   0 1 2 3 4 5 6 7
1: - - o - -
2:     - o

3つ目、4つ目も同様に線分として捉える直すとそれぞれ[5:6][5:6][0:1][0:1]です.
プロットすると、

   0 1 2 3 4 5 6 7
1: - - o - -
2:     - o
3:           - o
4: o -

線分が重なっていないものは直接繋げられるということなので、答えは1以外の3点となります.

これを許容計算量であるソート+ 1回の走査で求めるには、まずすべての線分を終点が小さい順にソートします.
これは一番終端が小さい線分ほど他の線分と重なる可能性が低いからです.
より多くの点と繋がることが出来る可能性が高いからです.

上記の例ではこのようになります.

   0 1 2 3 4 5 6 7
4: o -
2:     - o
1: - - o - -
3:           - o

あとは先頭から見ていって最後に追加した線分と重なっていたら捨てる.
重なっていなかったら線分として追加する、ということを繰り返すだけです.

この走査が終わったら追加した線分の数が答えになります.

# 計算量

貪欲法のループはO(N)O(N)ですが、終端点による線分のソートにO(NlogN)O(N \log{N})かかります.

O(NlogN)O(N \log{N})

# 解答

// C++ 14
#include <iostream>
#include <string>
#include <vector>
#include <list>
#include <algorithm>
#include <queue>
#include <stack>
#include <set>
#include <map>
#include <unordered_map>
#include <math.h>

#define ll long long
#define Int int
#define loop(x, start, end) for(Int x = start; x < end; x++)
#define loopdown(x, start, end) for(int x = start; x > end; x--)
#define rep(n) for(int x = 0; x < n; x++)
#define span(a,x,y) a.begin()+x,a.begin()+y
#define span_all(a) a.begin(),a.end()
#define len(x) (x.size())
#define last(x) (*(x.end()-1))

using namespace std;
#define MAX_N 200000
#define MAX_X 1000000000
Int N;
vector<Segment> segments;

void input() {
  cin >> N;
  Int x_, w_;
  loop(n,0,N) {
    cin >> x_ >> w_;
    Int a = x_ - w_, b = x_ + w_ - 1; // [a:b)
    if (a < 0) a = 0;
    if (b > MAX_X) b = MAX_X;
    segments.push_back(Segment(a, b));
  }
}

void solve() {
  sort(span_all(segments));
  Int count = 0, right = -1;
  for (auto s: segments) {
    if (right < s.start) {
      count++;
      right = s.end;
    }
  }
  cout << count << endl;
}

int main() {
  input();
  solve();
  return 0;
}

リモートフリーランス。ウェブサービス、スマホアプリエンジニア。
東アジアを拠点に世界を移動しながら活動してます!

お仕事のご依頼・お問い合わせはこちら

コメント