Submission #852214

# Submission time Handle Problem Language Result Execution time Memory
852214 2023-09-21T12:33:27 Z MilosMilutinovic Tri (CEOI09_tri) C++14
90 / 100
2000 ms 13928 KB
#include <bits/stdc++.h>
using namespace std;
#define rep(i,a,n) for (int i=a;i<n;i++)
#define per(i,a,n) for (int i=n-1;i>=a;i--)
#define pb push_back
#define eb emplace_back
#define mp make_pair
#define all(x) (x).begin(),(x).end()
#define fi first
#define se second
#define SZ(x) ((int)(x).size())
typedef vector<int> VI;
typedef basic_string<int> B;
typedef long long ll;
typedef pair<int,int> PII;
typedef double db;
mt19937 mrand(random_device{}()); 
const ll mod=1000000007;
int rnd(int x) { return mrand() % x;}
ll powmod(ll a,ll b) {ll res=1;a%=mod; assert(b>=0); for(;b;b>>=1){if(b&1)res=res*a%mod;a=a*a%mod;}return res;}
ll gcd(ll a,ll b) { return b?gcd(b,a%b):a;}
// head

const int N=100100;
const int L=20;
const int S=sqrt(N);
int n,m,st[N][L],logs[N];
ll d[N];

mt19937 rng(time(0));

struct pt {
	int x,y;
	pt(int _x=0,int _y=0) : x(_x),y(_y) {}
}p[N];
vector<pt> ch[N];
int orient(pt a,pt b,pt c) {
	ll v=(b.y-a.y)*1ll*(c.x-b.x)-(b.x-a.x)*1ll*(c.y-b.y);
	return v==0?0:(v<0?-1:1);
}
bool const operator<(const pt& a,const pt& b) {
	if (orient(pt(0,0),a,b)==0) {
		return (a.x==b.x?a.y<b.y:a.x<b.x);
	}
	return orient(pt(0,0),a,b)<0;
}
bool in_tri(pt a,pt b,pt c,pt p) {
	return orient(a,b,p)!=1&&orient(b,c,p)!=1&&orient(c,a,p)!=1;
}
int cmp(int i,int j) {
	return d[i]<d[j]?i:j;
}
int query(int l,int r) {
	int k=logs[r-l+1];
	return cmp(st[l][k],st[r-(1<<k)+1][k]);
}
int main() {
	scanf("%d%d",&n,&m);
	rep(i,0,n) {
		scanf("%d%d",&p[i].x,&p[i].y);
	}
	sort(p,p+n);
	rep(i,0,n) {
		st[i][0]=i;
		d[i]=p[i].x*1ll*p[i].x+p[i].y*1ll*p[i].y;
	}
	rep(i,2,n+1) logs[i]=logs[i>>1]+1;
	rep(j,1,L) {
		rep(i,0,n-(1<<j)+1) {
			st[i][j]=cmp(st[i][j-1],st[i+(1<<(j-1))][j-1]);
		}
	}
	rep(i,0,n) {
		ch[i/S].pb(p[i]);
	}
	rep(i,0,n+1) {
		if (ch[i].empty()) continue;
		vector<pt> hull;
		sort(all(ch[i]),[&](pt a,pt b) { return a.x<b.x||(a.x==b.x&&a.y<b.y); });
		for (auto& p:ch[i]) {
			while (SZ(hull)>=2&&orient(hull[SZ(hull)-2],hull[SZ(hull)-1],p)>0) {
				hull.pop_back();
			}
			hull.pb(p);
		}
		ch[i]=hull;
	}
	while (m--) {
		pt a,b;
		scanf("%d%d%d%d",&a.x,&a.y,&b.x,&b.y);
		if (orient(b,pt(0,0),a)>=0) swap(a,b);
		int L,R;
		{
			L=n-1;
			int low=0,high=n-1;
			while (low<=high) {
				int mid=low+high>>1;
				if (orient(p[mid],pt(0,0),a)<=0) {
					L=mid;
					high=mid-1;
				} else {
					low=mid+1;
				}
			}
		}
		{
			R=0;
			int low=0,high=n-1;
			while (low<=high) {
				int mid=low+high>>1;
				if (orient(p[mid],pt(0,0),b)>=0) {
					R=mid;
					low=mid+1;
				} else {
					high=mid-1;
				}
			}
		}
		if (L>R) {
			printf("N\n");
			continue;
		}
		if (R-L<=S) {
			bool ok=false;
			rep(i,L,R+1) if (in_tri(a,b,pt(0,0),p[i])) ok=true;
			puts(ok?"Y":"N");
			continue;
		}
		bool ok=false;
		while (L%S!=0) {
			if (in_tri(a,b,pt(0,0),p[L])) ok=true;
			L++;
		}
		while (R>=L&&R%S!=S-1) {
			if (in_tri(a,b,pt(0,0),p[R])) ok=true;
			R--;
		}
		VI ord;
		rep(i,L/S,R/S+1) ord.pb(i);
		shuffle(ord.begin(),ord.end(),rng);
		for (int i:ord) {
			for (auto& p:ch[i]) {
				if (in_tri(a,b,pt(0,0),p)) ok=true;
			}
			if (ok) break;
		}
		puts(ok?"Y":"N");
	}
}

Compilation message

tri.cpp: In function 'int main()':
tri.cpp:97:16: warning: suggest parentheses around '+' inside '>>' [-Wparentheses]
   97 |     int mid=low+high>>1;
      |             ~~~^~~~~
tri.cpp:110:16: warning: suggest parentheses around '+' inside '>>' [-Wparentheses]
  110 |     int mid=low+high>>1;
      |             ~~~^~~~~
tri.cpp:58:7: warning: ignoring return value of 'int scanf(const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
   58 |  scanf("%d%d",&n,&m);
      |  ~~~~~^~~~~~~~~~~~~~
tri.cpp:60:8: warning: ignoring return value of 'int scanf(const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
   60 |   scanf("%d%d",&p[i].x,&p[i].y);
      |   ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~
tri.cpp:90:8: warning: ignoring return value of 'int scanf(const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
   90 |   scanf("%d%d%d%d",&a.x,&a.y,&b.x,&b.y);
      |   ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Verdict Execution time Memory Grader output
1 Correct 2 ms 6744 KB Output is correct
2 Correct 3 ms 6748 KB Output is correct
3 Correct 123 ms 9256 KB Output is correct
4 Correct 1347 ms 9552 KB Output is correct
5 Execution timed out 2035 ms 13928 KB Time limit exceeded
6 Correct 651 ms 12024 KB Output is correct
7 Correct 837 ms 13904 KB Output is correct
8 Correct 187 ms 12072 KB Output is correct
9 Correct 201 ms 13028 KB Output is correct
10 Correct 234 ms 13924 KB Output is correct