QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#87610 | #4903. 细菌 | Pure_Furies | 90 | 1431ms | 17276kb | C++14 | 2.4kb | 2023-03-13 21:14:28 | 2023-03-13 21:14:30 |
Judging History
answer
#include<bits/stdc++.h>
using namespace std;
const int mod=998244353,ggg=3,invggg=(mod+1)/ggg;
const int maxn=524288;
int mypow(int x,int y){
int ret=1;
while(y){
if(y&1)ret=1ll*ret*x%mod;
y>>=1;x=1ll*x*x%mod;
}return ret;
}
int revs[maxn];
void initrev(int Len){
for(int i=1;i<Len;++i){
revs[i]=revs[i>>1]>>1;
if(i&1)revs[i]|=(Len>>1);
}
}
void DFT(int *f,int n,int rev){
int g0=rev==1?ggg:invggg;
initrev(n);
for(register int i=0;i<n;++i)
if(i<revs[i])
f[i]^=f[revs[i]]^=f[i]^=f[revs[i]];
for(register int h=2;h<=n;h<<=1){
int gn=mypow(g0,(mod-1)/h);
for(register int j=0;j<n;j+=h){
int gk=1;
for(register int k=j;k<j+(h>>1);++k){
int e=f[k],o=1ll*gk*f[k+(h>>1)]%mod;
f[k]=(e+o)%mod,f[k+(h>>1)]=(e+mod-o)%mod;
gk=1ll*gk*gn%mod;
}
}
}
if(rev==-1){
int invv=mypow(n,mod-2);
for(register int i=0;i<n;++i)
f[i]=1ll*f[i]*invv%mod;
}
}
int fac[maxn],inv[maxn];
int C(int x,int y){
if(x<y||y<0)return 0;
return 1ll*fac[x]*inv[y]%mod*inv[x-y]%mod;
}
struct node{
int p,l,r,w;
void init(){p=0;l=0;r=0;w=1;}
int calc(int qp,int ql,int qr){
while(p<qp)w=(2ll*w+C(p,l-1)-C(p,r)+mod)%mod,p++;
while(l<ql)w=(w-C(p,l++)+mod)%mod;
while(l>ql)w=(w+C(p,--l))%mod;
while(r>qr)w=(w-C(p,r--)+mod)%mod;
while(r<qr)w=(w+C(p,++r))%mod;
return w;
}
}w[120003];
const int B=700;
int dp[2][B+3],d;
void calc(int *f,int n,int p){
if(n<=B){
memset(dp,0,sizeof(dp));
for(int i=0;i<n;i++)
dp[0][i]=1;
f[0]=1;
for(int i=1;i<=d;i++){
for(int j=0;j<n;j++){
dp[i&1][j]=(j?dp[i-1&1][j-1]:0)+(j<n-1?dp[i-1&1][j+1]:0);
if(dp[i&1][j]>=mod)dp[i&1][j]-=mod;
}f[i]=dp[i&1][p-1];
}return;
}else{
int l=(p-d-n)/(n+1),r=(d+p-1)/(n+1);
for(int i=l;i<=r;i++)
w[i-l].init();
for(int i=0;i<=d;i++)
for(int j=l;j<=r;j++)
f[i]=(f[i]+(j%2?mod-1ll:1ll)*w[j-l].calc(i,(i-p+2+j*(n+1))>>1,(i-p+n)+j*(n+1)>>1))%mod;
}
}
int n,m,k,a,b,c;
int f[maxn],g[maxn],h[maxn];
int main(){
cin>>d>>n>>m>>k>>a>>b>>c;
fac[0]=1;
for(int i=1;i<maxn;i++)
fac[i]=1ll*fac[i-1]*i%mod;
for(int i=0;i<maxn;i++)
inv[i]=mypow(fac[i],mod-2);
calc(f,n,a);
calc(g,m,b);
calc(h,k,c);
for(int i=0;i<=d;i++)
f[i]=1ll*f[i]*inv[i]%mod,
g[i]=1ll*g[i]*inv[i]%mod,
h[i]=1ll*h[i]*inv[i]%mod;
DFT(f,maxn,1);
DFT(g,maxn,1);
DFT(h,maxn,1);
for(int i=0;i<maxn;i++)
f[i]=1ll*f[i]*g[i]%mod*h[i]%mod;
DFT(f,maxn,-1);
cout<<1ll*f[d]*fac[d]%mod;
}
Details
Tip: Click on the bar to expand more detailed information
Subtask #1:
score: 5
Accepted
Test #1:
score: 5
Accepted
time: 193ms
memory: 17276kb
input:
50 41 46 42 8 20 21
output:
791406134
result:
ok 1 number(s): "791406134"
Test #2:
score: 0
Accepted
time: 200ms
memory: 15616kb
input:
50 49 44 48 49 15 25
output:
544847893
result:
ok 1 number(s): "544847893"
Subtask #2:
score: 10
Accepted
Dependency #1:
100%
Accepted
Test #3:
score: 10
Accepted
time: 200ms
memory: 16636kb
input:
5000 4994 4990 4990 976 2257 2505
output:
836390717
result:
ok 1 number(s): "836390717"
Test #4:
score: 0
Accepted
time: 200ms
memory: 17016kb
input:
5000 4994 4997 4994 4399 1826 1332
output:
65414465
result:
ok 1 number(s): "65414465"
Subtask #3:
score: 15
Accepted
Test #5:
score: 15
Accepted
time: 393ms
memory: 17272kb
input:
120000 300 1 1 141 1 1
output:
637174
result:
ok 1 number(s): "637174"
Test #6:
score: 0
Accepted
time: 685ms
memory: 15720kb
input:
120000 994 1 1 15 1 1
output:
218803691
result:
ok 1 number(s): "218803691"
Test #7:
score: 0
Accepted
time: 204ms
memory: 17040kb
input:
120000 119999 1 1 19896 1 1
output:
68846585
result:
ok 1 number(s): "68846585"
Subtask #4:
score: 10
Accepted
Test #8:
score: 10
Accepted
time: 210ms
memory: 16496kb
input:
119000 119991 119991 1 23819 52139 1
output:
944500934
result:
ok 1 number(s): "944500934"
Subtask #5:
score: 15
Accepted
Dependency #4:
100%
Accepted
Test #9:
score: 15
Accepted
time: 269ms
memory: 16056kb
input:
120000 13997 13996 1 42 9266 1
output:
775637357
result:
ok 1 number(s): "775637357"
Test #10:
score: 0
Accepted
time: 280ms
memory: 16832kb
input:
120000 13997 13997 1 9768 6131 1
output:
151873213
result:
ok 1 number(s): "151873213"
Subtask #6:
score: 35
Accepted
Dependency #3:
100%
Accepted
Dependency #5:
100%
Accepted
Test #11:
score: 35
Accepted
time: 557ms
memory: 17264kb
input:
120000 294 296 1 142 273 1
output:
889786082
result:
ok 1 number(s): "889786082"
Test #12:
score: 0
Accepted
time: 672ms
memory: 16288kb
input:
120000 395 390 1 370 185 1
output:
884557050
result:
ok 1 number(s): "884557050"
Test #13:
score: 0
Accepted
time: 788ms
memory: 16640kb
input:
120000 491 495 1 430 25 1
output:
272929924
result:
ok 1 number(s): "272929924"
Test #14:
score: 0
Accepted
time: 924ms
memory: 16804kb
input:
120000 590 593 1 418 459 1
output:
446962505
result:
ok 1 number(s): "446962505"
Test #15:
score: 0
Accepted
time: 1039ms
memory: 15804kb
input:
120000 697 690 1 166 450 1
output:
216092107
result:
ok 1 number(s): "216092107"
Test #16:
score: 0
Accepted
time: 1417ms
memory: 16744kb
input:
120000 793 799 1 416 61 1
output:
661260042
result:
ok 1 number(s): "661260042"
Test #17:
score: 0
Accepted
time: 1182ms
memory: 16068kb
input:
120000 1000 1000 1 613 547 1
output:
429669083
result:
ok 1 number(s): "429669083"
Test #18:
score: 0
Accepted
time: 690ms
memory: 16568kb
input:
120000 1993 1995 1 493 565 1
output:
609392900
result:
ok 1 number(s): "609392900"
Test #19:
score: 0
Accepted
time: 408ms
memory: 15616kb
input:
120000 4995 4998 1 3166 3765 1
output:
394497547
result:
ok 1 number(s): "394497547"
Test #20:
score: 0
Accepted
time: 295ms
memory: 16060kb
input:
120000 9994 9991 1 6099 691 1
output:
795651799
result:
ok 1 number(s): "795651799"
Test #21:
score: 0
Accepted
time: 222ms
memory: 16104kb
input:
120000 49990 49990 1 41933 2862 1
output:
360787513
result:
ok 1 number(s): "360787513"
Test #22:
score: 0
Accepted
time: 202ms
memory: 15680kb
input:
120000 119996 119994 1 58014 49559 1
output:
682979057
result:
ok 1 number(s): "682979057"
Subtask #7:
score: 0
Time Limit Exceeded
Dependency #1:
100%
Accepted
Dependency #2:
100%
Accepted
Dependency #3:
100%
Accepted
Dependency #4:
100%
Accepted
Dependency #5:
100%
Accepted
Dependency #6:
100%
Accepted
Test #23:
score: 10
Accepted
time: 735ms
memory: 15648kb
input:
120000 296 300 297 89 130 280
output:
702788425
result:
ok 1 number(s): "702788425"
Test #24:
score: 0
Accepted
time: 894ms
memory: 15836kb
input:
120000 392 397 391 222 280 10
output:
322470184
result:
ok 1 number(s): "322470184"
Test #25:
score: 0
Accepted
time: 1091ms
memory: 16720kb
input:
120000 499 498 500 263 315 367
output:
449848666
result:
ok 1 number(s): "449848666"
Test #26:
score: 0
Accepted
time: 1255ms
memory: 17112kb
input:
120000 598 591 594 497 474 400
output:
35486519
result:
ok 1 number(s): "35486519"
Test #27:
score: 0
Accepted
time: 1431ms
memory: 15968kb
input:
120000 694 692 695 625 173 477
output:
785203749
result:
ok 1 number(s): "785203749"
Test #28:
score: -10
Time Limit Exceeded
input:
120000 794 796 800 395 465 507