Submission #119384

# Submission time Handle Problem Language Result Execution time Memory
119384 2019-06-21T06:42:08 Z tmwilliamlin168 Printed Circuit Board (CEOI12_circuit) C++14
95 / 100
100 ms 19808 KB
#include <bits/stdc++.h>
using namespace std;

int in() {
	char c=0;
	while(c<'0'||c>'9')
		c=getchar_unlocked();
	int r=0;
	while(c>='0'&&c<='9') {
		r=r*10+c-'0';
		c=getchar_unlocked();
	}
	return r;
}

void out(int x) {
	int c=0, r=0;
	for(; x%10==0; ++c, x/=10);
	for(; x; x/=10)
		r=r*10+x%10;
	for(; r; r/=10)
		putchar_unlocked(r%10+'0');
	while(c--)
		putchar_unlocked('0');
}

const int mxN=2e5, mxA=1e6;
int n, b[2][mxN], c[mxA+1], d[mxN+1], e1[mxN], e2[mxN], a1, t[2][mxN], tp[mxN];
long long x[mxN+1], y[mxN+1];
array<int, 3> a[mxN];
vector<array<int, 2>> r[18];
bool ans[mxN];

bool cw(int a, int b, int c) {
	return (y[c]-y[a])*(x[b]-x[a])<(y[b]-y[a])*(x[c]-x[a]);
}

bool scmp(const int &i, const int &j) {
	if(i<0)
		return 0;
	if(j<0)
		return 1;
	if(e1[i]==e1[j])
		return d[e2[i]]<d[e2[j]]?!cw(e1[j], e2[j], e2[i]):cw(e1[i], e2[i], e2[j]);
	return d[e1[i]]<d[e1[j]]?cw(e1[i], e2[i], e1[j]):!cw(e1[j], e2[j], e1[i]);
}

int main() {
	n=in();
	for(int i=0; i<n; ++i) {
		x[i]=in(), y[i]=in();
		long long z1=y[i]*mxA*mxA/x[i]-mxA, z2=z1/mxA, z3=z2/mxA;
		a[i]={z1-z2*mxA, z2-z3*mxA, z3};
	}
	x[n]=x[0];
	y[n]=y[0];
	iota(b[0], b[0]+n, 0);
	for(int k=0; k<3; ++k) {
		memset(c+1, 0, 4*mxA);
		for(int i=0; i<n; ++i)
			++c[a[i][k]+1];
		for(int i=0; i<mxA; ++i)
			c[i+1]+=c[i];
		for(int i=n-1; ~i; --i)
			b[k&1^1][--c[a[b[k&1][i]][k]+1]]=b[k&1][i];
	}
	memset(tp, -1, 4*n);
	for(int i=0, j=0; i<n; i=j) {
		for(; j<n&&a[b[1][j]]==a[b[1][i]]; ++j) {
			d[b[1][j]]=i;
			if(tp[i]<0||x[b[1][j]]+y[b[1][j]]<x[tp[i]]+y[tp[i]])
				tp[i]=b[1][j];
		}
	}
	d[n]=d[0];
	for(int i=0; i<n; ++i) {
		e1[i]=i;
		e2[i]=i+1;
		if(d[e2[i]]<d[e1[i]])
			swap(e1[i], e2[i]);
		if(d[e1[i]]<d[e2[i]]) {
			int k=31-__builtin_clz(d[e2[i]]-d[e1[i]]);
			r[k].push_back({d[e1[i]], i});
			r[k].push_back({d[e2[i]]-(1<<k), i});
		}
	}
	memset(t[1], -1, 4*n);
	for(int i=17; ~i; --i) {
		for(array<int, 2> a : r[i])
			t[i&1][a[0]]=min(a[1], t[i&1][a[0]], scmp);
		if(i) {
			memset(t[i&1^1], -1, 4*n);
			for(int j=0; j<=n-(1<<i); ++j) {
				t[i&1^1][j]=min(t[i&1][j], t[i&1^1][j], scmp);
				t[i&1^1][j+(1<<i-1)]=min(t[i&1][j], t[i&1^1][j+(1<<i-1)], scmp);
			}
		}
	}
	for(int i=0; i<n; ++i)
		if(~tp[i]&&(t[0][i]<0||!cw(e1[t[0][i]], e2[t[0][i]], tp[i])))
			ans[tp[i]]=1;
	for(int i=0; i<n; ++i)
		a1+=ans[i];
	out(a1);
	putchar_unlocked('\n');
	for(int i=0; i<n; ++i) {
		if(ans[i]) {
			out(i+1);
			putchar_unlocked(' ');
		}
	}
}

