Submission #429666

# Submission time Handle Problem Language Result Execution time Memory
429666 2021-06-16T08:35:34 Z 송준혁(#7525) Power plants (CPSPC17_power) C++17
35 / 100
6000 ms 125592 KB
#include <bits/stdc++.h>
#define fi first
#define se second
#define pb push_back
#define lb lower_bound
#define MOD 1000000007
#define HASH 2979659
#define INF (1ll<<61)
using namespace std;
typedef long long LL;
typedef pair<int,int> pii;
typedef pair<LL,int> pli;

int N, M, Q, num;
int chk[101010];
pii B[101010];
LL ans, X[101010], Y[101010];
vector<int> V[101010], ansX, ansY;
vector<pli> H[3030303];

int bid(LL x, LL k){
	int lo=0, hi=MOD, ret=0;
	while (lo<=hi){
		LL m=lo+hi>>1;
		if (m*m <= 2*x*x/k) ret=m, lo=m+1;
		else hi=m-1;
	}
	return ret;
}

inline int Hash(int x, int y){return (((LL)(x+5)*MOD + (y+5)) ^ 19059) % HASH;}

int id(int x, int y){
	int t = Hash(x, y);
	for (pli it : H[t]) if (it.fi == (LL)x*MOD+y) return it.se;
	return -1;
}

bool dfs(int u, int c, LL k){
	if (chk[u]) return chk[u]*c < 0;
	chk[u] = c;
	for (int dx=-2; dx<=2; dx++) for (int dy=-2; dy<=2; dy++){
		if (dx*dy == 4 || dx*dy == -4) continue;
		int t = id(B[u].fi+dx, B[u].se+dy);
		if (t < 0) continue;
		for (int v : V[t]) if (v != u && (X[u]-X[v])*(X[u]-X[v])+(Y[u]-Y[v])*(Y[u]-Y[v]) <= k){
			if (dfs(v, -c, k)) return true;
		}
	}
	return false;
}

bool f(LL k){
	for (int i=0; i<num; i++) V[i].clear();
	for (int i=1; i<=N; i++) H[Hash(B[i].fi, B[i].se)].clear();
	for (int i=1; i<=N; i++) chk[i] = 0;
	num = 0;
	for (int i=1; i<=N; i++){
		B[i].fi = bid(X[i], k);
		B[i].se = bid(Y[i], k);
		if (id(B[i].fi, B[i].se) < 0){
			H[Hash(B[i].fi, B[i].se)].pb(pli((LL)MOD*B[i].fi + B[i].se, num));
			num++;
		}
	}
	for (int i=1; i<=N; i++){
		int t = id(B[i].fi, B[i].se);
		V[t].pb(i);
		if (V[t].size() > 2) return false;
	}
	for (int i=1; i<=N; i++){
		if (chk[i]) continue;
		if (dfs(i, 1, k)) return false;
	}
	return true;
}

int main(){
	scanf("%d", &N);
	for (int i=1; i<=N; i++) scanf("%lld %lld", &X[i], &Y[i]);
	LL lo=0, hi=INF;
	while (lo<=hi){
		LL m = lo+hi>>1;
		if (f(m)) ans = m, lo = m + 1;
		else hi = m - 1;
	}
	printf("%lld\n", ans+1);
	f(ans);
	for (int i=1; i<=N; i++){
		if (chk[i] > 0) ansX.pb(i);
		else ansY.pb(i);
	}
	printf("%d\n", ansX.size());
	for (int x : ansX) printf("%d ", x);
	printf("\n");
	printf("%d\n", ansY.size());
	for (int x : ansY) printf("%d ", x);
	printf("\n");
	return 0;
}

Compilation message

Main.cpp: In function 'int bid(LL, LL)':
Main.cpp:24:10: warning: suggest parentheses around '+' inside '>>' [-Wparentheses]
   24 |   LL m=lo+hi>>1;
      |        ~~^~~
Main.cpp: In function 'int main()':
Main.cpp:83:12: warning: suggest parentheses around '+' inside '>>' [-Wparentheses]
   83 |   LL m = lo+hi>>1;
      |          ~~^~~
Main.cpp:93:11: warning: format '%d' expects argument of type 'int', but argument 2 has type 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wformat=]
   93 |  printf("%d\n", ansX.size());
      |          ~^     ~~~~~~~~~~~
      |           |              |
      |           int            std::vector<int>::size_type {aka long unsigned int}
      |          %ld
Main.cpp:96:11: warning: format '%d' expects argument of type 'int', but argument 2 has type 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wformat=]
   96 |  printf("%d\n", ansY.size());
      |          ~^     ~~~~~~~~~~~
      |           |              |
      |           int            std::vector<int>::size_type {aka long unsigned int}
      |          %ld
Main.cpp:79:7: warning: ignoring return value of 'int scanf(const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
   79 |  scanf("%d", &N);
      |  ~~~~~^~~~~~~~~~
Main.cpp:80:32: warning: ignoring return value of 'int scanf(const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
   80 |  for (int i=1; i<=N; i++) scanf("%lld %lld", &X[i], &Y[i]);
      |                           ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~
# Verdict Execution time Memory Grader output
1 Correct 46 ms 73868 KB Output is correct
2 Correct 46 ms 73828 KB Output is correct
3 Correct 44 ms 73788 KB Output is correct
4 Correct 44 ms 73872 KB Output is correct
5 Correct 47 ms 73772 KB Output is correct
6 Correct 45 ms 73776 KB Output is correct
7 Correct 44 ms 73748 KB Output is correct
8 Correct 58 ms 73744 KB Output is correct
9 Correct 44 ms 73800 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 46 ms 73868 KB Output is correct
2 Correct 46 ms 73828 KB Output is correct
3 Correct 44 ms 73788 KB Output is correct
4 Correct 44 ms 73872 KB Output is correct
5 Correct 47 ms 73772 KB Output is correct
6 Correct 45 ms 73776 KB Output is correct
7 Correct 44 ms 73748 KB Output is correct
8 Correct 58 ms 73744 KB Output is correct
9 Correct 44 ms 73800 KB Output is correct
10 Correct 90 ms 74808 KB Output is correct
11 Correct 76 ms 74604 KB Output is correct
12 Correct 72 ms 74564 KB Output is correct
13 Correct 77 ms 74884 KB Output is correct
14 Correct 79 ms 74840 KB Output is correct
15 Correct 63 ms 74232 KB Output is correct
16 Correct 87 ms 74524 KB Output is correct
17 Correct 87 ms 74644 KB Output is correct
18 Correct 79 ms 74772 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 46 ms 73868 KB Output is correct
2 Correct 46 ms 73828 KB Output is correct
3 Correct 44 ms 73788 KB Output is correct
4 Correct 44 ms 73872 KB Output is correct
5 Correct 47 ms 73772 KB Output is correct
6 Correct 45 ms 73776 KB Output is correct
7 Correct 44 ms 73748 KB Output is correct
8 Correct 58 ms 73744 KB Output is correct
9 Correct 44 ms 73800 KB Output is correct
10 Correct 90 ms 74808 KB Output is correct
11 Correct 76 ms 74604 KB Output is correct
12 Correct 72 ms 74564 KB Output is correct
13 Correct 77 ms 74884 KB Output is correct
14 Correct 79 ms 74840 KB Output is correct
15 Correct 63 ms 74232 KB Output is correct
16 Correct 87 ms 74524 KB Output is correct
17 Correct 87 ms 74644 KB Output is correct
18 Correct 79 ms 74772 KB Output is correct
19 Execution timed out 6059 ms 125592 KB Time limit exceeded
20 Halted 0 ms 0 KB -