QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#830957 | #9611. 木桶效应 | fzx | 100 ✓ | 861ms | 44808kb | C++20 | 4.8kb | 2024-12-25 07:52:49 | 2024-12-25 07:52:49 |
Judging History
answer
#include <bits/stdc++.h>
#define int long long
using namespace std;
const int INF=1e6+5;
const int Mod=998244353;
const int INFN=55;
const int INFS=1035;
struct P3{
int x,y,z;
}aa[INF];
int n,m,q,b[INF],c[INF],sum3,ans,bb[INF],L;
int fav[INF],ifav[INF];
int ksm(int x,int y) {
int ba=x%Mod,ans=1;
while (y) {
if (y&1) ans=(ans*ba)%Mod;
ba=(ba*ba)%Mod;y>>=1;
}
return ans;
}
void init() {
fav[0]=1;int N=1e6;
for (int i=1;i<=N;i++) fav[i]=fav[i-1]*i%Mod;
ifav[N]=ksm(fav[N],Mod-2);
for (int i=N-1;~i;i--) ifav[i]=ifav[i+1]*(i+1)%Mod;
}
int C(int x,int y) {
if (x<y) return 0;
return fav[x]*ifav[y]%Mod*ifav[x-y]%Mod;
}
int vis6[INFN],L6[INFN],B[INFN][INFN];
int tot3,num[INFN][INFN],s5[INFN],S5[INFS],num2[INFN][INFN],sum4i[INFN][INFN],sum4[INFN][INFN],len[INFN];
int f[INFN][INFN][INFS],vis0[INFN][INFS],sum5[INFN][13][INFN],sum5i[INFN][13][INFN];
int C6[INFN*2][INFN*2];
signed main()
{
ios::sync_with_stdio(false);
cin>>n>>m>>q;init();
C6[0][0]=1;
for (int i=1;i<=100;i++) {
C6[i][0]=C6[i][i]=1;
for (int j=1;j<i;j++)
C6[i][j]=C6[i-1][j-1]+C6[i-1][j],C6[i][j]%=Mod;
}
for (int i=1;i<=q;i++)
cin>>aa[i].x>>aa[i].y>>aa[i].z,vis6[aa[i].y]=1;
int sum6=0;
for (int w=1;w<=n;w++)
if (vis6[w]) sum6+=vis6[w],L6[++L6[0]]=w,L+=vis6[w];
sort(aa+1,aa+1+q,[](P3 xx,P3 yy){return xx.x<yy.x;});
for (int i=1;i<=q;i++) {
int r=i;
while (r+1<=q && aa[r+1].x==aa[r].x) r++;
++tot3;
for (int w=i;w<=r;w++) {
num[tot3][aa[w].y]=aa[w].z,len[tot3]++;
// cerr<<num[tot3][aa[w].y]<<" "<<tot3<<" "<<aa[w].y<<" QAQWQ?\n";
for (int t=1;t<=L6[0];t++)
if (L6[t]==aa[w].y) s5[tot3]|=(1ll<<(t-1));
}
for (int w=i;w<=r;w++)
for (int q=aa[w].z;q;q--)
num2[tot3][q]++;
i=r;
}
sum3=0;
for (int i=1;i<=q;i++) {
if (aa[i].x!=aa[i-1].x) sum3++;
}
for (int i=0;i<=n;i++) {
sum4[i][n+1]=1;
for (int l=n;l;l--)
sum4[i][l]=sum4[i][l+1]*ksm(l-(i+1)+1,m-sum3),sum4[i][l]%=Mod;
for (int l=1;l<=n+1;l++)
sum4i[i][l]=ksm(sum4[i][l],Mod-2);
// for (int l=1;l<=n;l++) {
// for (int r=l;r<=n;r++) {
// if (l==r) sum4[i][l][r]=ksm(l-(i+1)+1,m-sum3);
// else sum4[i][l][r]=sum4[i][l][r-1]*ksm(r-(i+1)+1,m-sum3)%Mod;
// }
// }
}
for (int i=0;i<n;i++) {
S5[i]=(1<<n)-1;
for (int t=0;t<(1<<L);t++) {
for (int u=0;u<L6[0];u++)
if ((t>>u)&1) {
for (int w=1;w<=tot3;w++) {
if (!((s5[w]>>u)&1)) continue;
// num[w][L6[u+1]]
if (num[w][L6[u+1]]<i+1) {
vis0[i][t]=1;
S5[i]&=~(1<<u);
}
}
}
}
}
for (int i=0;i<=n;i++) {
for (int u=1;u<=tot3;u++) {
// for (int l=1;l<=n;l++) {
// for (int r=l;r<=n;r++) {
// if (l==r) sum5[i][u][l][r]=(r+len[u]-(i+1)+1-num2[u][i+1]);
// else sum5[i][u][l][r]=sum5[i][u][l][r-1]*(r+len[u]-(i+1)+1-num2[u][i+1])%Mod;
// }
// }
sum5[i][u][n+1]=(n+1+len[u]-(i+1)+1-num2[u][i+1]);
for (int l=n;l;l--)
sum5[i][u][l]=sum5[i][u][l+1]*(l+len[u]-(i+1)+1-num2[u][i+1])%Mod;
for (int l=1;l<=n+1;l++)
sum5i[i][u][l]=ksm(sum5[i][u][l],Mod-2);
}
}
f[0][0][0]=1;
for (int i=0;i<n;i++) {
for (int j=0;j<=n-L;j++) {
for (int s=0;s<(1<<L);s++) {
if (!f[i][j][s]) continue;
// 钦定填了 k 列,t 集合列
for (int k=0;k<=n-L-j;k++) {
int S=((1<<L)-s-1)&S5[i];
int ba=C6[n-L-j][k];
int P2=__builtin_popcount(s);
ba*=sum4[i][j+P2+1];ba%=Mod;
for (int u=1;u<=tot3;u++) {
int l=j+__builtin_popcount(s&(~s5[u]))+1;
ba*=sum5[i][u][l],ba%=Mod;
}
f[i+1][j+k][s]+=f[i][j][s]*ba;f[i+1][j+k][s]%=Mod;
// for (int t=S;;t=(t-1)&S) {
// int ba=C6[n-L-j][k],P=__builtin_popcount(t);
// int P2=__builtin_popcount(s);
// ba*=sum4[i][j+P2+1]*sum4i[i][j+P2+k+P+1]%Mod,ba%=Mod;
// if (!ba) {if (!t) break;continue;}
// for (int u=1;u<=tot3;u++) {
// int l=j+__builtin_popcount(s&(~s5[u]))+1,r=j+__builtin_popcount((s+t)&(~s5[u]))+k;
// ba*=sum5[i][u][l]*sum5i[i][u][r+1]%Mod,ba%=Mod;
// }
// if (!ba) {if (!t) break;continue;}
// f[i+1][j+k][s+t]+=f[i][j][s]*ba;
// f[i+1][j+k][s+t]%=Mod;
// if (!t) break;
// }
}
}
}
for (int j=0;j<=n-L;j++) {
for (int p=0;p<L;p++) {
if ((S5[i]>>p)&1) {
for (int s=0;s<(1<<L);s++)
if ((s>>p)&1)
f[i+1][j][s]+=f[i+1][j][s^(1<<p)],f[i+1][j][s]%=Mod;
}
}
}
for (int j=0;j<=n-L;j++)
for (int s=0;s<(1<<L);s++) {
int P=__builtin_popcount(s);
f[i+1][j][s]*=sum4i[i][j+P+1];f[i+1][j][s]%=Mod;
for (int u=1;u<=tot3;u++) {
int r=j+__builtin_popcount(s&(~s5[u]));
f[i+1][j][s]*=sum5i[i][u][r+1]%Mod,f[i+1][j][s]%=Mod;
}
}
}
cout<<f[n][n-L][(1<<L)-1]<<"\n";
return 0;
}
Details
Tip: Click on the bar to expand more detailed information
Subtask #1:
score: 4
Accepted
Test #1:
score: 4
Accepted
time: 13ms
memory: 26220kb
input:
5 3 0
output:
21412920
result:
ok single line: '21412920'
Test #2:
score: 4
Accepted
time: 9ms
memory: 26532kb
input:
5 3 2 3 4 4 2 5 3
output:
847674
result:
ok single line: '847674'
Test #3:
score: 4
Accepted
time: 9ms
memory: 26488kb
input:
5 3 3 3 5 5 1 2 5 2 5 3
output:
168780
result:
ok single line: '168780'
Subtask #2:
score: 8
Accepted
Dependency #1:
100%
Accepted
Test #4:
score: 8
Accepted
time: 6ms
memory: 26432kb
input:
7 3 0
output:
160221085
result:
ok single line: '160221085'
Test #5:
score: 8
Accepted
time: 9ms
memory: 28624kb
input:
7 3 2 1 1 5 1 6 2
output:
598007855
result:
ok single line: '598007855'
Test #6:
score: 8
Accepted
time: 10ms
memory: 26444kb
input:
7 3 3 1 1 5 2 6 3 2 1 7
output:
950880222
result:
ok single line: '950880222'
Subtask #3:
score: 8
Accepted
Test #7:
score: 8
Accepted
time: 7ms
memory: 38360kb
input:
50 1 0
output:
263941435
result:
ok single line: '263941435'
Test #8:
score: 8
Accepted
time: 15ms
memory: 35536kb
input:
43 2 0
output:
136378346
result:
ok single line: '136378346'
Test #9:
score: 8
Accepted
time: 4ms
memory: 36592kb
input:
50 2 0
output:
489087596
result:
ok single line: '489087596'
Subtask #4:
score: 12
Accepted
Dependency #3:
100%
Accepted
Test #10:
score: 12
Accepted
time: 15ms
memory: 36544kb
input:
50 292247015 0
output:
226872193
result:
ok single line: '226872193'
Test #11:
score: 12
Accepted
time: 11ms
memory: 36404kb
input:
50 873009728 0
output:
63648791
result:
ok single line: '63648791'
Subtask #5:
score: 16
Accepted
Dependency #2:
100%
Accepted
Dependency #4:
100%
Accepted
Test #12:
score: 16
Accepted
time: 11ms
memory: 27568kb
input:
20 728836785 5 248289571 15 16 439110385 8 15 339467267 12 7 585491339 7 9 605518440 4 1
output:
761275883
result:
ok single line: '761275883'
Test #13:
score: 16
Accepted
time: 15ms
memory: 27688kb
input:
20 288197925 5 31379347 9 1 222153278 1 1 24531748 20 1 106427339 18 1 212338331 17 1
output:
877851586
result:
ok single line: '877851586'
Test #14:
score: 16
Accepted
time: 3ms
memory: 27772kb
input:
20 880439688 5 563540321 17 20 395236367 3 20 300712779 6 20 121406689 18 20 25890496 9 20
output:
186649553
result:
ok single line: '186649553'
Test #15:
score: 16
Accepted
time: 10ms
memory: 29788kb
input:
20 311402369 5 311293636 14 13 306211116 19 19 36858994 5 11 36858994 19 10 306211116 14 18
output:
415461922
result:
ok single line: '415461922'
Test #16:
score: 16
Accepted
time: 11ms
memory: 27656kb
input:
20 98953332 5 1075868 12 5 31161114 8 12 46790018 9 10 39214697 15 7 46790018 8 8
output:
204149614
result:
ok single line: '204149614'
Subtask #6:
score: 12
Accepted
Dependency #5:
100%
Accepted
Test #17:
score: 12
Accepted
time: 33ms
memory: 36940kb
input:
50 915702052 5 541920465 39 16 833447607 49 14 326677362 14 34 412319566 10 36 206682128 46 28
output:
783441394
result:
ok single line: '783441394'
Test #18:
score: 12
Accepted
time: 26ms
memory: 36736kb
input:
50 47879239 5 30666754 20 1 23945845 17 1 27229141 25 1 40551703 9 1 15723198 18 1
output:
546399382
result:
ok single line: '546399382'
Test #19:
score: 12
Accepted
time: 15ms
memory: 36784kb
input:
50 334191703 5 167838076 19 50 95127599 33 50 221030062 4 50 89223523 36 50 40736662 39 50
output:
242104112
result:
ok single line: '242104112'
Test #20:
score: 12
Accepted
time: 3ms
memory: 35912kb
input:
50 33538004 5 1341436 40 26 19132404 4 13 22271562 24 25 19132404 40 18 19132404 24 38
output:
509216038
result:
ok single line: '509216038'
Test #21:
score: 12
Accepted
time: 23ms
memory: 35820kb
input:
50 453197583 5 304895210 17 41 450727109 33 20 129855576 4 24 214262208 12 46 214262208 4 27
output:
929330226
result:
ok single line: '929330226'
Subtask #7:
score: 20
Accepted
Dependency #6:
100%
Accepted
Test #22:
score: 20
Accepted
time: 96ms
memory: 36964kb
input:
50 733592974 7 204969642 20 32 532101473 38 22 428067108 29 19 699479674 23 35 403248947 11 33 523610998 39 50 207894041 37 47
output:
878602112
result:
ok single line: '878602112'
Test #23:
score: 20
Accepted
time: 101ms
memory: 37384kb
input:
50 363989184 7 275675456 32 1 14440487 41 1 84675645 18 1 112753584 7 1 259289962 30 1 244964889 42 1 34667760 27 1
output:
125183718
result:
ok single line: '125183718'
Test #24:
score: 20
Accepted
time: 85ms
memory: 37736kb
input:
50 987882639 7 17913003 22 50 440453368 26 50 849090802 8 50 106136773 13 50 491841612 34 50 131376512 15 50 231961452 39 50
output:
573542486
result:
ok single line: '573542486'
Test #25:
score: 20
Accepted
time: 19ms
memory: 36704kb
input:
50 739368354 7 163377222 20 27 192827369 29 10 688310083 50 13 20398047 27 45 688310083 29 11 192827369 20 5 688310083 27 27
output:
528550901
result:
ok single line: '528550901'
Test #26:
score: 20
Accepted
time: 32ms
memory: 35984kb
input:
50 691852384 7 552436456 27 11 103126750 8 13 122888842 49 34 306609066 7 11 553011271 38 32 553011271 27 47 553011271 49 50
output:
407673475
result:
ok single line: '407673475'
Test #27:
score: 20
Accepted
time: 43ms
memory: 37372kb
input:
50 828246844 7 528561184 35 4 15522388 43 35 679004541 17 20 161133582 6 20 289045078 32 49 705209075 24 38 705209075 17 44
output:
734432988
result:
ok single line: '734432988'
Subtask #8:
score: 12
Accepted
Dependency #7:
100%
Accepted
Test #28:
score: 12
Accepted
time: 400ms
memory: 43824kb
input:
50 818455715 9 460012830 5 33 118588510 47 47 281020207 23 25 175303600 1 18 234157803 48 32 256460906 19 46 287657461 24 13 81772046 14 22 114821805 44 41
output:
863394340
result:
ok single line: '863394340'
Test #29:
score: 12
Accepted
time: 417ms
memory: 44564kb
input:
50 24969517 9 23565097 13 1 11005666 43 1 21727790 47 1 12763778 6 1 6881701 34 1 2985794 23 1 13971040 15 1 13187518 29 1 1995154 12 1
output:
755357189
result:
ok single line: '755357189'
Test #30:
score: 12
Accepted
time: 397ms
memory: 43564kb
input:
50 429772980 9 414497820 4 50 141139445 19 50 407732495 7 50 83785413 1 50 164351908 8 50 213988430 26 50 52240611 50 50 352159955 27 50 279326344 29 50
output:
247413078
result:
ok single line: '247413078'
Test #31:
score: 12
Accepted
time: 49ms
memory: 36364kb
input:
50 919810726 9 403069670 48 3 193931638 39 28 821686776 22 11 432351038 8 49 466123383 14 34 254851369 2 48 193931638 48 38 821686776 8 19 254851369 8 20
output:
481094724
result:
ok single line: '481094724'
Test #32:
score: 12
Accepted
time: 88ms
memory: 37972kb
input:
50 619619808 9 140356219 25 31 289966305 18 15 308904151 23 44 530131662 15 2 174186843 35 5 452736627 42 7 12231110 21 2 308904151 18 33 140356219 18 30
output:
194056785
result:
ok single line: '194056785'
Test #33:
score: 12
Accepted
time: 183ms
memory: 40860kb
input:
50 405638475 9 189173757 25 33 313860932 22 8 196599730 26 6 273814431 3 2 52701532 15 14 335951766 19 3 313077524 29 13 349363426 7 15 313860932 19 19
output:
435209021
result:
ok single line: '435209021'
Subtask #9:
score: 8
Accepted
Dependency #8:
100%
Accepted
Test #34:
score: 8
Accepted
time: 828ms
memory: 42228kb
input:
50 946059042 10 160844820 11 40 246397981 49 28 373026730 1 18 602144414 9 9 284028279 39 44 22155194 16 2 851925494 6 17 862159117 27 32 773018252 10 20 462257810 48 19
output:
680391855
result:
ok single line: '680391855'
Test #35:
score: 8
Accepted
time: 859ms
memory: 44808kb
input:
50 260661857 10 70899239 47 1 248286199 34 1 49459370 23 1 110741847 3 1 13147285 19 1 75609139 10 1 16001621 49 1 258636033 17 1 41404092 36 1 259336018 9 1
output:
621578303
result:
ok single line: '621578303'
Test #36:
score: 8
Accepted
time: 861ms
memory: 44500kb
input:
50 444766827 10 138606634 2 50 356874907 48 50 186818234 44 50 328471606 8 50 57922803 45 50 442175507 50 50 45412871 41 50 111164719 24 50 278665393 39 50 49548102 15 50
output:
476500313
result:
ok single line: '476500313'
Test #37:
score: 8
Accepted
time: 99ms
memory: 37980kb
input:
50 471185620 10 330350191 32 46 19091103 17 9 383913632 38 25 125638511 4 13 227803510 12 49 140314389 27 34 97824926 22 46 19091103 27 16 140314389 38 4 97824926 12 1
output:
746959227
result:
ok single line: '746959227'
Test #38:
score: 8
Accepted
time: 181ms
memory: 40188kb
input:
50 534235345 10 505481965 11 39 325061903 12 35 443938537 40 20 144968487 29 31 292692014 18 39 423914803 14 36 221863124 46 41 413488879 5 39 413488879 14 5 325061903 11 42
output:
83413094
result:
ok single line: '83413094'
Test #39:
score: 8
Accepted
time: 399ms
memory: 42656kb
input:
50 268783055 10 164040236 5 3 77416439 15 8 215664484 32 34 84807773 8 44 242418459 20 43 132116147 31 13 11369559 50 23 264864429 6 42 153908131 3 27 132116147 20 50
output:
858392519
result:
ok single line: '858392519'
Test #40:
score: 8
Accepted
time: 46ms
memory: 37172kb
input:
50 680178330 10 381317639 8 3 97592656 22 6 66456445 4 9 277935856 2 13 629083753 35 33 86588466 7 25 86588466 8 6 66456445 22 11 97592656 8 23 277935856 7 24
output:
274256402
result:
ok single line: '274256402'