QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#656046#5169. 夹娃娃QBF3 1889ms15432kbC++145.6kb2024-10-19 11:02:442024-10-19 11:02:45

Judging History

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

  • [2024-10-19 11:02:45]
  • 评测
  • 测评结果:3
  • 用时:1889ms
  • 内存:15432kb
  • [2024-10-19 11:02:44]
  • 提交

answer

#include<bits/stdc++.h>
#define ci const int
#define ll long long
#define ms(x) memset(x,0,sizeof(x))
using namespace std;
char buf[1<<20],*p1=buf,*p2=buf;
// #define getchar() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<20,stdin),p1==p2)?EOF:*p1++)
int read(){int res(0);char ch(getchar());while(ch<48||ch>57)ch=getchar();while(ch>=48&&ch<=57)res=(res<<1)+(res<<3)+(ch^48),ch=getchar();return res;}
void out(ci x){
	if(x>9)out(x/10);
	putchar(x%10+48);
}
namespace NTT{
	ci N=525,mod=998244353,g=3,ig=(mod+1)/3;
	inline void add(int &x,ci v){
		x+=v,x-=x<mod?0:mod;
	}
	inline void sub(int &x,ci v){
		x-=v,x+=x<0?mod:0;
	}
	ll qk(ll x,int y=mod-2){
		ll ans=1;
		while(y){
			if(y&1)ans=ans*x%mod;
			x=x*x%mod,y>>=1;
		}
		return ans;
	}
	int R[N<<2];
	void NTT(int *a,ci n,ci tp){
		for(int i=0;i<n;++i)R[i]=(R[i>>1]>>1)|(i&1?(n>>1):0);
		for(int i=0;i<n;++i)if(i<R[i])swap(a[i],a[R[i]]);
		for(int len=2;len<=n;len<<=1){
			ci gn=qk(tp==1?g:ig,(mod-1)/len);
			for(int l=0;l<n;l+=len)
				for(int k=l,al,ar,cur=1;k<l+(len>>1);++k)
					al=a[k],ar=a[k+(len>>1)],
					a[k]=(al+(ll)ar*cur)%mod,
					a[k+(len>>1)]=(al-(ll)ar*cur%mod+mod)%mod,
					cur=(ll)cur*gn%mod;
		}
		if(tp==-1){
			ci inv=qk(n);
			for(int i=0;i<n;++i)a[i]=(ll)a[i]*inv%mod;
		}
	}
	int a[N<<2],b[N<<2],c[N];
	void Mul(int *A,int *B,int *C,ci n){
		int lim=1;
		while(lim<=min(n,510)*2)lim<<=1;
		for(int i=0;i<lim;++i)a[i]=b[i]=0;
		for(int i=0;i<=min(n,510);++i)a[i]=A[i],b[i]=B[i];
		NTT(a,lim,1),NTT(b,lim,1);
		for(int i=0;i<lim;++i)a[i]=(ll)a[i]*b[i]%mod;
		NTT(a,lim,-1);
		for(int i=0;i<=n;++i)c[i]=a[i];
		for(int i=511;i<=n;++i)
			for(int j=0;i+j<=n;++j)
				add(c[i+j],(ll)A[i]*B[j]%mod);
		for(int j=511;j<=n;++j)
			for(int i=0;i+j<=n;++i)
				add(c[i+j],(ll)A[i]*B[j]%mod);
		for(int i=0;i<=n;++i)C[i]=c[i];
	}
}
ci N=525,V=520,H=1.5e5;
int n,q,mod;
int pre[N][N],tmp[N],sum[N];
vector<int>vec[N];
inline void add(int &x,ci v){
	x+=v,x-=x<mod?0:mod;
}
struct DAT{
	map<vector<int>,int>mp;
	int n;
	int a[N],tot,dp[H][N];
	int cnt;
	int pre[N][N];
	int z[1<<8][N];
	void check(ci k){
		int sum=0;
		for(int i=1;i<=n;++i)sum+=a[i];
		if(sum>V)return;
		vector<int>vec;
		for(int i=1;i<=n;++i)vec.push_back(a[i]);
		if(mp.find(vec)!=mp.end())return;
		mp[vec]=++tot;
		int pos=0;
		for(int i=1;i<=n;++i)
			if(a[i]==k){
				pos=i;
				break;
			}
		assert(pos!=0);
		vec[pos-1]=k+1;
		if(mp.find(vec)!=mp.end()){
			ci cur=mp[vec];
			for(int i=0;i<=V;++i)add(dp[tot][i],dp[cur][i]);
		}
		vec[pos-1]=0;
		if(mp.find(vec)!=mp.end()){
			ci cur=mp[vec];
			for(int i=k;i<=V;++i)add(dp[tot][i],(ll)dp[cur][i-k]*pre[pos][k]%mod);
		}
		// bool ok=1;for(int i=1;i<=n;++i)ok&=a[i]!=k+1;
		// if(ok){
		// 	mp2[vec]=++cnt;
		// }
	}
	void dfs(ci d,ci k,ci tp){
		if(d==n+1){
			check(k);
			return;
		}
		a[d]=0,dfs(d+1,k,tp);
		if(k!=0){
			if(tp==0){
				a[d]=k+1,dfs(d+1,k,tp);
				a[d]=k,dfs(d+1,k,1);
			}
			else{
				a[d]=k,dfs(d+1,k,1);
			}
		}
		else{
			for(int k=V;k;--k)
				a[d]=k+1,dfs(d+1,k,0),
				a[d]=k,dfs(d+1,k,1);
		}
	}
	void Solve(ci _n){
		n=_n;
		vector<int>vec;
		for(int i=1;i<=n;++i)vec.push_back(0);
		mp[vec]=++tot;
		dp[1][0]=1;
		dfs(1,0,0);
		z[0][0]=1;
		for(int s=1;s<1<<n;++s){
			ci x=__builtin_ctz(s&-s)+1;
			NTT::Mul(z[s^(s&-s)],pre[x],z[s],V);
			// for(int i=0;i<=V;++i)
			// 	for(int j=0;i+j<=V;++j)
			// 		add(z[s][i+j],(ll)z[s^(s&-s)][i]*pre[x][j]%mod);
		}
	}
}L,R;
int f[N],trs[N],g[N];