Compilation message

circuit.cpp: In function 'int main()':
circuit.cpp:53:11: warning: narrowing conversion of '(z1 - (z2 * 1000000))' from 'long long int' to 'int' inside { } [-Wnarrowing]
   a[i]={z1-z2*mxA, z2-z3*mxA, z3};
         ~~^~~~~~~
circuit.cpp:53:22: warning: narrowing conversion of '(z2 - (z3 * 1000000))' from 'long long int' to 'int' inside { } [-Wnarrowing]
   a[i]={z1-z2*mxA, z2-z3*mxA, z3};
                    ~~^~~~~~~
circuit.cpp:53:33: warning: narrowing conversion of 'z3' from 'long long int' to 'int' inside { } [-Wnarrowing]
   a[i]={z1-z2*mxA, z2-z3*mxA, z3};
                                 ^
circuit.cpp:65:7: warning: suggest parentheses around arithmetic in operand of '^' [-Wparentheses]
    b[k&1^1][--c[a[b[k&1][i]][k]+1]]=b[k&1][i];
      ~^~
circuit.cpp:92:14: warning: suggest parentheses around arithmetic in operand of '^' [-Wparentheses]
    memset(t[i&1^1], -1, 4*n);
             ~^~
circuit.cpp:94:8: warning: suggest parentheses around arithmetic in operand of '^' [-Wparentheses]
     t[i&1^1][j]=min(t[i&1][j], t[i&1^1][j], scmp);
       ~^~
circuit.cpp:94:35: warning: suggest parentheses around arithmetic in operand of '^' [-Wparentheses]
     t[i&1^1][j]=min(t[i&1][j], t[i&1^1][j], scmp);
                                  ~^~
circuit.cpp:95:8: warning: suggest parentheses around arithmetic in operand of '^' [-Wparentheses]
     t[i&1^1][j+(1<<i-1)]=min(t[i&1][j], t[i&1^1][j+(1<<i-1)], scmp);
       ~^~
circuit.cpp:95:21: warning: suggest parentheses around '-' inside '<<' [-Wparentheses]
     t[i&1^1][j+(1<<i-1)]=min(t[i&1][j], t[i&1^1][j+(1<<i-1)], scmp);
                    ~^~
circuit.cpp:95:44: warning: suggest parentheses around arithmetic in operand of '^' [-Wparentheses]
     t[i&1^1][j+(1<<i-1)]=min(t[i&1][j], t[i&1^1][j+(1<<i-1)], scmp);
                                           ~^~
circuit.cpp:95:57: warning: suggest parentheses around '-' inside '<<' [-Wparentheses]
     t[i&1^1][j+(1<<i-1)]=min(t[i&1][j], t[i&1^1][j+(1<<i-1)], scmp);
                                                        ~^~
# Verdict Execution time Memory Grader output
1 Correct 7 ms 4352 KB Output is correct
2 Correct 7 ms 4480 KB Output is correct
3 Correct 8 ms 4608 KB Output is correct
4 Correct 8 ms 4608 KB Output is correct
5 Correct 14 ms 5120 KB Output is correct
6 Correct 11 ms 5248 KB Output is correct
7 Correct 14 ms 6016 KB Output is correct
8 Correct 10 ms 5120 KB Output is correct
9 Correct 10 ms 5120 KB Output is correct
10 Correct 12 ms 5120 KB Output is correct
11 Correct 12 ms 5248 KB Output is correct
12 Correct 13 ms 5988 KB Output is correct
13 Correct 19 ms 6772 KB Output is correct
14 Correct 21 ms 7552 KB Output is correct
15 Correct 23 ms 8308 KB Output is correct
16 Correct 41 ms 12136 KB Output is correct
17 Correct 58 ms 12136 KB Output is correct
18 Correct 78 ms 19736 KB Output is correct
19 Correct 78 ms 19808 KB Output is correct
20 Execution timed out 128 ms 19804 KB Time limit exceeded