QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#411928#8669. 正方形计数grass8cow100 ✓625ms66660kbC++174.1kb2024-05-15 21:41:102024-05-15 21:41:11

Judging History

你现在查看的是最新测评结果

  • [2024-05-15 21:41:11]
  • 评测
  • 测评结果:100
  • 用时:625ms
  • 内存:66660kb
  • [2024-05-15 21:41:10]
  • 提交

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"