int main(){
	// freopen("doll.in","r",stdin);
	// freopen("doll.out","w",stdout);
	int st=clock();
	n=read(),q=read(),mod=read();
	for(int t=1;t<=n;++t){
		pre[t][0]=1;
		for(int i=1;i<=V;++i)vec[i].clear();
		for(int k=read();k;--k){
			ci b=read(),c=read();
			vec[b].push_back(c);
		}
		for(int i=1;i<=V;++i){
			for(int j=0;j<=V;++j)tmp[j]=0;
			tmp[0]=1;
			for(int x:vec[i]){
				for(int k=0;k<=V;++k)sum[k]=((k?sum[k-1]:0)+tmp[k])%mod;
				for(int k=0;k<=V;++k)tmp[k]=(sum[k]-((k-x-1>=0)?sum[k-x-1]:0)+mod)%mod;
			}
			for(int j=V;~j;--j)
				for(int k=1;j+k*i<=V;++k)
					add(pre[t][j+k*i],(ll)pre[t][j]*tmp[k]%mod);
		}
	}
	ci nl=n>>1,nr=n-nl;
	for(int i=1;i<=nl;++i)for(int j=0;j<=V;++j)L.pre[i][j]=pre[i][j];
	for(int i=1;i<=nr;++i)for(int j=0;j<=V;++j)R.pre[i][j]=pre[i+nl][j];
	L.Solve(nl),R.Solve(nr);
	while(q--){
		int st=0,st0=0;
		for(int i=1;i<=n;++i){
			char c=getchar();
			while(c!='0'&&c!='1')c=getchar();
			if(c=='1')st|=1<<i-1;
			else st0|=1<<i-1;
		}
		ci m=read(),k=read();
		for(int i=0;i<=m;++i)f[i]=L.z[st0&((1<<nl)-1)][i];
		NTT::Mul(f,R.z[st0>>nl],f,m);
		vector<int>vec;
		for(int i=1;i<=nl;++i)
			if(st>>(i-1)&1)vec.push_back(k);
			else vec.push_back(0);
		if(L.mp.find(vec)==L.mp.end()){
			puts("0");
			continue;
		}
		else{
			ci cur=L.mp[vec];
			NTT::Mul(f,L.dp[cur],f,m);
			// for(int i=0;i<=m;++i)g[i]=0;
			// for(int i=0;i<=m;++i)trs[i]=L.dp[cur][i];
			// for(int i=0;i<=m;++i)
			// 	for(int j=0;i+j<=m;++j)
			// 		add(g[i+j],(ll)f[i]*trs[j]%mod);
			// for(int i=0;i<=m;++i)f[i]=g[i];
		}
		vec.clear();
		for(int i=1;i<=nr;++i)
			if(st>>(i+nl-1)&1)vec.push_back(k);
			else vec.push_back(0);
		if(R.mp.find(vec)==R.mp.end()){
			puts("0");
			continue;
		}
		else{
			ci cur=R.mp[vec];
			NTT::Mul(f,R.dp[cur],f,m);
			// for(int i=0;i<=m;++i)g[i]=0;
			// for(int i=0;i<=m;++i)trs[i]=R.dp[cur][i];
			// for(int i=0;i<=m;++i)
			// 	for(int j=0;i+j<=m;++j)
			// 		add(g[i+j],(ll)f[i]*trs[j]%mod);
			// for(int i=0;i<=m;++i)f[i]=g[i];
		}
		int ans=0;
		for(int i=0;i<=m;++i)add(ans,f[i]);
		printf("%d\n",ans);
	}
	// printf("%d ms.\n",clock()-st);
	return 0;
}

