QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#411928 | #8669. 正方形计数 | grass8cow | 100 ✓ | 625ms | 66660kb | C++17 | 4.1kb | 2024-05-15 21:41:10 | 2024-05-15 21:41:11 |
Judging History
answer
#include<bits/stdc++.h>
using namespace std;
const int V=2e3;
#define pi pair<int,int>
#define mp make_pair
#define fi first
#define se second
#define ll long long
#define db long double
const int I=1e9;
int p[4],q[4];
const db eps=1e-13;
ll ans;
ll xq(ll x,ll y){
if(x>=0)return x/y;
return (x+1)/y-1;
}
struct qb{
ll x,y;
inline ll fk(){return xq(x,y);}
};
struct TB{
bool os;
pi a[20];int m,ax[20],ay[20];
ll as[2010][2010];
int cp[20];
inline void build(bool F){
os=F;sort(a+1,a+m+1);m=unique(a+1,a+m+1)-a-1;
m--;
for(int i=1;i<=m;i++){
ax[i]=a[i+1].fi-a[i].fi,ay[i]=a[i+1].se-a[i].se;
cp[i]=cp[i-1]+ax[i];
for(int j=0;j<ax[i];j++){
for(int k=0;k<=V;k++){
as[j+cp[i-1]][k]=xq(j+ay[i]*k,ax[i]);
if(k)as[j+cp[i-1]][k]+=as[j+cp[i-1]][k-1];
}
}
}
}
inline void ig(int x,int y){a[++m]=mp(x,y);}
int sta[20],tp,b[20];
inline qb ik(int i,int j){
qb e=(qb){1ll*b[i]*ax[j]-1ll*b[j]*ax[i],1ll*ay[j]*ax[i]-ay[i]*ax[j]};
if(e.y<0)e.x=-e.x,e.y=-e.y;
return e;
}
inline bool chk(qb a,qb b){
return a.x*b.y>=b.x*a.y;
}
int iq[20];
db B[20],K[20];db ix[20];
inline void gt(){
tp=0;
for(int i=1;i<=m;i++){
b[i]=(os?I:-I);
for(int j=0;j<4;j++){
int x=a[i].fi+p[j],y=a[i].se+q[j],vb=y*ax[i]-x*ay[i];
b[i]=(os?min(b[i],vb):max(b[i],vb));
//(x,y)
}
if(!os)b[i]--;
B[i]=(db)b[i]/ax[i],K[i]=(db)ay[i]/ax[i];
//b[i]*ax[j]+ay[i]*ax[j]*x=b[j]*ax[i]+ay[j]*ax[i]*x
//该直线:y=(b[i]+ay[i]*x)/ax[i]
while(tp>1&&chk(ik(sta[tp-1],sta[tp]),ik(sta[tp-1],i)))tp--;
sta[++tp]=i;
}
for(int i=1;i<tp;i++){
qb A=ik(sta[i],sta[i+1]);
iq[i]=A.fk(),ix[i]=1.0*A.x/A.y;
}
}
inline void sol(int l,int r){
for(int i=1;i<=tp;i++){
int u=sta[i],xl=l,xr=r;
if(i>1)xl=max(xl,iq[i-1]+1);
if(i<tp)xr=min(xr,iq[i]);
if(xl>xr)continue;
ll ct=xq(b[u],ax[u]),z=b[u]-ax[u]*ct+cp[u-1];
if(os)ans+=(xr-xl+1)*ct,ans+=as[z][xr]-(xl?as[z][xl-1]:0);
else ans-=(xr-xl+1)*ct,ans-=as[z][xr]-(xl?as[z][xl-1]:0);
/*for(int j=xl;j<=xr;j++){
ll cj=xq(b[u]+ay[u]*j,ax[u]);
if(os)ans+=cj;else ans-=cj;
}*/
}
}
}A,B;
void sol(int a,int b,int l,int r){
p[0]=q[0]=0,p[1]=-a,q[1]=-b,p[2]=b-a,q[2]=-a-b,p[3]=b,q[3]=-a;
A.gt(),B.gt();
db tl=V+1,tr=-1,lst=-1e9;
int x=1,y=1;
while(1){
db oe=1e9;
if(x<A.tp)oe=min(oe,A.ix[x]);if(y<B.tp)oe=min(oe,B.ix[y]);
int u=A.sta[x],v=B.sta[y];
db al=A.K[u]*lst+A.B[u],ar=A.K[u]*oe+A.B[u];
db bl=B.K[v]*lst+B.B[v],br=B.K[v]*oe+B.B[v];
if(al+eps>=bl&&ar+eps>=br)tl=min(tl,lst),tr=max(tr,oe);
else if(al+eps<bl&&ar+eps<br);
else{
db sd=(B.B[v]-A.B[u])/(A.K[u]-B.K[v]);
if(A.K[u]>B.K[v])tr=max(tr,oe),tl=min(tl,sd);
else tl=min(tl,lst),tr=max(tr,sd);
}
lst=oe;
if(x==A.tp&&y==B.tp)break;
if(x<A.tp&&(y==B.tp||A.ix[x]<B.ix[y]))x++;else y++;
}
tl=max(tl,(db)l),tr=min(tr,(db)r);
int kl=(int)ceil(tl),kr=(int)floor(tr);
l=kl,r=kr;if(l>r)return;
A.sol(l,r),B.sol(l,r);
}
int x_[20],y_[20],n;
int main(){
scanf("%d",&n);
int l_=V+1,r_=-1;
for(int i=1;i<=n;i++)scanf("%d%d",&x_[i],&y_[i]),
l_=min(l_,x_[i]),r_=max(r_,x_[i]);
for(int i=1;i<=n;i++)if(x_[i]<x_[i%n+1])A.ig(x_[i],y_[i]),A.ig(x_[i%n+1],y_[i%n+1]);
for(int i=1;i<=n;i++)if(x_[i]>x_[i%n+1])B.ig(x_[i],y_[i]),B.ig(x_[i%n+1],y_[i%n+1]);
A.build(1),B.build(0);
for(int a=1;a<=V;a++)for(int b=0;a+b<=V;b++){
int L=l_+b,R=r_-a;if(L>R)continue;
sol(a,b,L,R);
}
printf("%lld",ans);
return 0;
}
Details
Tip: Click on the bar to expand more detailed information
Subtask #1:
score: 10
Accepted
Test #1:
score: 10
Accepted
time: 123ms
memory: 60044kb
input:
4 131 603 131 1828 1919 1828 1919 603
output:
361182910200
result:
ok 1 number(s): "361182910200"
Test #2:
score: 10
Accepted
time: 2ms
memory: 4456kb
input:
4 239 211 239 962 261 962 261 211
output:
1498772
result:
ok 1 number(s): "1498772"
Test #3:
score: 10
Accepted
time: 190ms
memory: 66632kb
input:
4 0 0 0 2000 2000 2000 2000 0
output:
1336001667000
result:
ok 1 number(s): "1336001667000"
Test #4:
score: 10
Accepted
time: 24ms
memory: 23816kb
input:
4 36 771 36 786 672 786 672 771
output:
427720
result:
ok 1 number(s): "427720"
Test #5:
score: 10
Accepted
time: 9ms
memory: 9912kb
input:
4 0 100 100 200 200 100 100 0
output:
34001650
result:
ok 1 number(s): "34001650"
Subtask #2:
score: 25
Accepted
Test #6:
score: 25
Accepted
time: 135ms
memory: 59948kb
input:
3 131 603 131 1828 1919 603
output:
63739309181
result:
ok 1 number(s): "63739309181"
Test #7:
score: 25
Accepted
time: 2ms
memory: 4588kb
input:
3 239 211 239 962 261 211
output:
353073
result:
ok 1 number(s): "353073"
Test #8:
score: 25
Accepted
time: 149ms
memory: 66660kb
input:
3 0 0 0 2000 2000 0
output:
222889277611
result:
ok 1 number(s): "222889277611"
Test #9:
score: 25
Accepted
time: 15ms
memory: 23796kb
input:
3 36 771 36 786 672 771
output:
98847
result:
ok 1 number(s): "98847"
Test #10:
score: 25
Accepted
time: 4ms
memory: 6936kb
input:
3 0 0 0 100 100 0
output:
1473186
result:
ok 1 number(s): "1473186"
Subtask #3:
score: 15
Accepted
Test #11:
score: 15
Accepted
time: 2ms
memory: 4244kb
input:
8 0 13 4 15 15 15 15 6 13 1 12 0 5 0 0 6
output:
4047
result:
ok 1 number(s): "4047"
Test #12:
score: 15
Accepted
time: 2ms
memory: 4308kb
input:
8 0 4 1 15 2 15 15 14 15 4 14 0 1 0 0 2
output:
4200
result:
ok 1 number(s): "4200"
Test #13:
score: 15
Accepted
time: 0ms
memory: 4312kb
input:
5 7 15 15 13 15 0 3 0 0 15
output:
3635
result:
ok 1 number(s): "3635"
Test #14:
score: 15
Accepted
time: 2ms
memory: 4316kb
input:
8 0 12 2 14 7 15 13 15 15 10 15 1 8 0 0 0
output:
4511
result:
ok 1 number(s): "4511"
Test #15:
score: 15
Accepted
time: 2ms
memory: 4244kb
input:
6 0 11 3 15 7 15 15 12 10 0 0 0
output:
3006
result:
ok 1 number(s): "3006"
Test #16:
score: 15
Accepted
time: 2ms
memory: 3912kb
input:
5 0 0 0 2 1 2 2 1 2 0
output:
4
result:
ok 1 number(s): "4"
Subtask #4:
score: 20
Accepted
Dependency #3:
100%
Accepted
Test #17:
score: 20
Accepted
time: 17ms
memory: 13164kb
input:
8 49 299 144 300 300 260 250 15 115 0 30 0 23 19 0 85
output:
443602646
result:
ok 1 number(s): "443602646"
Test #18:
score: 20
Accepted
time: 17ms
memory: 13316kb
input:
8 0 133 103 300 130 300 257 294 297 227 300 150 277 40 161 4
output:
351466521
result:
ok 1 number(s): "351466521"
Test #19:
score: 20
Accepted
time: 16ms
memory: 13168kb
input:
8 76 286 114 300 300 300 300 205 291 0 47 0 4 57 2 235
output:
605026927
result:
ok 1 number(s): "605026927"
Test #20:
score: 20
Accepted
time: 21ms
memory: 13256kb
input:
8 0 102 40 274 282 300 300 234 267 0 34 0 6 57 0 86
output:
497330741
result:
ok 1 number(s): "497330741"
Test #21:
score: 20
Accepted
time: 18ms
memory: 13260kb
input:
7 0 288 156 300 212 300 265 176 300 86 278 0 0 36
output:
446722651
result:
ok 1 number(s): "446722651"
Subtask #5:
score: 15
Accepted
Dependency #4:
100%
Accepted
Test #22:
score: 15
Accepted
time: 63ms
memory: 27548kb
input:
5 257 800 766 800 800 353 667 0 42 0
output:
18881369614
result:
ok 1 number(s): "18881369614"
Test #23:
score: 15
Accepted
time: 103ms
memory: 28372kb
input:
8 691 800 737 795 800 651 372 98 136 266 118 318 24 629 12 753
output:
8760058886
result:
ok 1 number(s): "8760058886"
Test #24:
score: 15
Accepted
time: 72ms
memory: 27232kb
input:
8 718 800 740 800 800 726 800 670 711 367 595 150 86 0 57 136
output:
3064355626
result:
ok 1 number(s): "3064355626"
Test #25:
score: 15
Accepted
time: 93ms
memory: 28768kb
input:
8 0 347 16 449 364 798 674 800 750 800 797 14 195 0 0 70
output:
23587042437
result:
ok 1 number(s): "23587042437"
Test #26:
score: 15
Accepted
time: 102ms
memory: 28916kb
input:
8 322 800 596 800 686 777 800 280 764 69 396 0 46 179 0 660
output:
23185884331
result:
ok 1 number(s): "23185884331"
Subtask #6:
score: 15
Accepted
Dependency #1:
100%
Accepted
Dependency #2:
100%
Accepted
Dependency #3:
100%
Accepted
Dependency #4:
100%
Accepted
Dependency #5:
100%
Accepted
Test #27:
score: 15
Accepted
time: 416ms
memory: 59680kb
input:
8 0 1150 314 2000 1101 2000 1617 1607 1778 551 738 0 607 10 0 1011
output:
577130875850
result:
ok 1 number(s): "577130875850"
Test #28:
score: 15
Accepted
time: 578ms
memory: 66012kb
input:
8 0 1841 1526 2000 1981 1680 1968 678 1893 26 973 0 616 315 524 434
output:
735496008519
result:
ok 1 number(s): "735496008519"
Test #29:
score: 15
Accepted
time: 406ms
memory: 66612kb
input:
6 0 258 10 2000 1730 2000 2000 1510 1973 0 0 129
output:
1203935109430
result:
ok 1 number(s): "1203935109430"
Test #30:
score: 15
Accepted
time: 560ms
memory: 66656kb
input:
7 200 2000 1686 2000 1951 1878 2000 863 1422 0 21 0 0 1015
output:
1100462975231
result:
ok 1 number(s): "1100462975231"
Test #31:
score: 15
Accepted
time: 625ms
memory: 66644kb
input:
8 701 2000 1449 2000 1847 1928 2000 1496 1987 668 1588 108 263 0 0 1985
output:
997591862206
result:
ok 1 number(s): "997591862206"
Test #32:
score: 15
Accepted
time: 574ms
memory: 65488kb
input:
8 15 2000 1235 2000 1545 1886 1970 1526 1828 427 1238 97 372 0 0 1786
output:
816089046494
result:
ok 1 number(s): "816089046494"
Test #33:
score: 15
Accepted
time: 421ms
memory: 66648kb
input:
7 0 1685 1331 2000 2000 1941 2000 1310 1757 631 21 113 0 575
output:
633230324466
result:
ok 1 number(s): "633230324466"
Test #34:
score: 15
Accepted
time: 441ms
memory: 66648kb
input:
8 0 650 0 1350 650 2000 1350 2000 2000 1350 2000 650 1350 0 650 0
output:
900037062925
result:
ok 1 number(s): "900037062925"