Submission #429659

# Submission time Handle Problem Language Result Execution time Memory
429659 2021-06-16T08:30:00 Z 송준혁(#7525) Power plants (CPSPC17_power) C++17
35 / 100
6000 ms 127856 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 <= 4*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++){
		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:82:12: warning: suggest parentheses around '+' inside '>>' [-Wparentheses]
   82 |   LL m = lo+hi>>1;
      |          ~~^~~
Main.cpp:92:11: warning: format '%d' expects argument of type 'int', but argument 2 has type 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wformat=]
   92 |  printf("%d\n", ansX.size());
      |          ~^     ~~~~~~~~~~~
      |           |              |
      |           int            std::vector<int>::size_type {aka long unsigned int}
      |          %ld
Main.cpp:95:11: warning: format '%d' expects argument of type 'int', but argument 2 has type 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wformat=]
   95 |  printf("%d\n", ansY.size());
      |          ~^     ~~~~~~~~~~~
      |           |              |
      |           int            std::vector<int>::size_type {aka long unsigned int}
      |          %ld
Main.cpp:78:7: warning: ignoring return value of 'int scanf(const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
   78 |  scanf("%d", &N);
      |  ~~~~~^~~~~~~~~~
Main.cpp:79:32: warning: ignoring return value of 'int scanf(const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
   79 |  for (int i=1; i<=N; i++) scanf("%lld %lld", &X[i], &Y[i]);
      |                           ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~
# Verdict Execution time Memory Grader output
1 Correct 45 ms 73864 KB Output is correct
2 Correct 50 ms 73756 KB Output is correct
3 Correct 45 ms 73796 KB Output is correct
4 Correct 51 ms 73796 KB Output is correct
5 Correct 45 ms 73840 KB Output is correct
6 Correct 47 ms 73772 KB Output is correct
7 Correct 44 ms 73796 KB Output is correct
8 Correct 44 ms 73844 KB Output is correct
9 Correct 49 ms 73808 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 45 ms 73864 KB Output is correct
2 Correct 50 ms 73756 KB Output is correct
3 Correct 45 ms 73796 KB Output is correct
4 Correct 51 ms 73796 KB Output is correct
5 Correct 45 ms 73840 KB Output is correct
6 Correct 47 ms 73772 KB Output is correct
7 Correct 44 ms 73796 KB Output is correct
8 Correct 44 ms 73844 KB Output is correct
9 Correct 49 ms 73808 KB Output is correct
10 Correct 86 ms 74984 KB Output is correct
11 Correct 82 ms 74688 KB Output is correct
12 Correct 76 ms 74588 KB Output is correct
13 Correct 76 ms 74908 KB Output is correct
14 Correct 83 ms 74948 KB Output is correct
15 Correct 64 ms 74220 KB Output is correct
16 Correct 81 ms 74636 KB Output is correct
17 Correct 104 ms 74772 KB Output is correct
18 Correct 93 ms 74844 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 45 ms 73864 KB Output is correct
2 Correct 50 ms 73756 KB Output is correct
3 Correct 45 ms 73796 KB Output is correct
4 Correct 51 ms 73796 KB Output is correct
5 Correct 45 ms 73840 KB Output is correct
6 Correct 47 ms 73772 KB Output is correct
7 Correct 44 ms 73796 KB Output is correct
8 Correct 44 ms 73844 KB Output is correct
9 Correct 49 ms 73808 KB Output is correct
10 Correct 86 ms 74984 KB Output is correct
11 Correct 82 ms 74688 KB Output is correct
12 Correct 76 ms 74588 KB Output is correct
13 Correct 76 ms 74908 KB Output is correct
14 Correct 83 ms 74948 KB Output is correct
15 Correct 64 ms 74220 KB Output is correct
16 Correct 81 ms 74636 KB Output is correct
17 Correct 104 ms 74772 KB Output is correct
18 Correct 93 ms 74844 KB Output is correct
19 Execution timed out 6018 ms 127856 KB Time limit exceeded
20 Halted 0 ms 0 KB -