详细

Subtask #1:

score: 3
Accepted

Test #1:

score: 3
Accepted
time: 130ms
memory: 13144kb

input:

1 521 998244353
39 520 520 11 22 414 8 95 18 229 356 26 407 316 10 24 26 19 61 11 130 482 476 420 15 192 193 208 24 19 233 494 217 275 294 26 28 439 20 272 277 28 198 5 335 22 8 28 17 154 78 6 13 175 17 2 5 477 256 200 4 1 36 427 371 439 23 10 65 426 25 24 27 121 29 28 13 12 453
0 520 1
1 519 1
1 51...

output:

38813347
922143638
98254957
38813343
922143633
38813338
98254946
922143620
98254933
922143604
38813302
38813288
922143562
38813247
38813220
38813188
38813150
98254715
38813047
922143273
98254516
38812814
922142999
98254191
922142723
38812257
38812058
98253436
922141847
38811240
922141173
38810463
38...

result:

ok 521 lines

Test #2:

score: 3
Accepted
time: 408ms
memory: 13424kb

input:

2 1561 998244353
151 520 520 511 30 121 396 25 16 113 11 6 175 242 20 8 5 61 13 518 447 404 8 220 177 4 19 18 15 70 233 9 14 26 512 17 9 9 19 30 8 495 20 13 27 277 22 396 14 4 29 345 442 19 25 14 5 16 295 19 65 134 10 10 296 245 6 7 30 253 15 187 26 482 454 28 414 170 404 11 27 27 25 13 509 1 5 291 ...

output:

883965618
144348435
762074635
112296779
385763651
821718611
673974966
879750066
927942969
136450507
436584627
612945970
768262217
613885343
39304132
852224740
215596261
151746110
965953558
969833936
664053020
458247365
881060255
878484499
781573019
616944059
850325449
296113117
674829177
887392623
6...

result:

ok 1561 lines

Subtask #2:

