QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#373100 | #4371. Spin Doctor | InfinityNS | WA | 81ms | 13812kb | C++14 | 4.9kb | 2024-04-01 02:25:17 | 2024-04-01 02:25:19 |
Judging History
answer
#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define pb push_back
#define mp make_pair
struct pt{
ll x,y;
pt():x(0),y(0){}
pt(ll a,ll b):x(a),y(b){}
};
pt operator - (pt a,pt b){return pt(a.x-b.x,a.y-b.y);}
pt operator - (pt a){return pt(-a.x,-a.y);}
ll cross(pt a,pt b){return a.x*b.y-a.y*b.x;}
ll dot(pt a,pt b){return a.x*b.x+a.y*b.y;}
int part(pt a){return a.y<0 || (a.y==0 && a.x<0);}
bool operator < (pt a,pt b){return mp(part(a),(ll)0)<mp(part(b),cross(a,b));}
bool operator == (pt a,pt b){return mp(a.x,a.y)==mp(b.x,b.y);}
int sgn(ll x){return x==0?0:(x<0?-1:1);}
const int N=250050;
pt pts[N];
int c[N],n;
bool cmp(pt a,pt b){return mp(a.x,a.y)<mp(b.x,b.y);}
vector<pt> ConvexHull(vector<pt> poly){
sort(poly.begin(),poly.end(),cmp);
poly.erase(unique(poly.begin(),poly.end()),poly.end());
if(poly.size()<2)return poly;
vector<pt> hull;
int sz=0;
for(int t=0;t<2;t++){
int was=sz;
for(auto p:poly){
while(sz>was+1 && cross(hull[sz-1]-hull[sz-2],p-hull[sz-2])<=0){
hull.pop_back();
sz--;
}
hull.pb(p);
sz++;
}
hull.pop_back();
sz--;
reverse(poly.begin(),poly.end());
}
return hull;
}
bool OnSeg(pt a,pt b,pt p){
pt v=b-a;
if(cross(v,p-a)!=0)return false;
ll da=dot(v,a);
ll dp=dot(v,p);
ll db=dot(v,b);
return da<=dp && dp<=db;
}
bool InTriangle(pt a,pt b,pt c,pt p){
vector<pt> poly={a,b,c};
bool lo=false,hi=false;
for(int i=0;i<3;i++){
ll cr=cross(poly[i]-p,poly[(i+1)%3]-p);
if(cr<0)lo=true;
if(cr>0)hi=true;
}
return !lo || !hi;
}
bool Inside(vector<pt>& poly, pt p){
if(poly.size()==1)return p==poly[0];
if(poly.size()==2)return OnSeg(poly[0],poly[1],p);
int bot=0,top=(int)poly.size()-3,ans=0;
while(bot<=top){
int mid=top+bot>>1;
if(cross(poly[mid]-poly.back(),p-poly.back())>=0){
bot=mid+1;
ans=mid;
}else{
top=mid-1;
}
}
//printf("Inside? %i\n",ans);
return InTriangle(poly[ans],poly[ans+1],poly.back(),p);
}
pt Solve(vector<pt>& poly,pt p,int l,int r,int sg){
/*int bot=l,top=r-1,ans=r;
while(bot<=top){
int mid=bot+top>>1;
if(sgn(cross(poly[mid]-p,poly[mid+1]-p))==sg){
bot=mid+1;
}else{
ans=mid;
top=mid-1;
}
}
//printf("Solve (%lld, %lld) l:%i r:%i -- %i (%lld, %lld)\n",p.x,p.y,l,r,ans,poly[ans].x,poly[ans].y);
return poly[ans]-p;*/
int ans=l;
while(ans<r && sgn(cross(poly[ans]-p,poly[ans+1]-p))==sg)ans++;
return poly[ans]-p;
}
pair<pt,pt> Tangents(vector<pt>& poly, pt p){
int l=0,r=0;
for(int i=0;i<poly.size();i++){
if(sgn(cross(poly[l]-p,poly[i]-p))==1){
l=i;
}
if(sgn(cross(poly[r]-p,poly[i]-p))==-1){
r=i;
}
}
return {poly[l]-p,poly[r]-p};
if(poly.size()==1)return {poly[0]-p,poly[0]-p};
if(poly.size()==2){
if(cross(poly[0]-p,poly[1]-p)<=0)return {poly[0]-p,poly[1]-p};
else return {poly[1]-p,poly[0]-p};
}
int bot=1,top=(int)poly.size()-2,ans=0;
int s1=sgn(cross(poly.back()-p,poly.front()-p));
while(top>=bot){
int mid=top+bot>>1;
if(sgn(cross(poly.back()-p,poly[mid]-p))==s1){
ans=mid;
bot=mid+1;
}else{
top=mid-1;
}
}
//printf("Split\n");
//for(int i=0;i<=ans;i++)printf("(%lld, %lld) ",poly[i].x,poly[i].y);printf("\n");
//for(int i=ans+1;i<poly.size();i++)printf("(%lld, %lld) ",poly[i].x,poly[i].y);printf("\n");
//if(n==100 && s1==0)printf("s1 before %i\n",s1);
if(s1==0)s1=-sgn(cross(poly.back()-p,poly[ans+1]-p));
pt L=Solve(poly,p,0,ans,s1);
pt R=Solve(poly,p,ans+1,(int)poly.size()-1,-s1);
if(s1==-1)swap(L,R);
return {L,R};
}
int main(){
scanf("%i",&n);
vector<pt> poly;
for(int i=1;i<=n;i++){
scanf("%lld %lld %i",&pts[i].x,&pts[i].y,&c[i]);
if(c[i]==1)poly.pb(pts[i]);
}
if(poly.size()==1){
printf("1\n");
return 0;
}
poly=ConvexHull(poly);
//if(n==100)printf("CH size: %i\n",poly.size());
//printf("Convex hull: ");
//for(pt p:poly)printf("(%lld, %lld) ",p.x,p.y);
//printf("\n");
ll mny=poly[0].y,mxy=poly[0].y;
for(int i=1;i<poly.size();i++){
mny=min(mny,poly[i].y);
mxy=max(mxy,poly[i].y);
}
int cnt=0;
vector<pt> cand;
for(int i=1;i<=n;i++){
if(Inside(poly,pts[i])){
//printf("%i(%lld, %lld) is inside\n",i,pts[i].x,pts[i].y);
cnt++;
}else{
cand.pb(pts[i]);
}
}
//if(n==100)printf("Inside %i\n",cnt);
vector<pair<pt,int>> evs;
for(int i=0;i<cand.size();i++){
//printf("Cand (%lld, %lld)\n",cand[i].x,cand[i].y);
pt L,R;
tie(L,R)=Tangents(poly,cand[i]);
//printf("Tangents L(%lld, %lld) R(%lld, %lld)\n",L.x,L.y,R.x,R.y);
if(part(L)==1)L=-L;
if(part(R)==1)R=-R;
pt S=R;
pt E=L;
evs.pb({S,-1});
evs.pb({E,1});
//printf("(%lld, %lld) enter(%lld, %lld) exit(%lld, %lld)\n",cand[i].x,cand[i].y,S.x,S.y,E.x,E.y);
//if(L<R)cnt++;
if(cand[i].y>=mny && cand[i].y<=mxy)cnt++;
}
int ans=cnt;
sort(evs.begin(),evs.end());
for(auto ev:evs){
cnt-=ev.second;
ans=min(ans,cnt);
}
printf("%i\n",ans);
return 0;
}
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 100
Accepted
time: 1ms
memory: 8572kb
input:
6 0 10 0 10 0 1 12 8 1 5 5 0 11 2 1 11 3 0
output:
4
result:
ok single line: '4'
Test #2:
score: 0
Accepted
time: 1ms
memory: 8436kb
input:
10 6 1 1 0 2 0 2 1 1 6 1 1 8 2 0 4 4 0 4 0 0 2 3 1 6 1 0 6 3 1
output:
8
result:
ok single line: '8'
Test #3:
score: 0
Accepted
time: 1ms
memory: 7644kb
input:
5 5 7 0 3 4 0 5 7 0 5 7 1 9 4 0
output:
1
result:
ok single line: '1'
Test #4:
score: 0
Accepted
time: 2ms
memory: 7952kb
input:
100 7487 4751 1 7499 5064 1 7471 5376 1 7404 5683 1 7300 5979 1 7159 6260 1 6984 6520 1 6777 6757 1 6543 6966 1 6284 7144 1 6006 7288 1 5711 7396 1 5405 7466 1 5092 7498 1 4780 7490 1 4469 7442 1 4167 7357 1 3878 7234 1 3607 7075 1 3358 6884 1 3135 6664 1 2941 6417 1 2780 6147 1 2653 5860 1 2564 555...
output:
61
result:
ok single line: '61'
Test #5:
score: 0
Accepted
time: 2ms
memory: 7716kb
input:
200 7489 1285 0 6851 8471 0 6122 2766 1 2413 9338 0 1725 7382 0 6984 6520 1 5080 8417 0 2604 5711 1 5833 2643 1 48 2810 0 5316 2439 0 900 7419 0 6809 867 0 6006 7288 1 5092 7498 1 5531 2558 1 7075 6393 1 9979 9313 0 7436 4441 1 4595 2534 1 2598 909 0 842 284 0 3358 6884 1 6522 7976 0 604 5833 0 3607...
output:
125
result:
ok single line: '125'
Test #6:
score: 0
Accepted
time: 2ms
memory: 7712kb
input:
300 2253 823 1 1865 9556 0 306 6720 1 8588 1519 1 2756 9468 1 5810 9933 1 1625 9138 0 27 932 0 8124 1097 1 8699 8363 1 7202 9489 1 9420 7337 1 9862 6164 1 908 2128 1 9341 2521 1 372 8140 0 9318 7520 1 1760 34 0 785 612 0 1445 1442 0 9694 3280 1 9496 2152 0 9249 1891 0 2162 9828 0 580 2663 1 1269 832...
output:
234
result:
ok single line: '234'
Test #7:
score: -100
Wrong Answer
time: 81ms
memory: 13812kb
input:
100000 7487 4751 1 7487 4751 1 7487 4752 1 7487 4752 1 7487 4752 1 7487 4752 1 7487 4753 1 7487 4753 1 7487 4753 1 7487 4754 1 7487 4754 1 7487 4754 1 7487 4755 1 7487 4755 1 7487 4755 1 7487 4756 1 7488 4756 1 7488 4756 1 7488 4757 1 7488 4757 1 7488 4757 1 7488 4757 1 7488 4758 1 7488 4758 1 7488 ...
output:
68839
result:
wrong answer 1st lines differ - expected: '68833', found: '68839'