score: 0
Time Limit Exceeded

Dependency #1:

100%
Accepted

Test #3:

score: 13
Accepted
time: 1027ms
memory: 15432kb

input:

3 4160 998244353
444 520 520 26 332 29 183 25 479 175 14 13 16 1 447 2 293 4 20 64 472 491 11 21 259 75 22 390 401 8 508 405 3 137 4 15 154 164 1 484 13 257 14 44 20 7 13 26 15 26 432 14 9 478 24 18 10 22 28 8 21 260 25 431 22 7 6 20 26 8 27 239 19 1 134 2 322 16 225 6 42 517 6 197 407 268 500 433 5...

output:

516056999
990096150
497048298
345860798
899328070
577475723
191997503
533625761
516056999
863614705
652318084
514747110
811600228
92531482
136793394
218097588
352553395
821305819
739754364
569418540
402235631
844207347
78271439
896568337
516056999
243958673
201200148
634787992
552693501
893938722
98...

result:

ok 4160 lines

Test #4:

score: 13
Accepted
time: 1889ms
memory: 14532kb

input:

4 8320 998244353
303 520 520 288 10 15 24 306 456 495 124 20 419 24 473 7 462 365 405 4 30 1 29 15 25 29 324 407 14 30 184 425 451 6 414 7 417 155 12 18 20 2 475 78 174 467 23 300 26 13 15 345 319 10 27 497 25 21 51 24 485 359 268 87 20 509 13 18 261 13 6 20 237 305 26 245 330 514 29 21 197 25 345 1...

output:

857239630
694514392
340827658
834331936
573150389
560202020
302111919
422193966
147386541
201821565
447255018
322990367
192787601
197802108
461775999
315804262
316164169
338416167
240429979
359914423
321666890
541700460
506123940
701447430
823947537
621301718
62107305
163486246
380210777
211911024
9...

result:

ok 8320 lines

Test #5:

score: 0
Time Limit Exceeded

input:

4 52099 998244353
103 520 520 12 485 23 1 337 514 374 486 210 29 29 1 19 299 3 11 11 22 282 14 12 9 341 286 18 501 3 3 29 364 264 477 22 6 434 14 11 117 22 8 30 268 22 28 10 12 311 58 14 15 234 177 17 238 71 64 14 1 396 23 492 4 1 13 6 8 197 4 7 27 11 370 19 242 12 13 20 185 432 399 24 32 2 516 36 4...

output:

51769626
700322830
226311543
862334239
622358370
748398901
344771107
59670026
558254404
668258250
91212841
493756321
360353830
36696310
321158867
563614481
998042994
565120563
709404804
783802088
511531306
396636746
513730575
451308648
594545675
544172685
900482622
791631384
368742309
537404993
6927...

result:


Subtask #3:

score: 0
Skipped

Dependency #2:

0%

Subtask #4:

score: 0
Time Limit Exceeded

Test #9:

score: 0
Time Limit Exceeded

input:

15 52099 998244353
1 9 3
1 9 4
1 9 2
1 8 10
1 4 4
1 3 1
1 2 5
1 4 9
1 1 4
1 9 4
1 7 6
1 1 6
1 2 5
1 5 2
1 3 5
101000000001010 516 1
010001001010101 520 2
000000101000001 519 2
101011111100011 518 1
010110001000111 520 2
000110111100111 516 1
000100101001011 519 3
000111001010011 518 1
00001110010111...

output:

993379058
496689529
866368587
797687294
481245176
481245176
39022588
269889529
552778235
769822588
331666941
99789529
903956470
112750588
756797435
519045176
870912000
361582588
594280447
494747647
597778941
178845176
435456000
493445999
461733882
308912117
271186941
496689529
919511294
85533882
894...

result:


Subtask #5:

score: 0
Skipped

Dependency #3:

0%

Subtask #6:

score: 0
Skipped

Dependency #4:

0%

Subtask #7:

score: 0
Skipped

Dependency #6:

0%