QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#67076#5169. 夹娃娃hjl6660 1821ms19264kbC++145.0kb2022-12-09 22:51:112022-12-09 22:51:12

Judging History

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

  • [2023-08-10 23:21:45]
  • System Update: QOJ starts to keep a history of the judgings of all the submissions.
  • [2022-12-09 22:51:12]
  • 评测
  • 测评结果:0
  • 用时:1821ms
  • 内存:19264kb
  • [2022-12-09 22:51:11]
  • 提交

answer

#include<bits/stdc++.h>
#define db double
//#define int ll
#define ll long long
#define ull unsigned long long
#define pb emplace_back
#define MP make_pair
#define pii pair<int, int>
#define fi first
#define se second
#define ls k<<1
#define rs k<<1|1
#define CLK (double)clock()/(double)CLOCKS_PER_SEC
using namespace std;
mt19937 rnd(time(0));
inline int read(){
	register int x=0,f=1;
	register char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}
inline void write(register int x){
	if(x<0){putchar('-');x=-x;}
	if(x>9)write(x/10);
	putchar(x%10+'0');
}
const int m=520,N=525,M=2080,B=3;
int n,q,mod,a[15][N],Q[52100][2];
vector<int>V[N];
char s[20];
int con[N][M];
ll pw[M][N],fac[M+5],inv[M+5];
ll Qpow(ll b,int p){
	ll res=1ll;
	while(p){if(p&1)(res*=b)%=mod;(b*=b)%=mod;p>>=1;}
	return res;
}
struct node{
	int n,l,r;
	int f[32][M+5],g[32][M+5],h[32][M+5];
	void init(int L,int R){
		l=L,r=R,n=R-L;
		f[0][0]=1;//cout<<L<<' '<<R<<"\n";
		for(int i=0;i<n;i++)memcpy(f[1<<i],a[L+i],sizeof a[L+i]);
		int tot=(1<<n)-1;
		for(int S=1;S<=tot;S++){
			int x=S&(-S),y=S^x;if(!y)continue;
			for(int j=0;j<=m;j++){
				__int128 res=0;
				for(int k=0;k<=j;k++)res+=f[x][k]*f[y][j-k];
				f[S][j]=res%mod;
			} 
//			for(int j=0;j<=10;j++)cout<<f[S][j]<<' ';cout<<"\n";
		}
		for(int S=0;S<=tot;S++){
			static int g[M+5];
			for(int j=0;j<M;j++){
				__int128 res=0;
				for(int k=0;k<=m;k++)res+=f[S][k]*pw[j][k];
				g[j]=res%mod;
			}
			memcpy(f[S],g,sizeof g);
//			for(int j=0;j<=10;j++)cout<<f[S][j]<<' ';cout<<"\n";
		}
	}
	void work(int k){
		int tot=(1<<n)-1,ok=0;
		for(int S=0;S<=tot;S++)if(__builtin_popcount(S)*k<=m)ok|=1<<S;
//		cout<<k<<' '<<ok<<"\n";
		for(int i=0;i<M;i++)g[0][i]=1;
		for(int i=0;i<n;i++){
			for(int j=0;j<M;j++){
				(h[i][j]+=a[l+i][k-1]*pw[j][k-1]%mod)%=mod;
				g[1<<i][j]=-h[i][j];
			}
		}
		for(int S=0;S<=tot;S++){
			if(ok>>S&1^1)continue;
			int x=S&(-S),y=S^x;if(!y)continue;
			for(int j=0;j<M;j++)g[S][j]=1ll*g[x][j]*g[y][j]%mod;
//			for(int j=0;j<=10;j++)cout<<g[S][j]<<' ';cout<<"\n";
		}
		for(int S=0;S<=tot;S++){
			for(int j=0;j<M;j++)g[S][j]=1ll*g[S][j]*f[tot^S][j]%mod;
//			if(S==3){for(int j=0;j<=10;j++)cout<<g[S][j]<<' ';cout<<"\n";}
		}
		for(int h=2;h<=(1<<n);h<<=1){
			for(int i=0;i<(1<<n);i+=h){
				for(int j=i;j<i+h/2;j++){
					if(ok>>(j+h/2)&1^1)continue;
					for(int l=0;l<M;l++)g[j+h/2][l]=(g[j+h/2][l]+g[j][l])%mod;
				}
			}
		}
//		for(int i=2;i<=2;i++){
//			for(int j=0;j<=10;j++)cout<<g[i][j]<<' ';cout<<"\n";
//		}
	}
}shop[B];
void init(){
	fac[0]=1;
	for(int i=1;i<=M;i++)fac[i]=fac[i-1]*i%mod;
	inv[M]=Qpow(fac[M],mod-2);
	for(int i=M-1;i>=0;i--)inv[i]=inv[i+1]*(i+1)%mod;
	for(int i=0;i<M;i++){
		pw[i][0]=1;
		for(int j=1;j<=m;j++)pw[i][j]=pw[i][j-1]*i%mod;
	}
}
void Solve(){
	init();static int p[m+5],q[m+5];
	p[1]=1;
	for(int i=1;i<M;i++){
		for(int j=min(i,m+1);j;j--)(p[j+1]+=p[j])%=mod,p[j]=1ll*p[j]*(mod-i)%mod;
	}
//	for(int i=m-10;i<=m;i++)cout<<i<<' '<<p[i]<<"\n";
	for(int i=0;i<M;i++){
		ll down=((M-1-i)&1?mod-1:1)*inv[i]%mod*inv[M-1-i]%mod;//cout<<down<<"\n";
		if(i==0)memcpy(q,p+1,sizeof p);
		else {
			memcpy(q,p,sizeof p);
			ll iv=-inv[i]*fac[i-1]%mod;
			for(int j=0;j<=m;j++)q[j+1]=(q[j+1]-(q[j]=q[j]*iv%mod))%mod;
		}
		for(int j=0;j<=m;j++)con[j][i]=((j>0?con[j-1][i]:0)+q[j]*down)%mod;
		for(int j=0;j<=m;j++)con[j][i]=(con[j][i]+mod)%mod;
//		for(int j=0;j<=m;j++)con[j+1][i]=con[j][i]=(con[j][i]+q[j]*down%mod)%mod;
	}
//	for(int i=0;i<=10;i++){
//		for(int j=0;j<=10;j++)cout<<con[i][j]<<' ';cout<<"\n";
//	}
	for(int i=0;i<B;i++)shop[i].init(i*n/B,(i+1)*n/B);
	for(int k=1;k<=m;k++){
		for(int i=0;i<B;i++)shop[i].work(k);
		for(auto t:V[k]){
			__int128 res=0;;
			int m=Q[t][0],sz[B];
			for(int i=0;i<B;i++)sz[i]=(Q[t][1]>>shop[i].l)&((1<<shop[i].n)-1);
//			cout<<m<<' '<<sz[0]<<' '<<sz[1]<<' '<<sz[2]<<"\n";
			for(int i=0;i<M;i++){
				res+=(__int128)con[m][i]*shop[0].g[sz[0]][i]*shop[1].g[sz[1]][i]*shop[2].g[sz[2]][i];
//				if(i<=3)cout<<con[m][i]<<' '<<shop[0].g[sz[0]][i]%mod*shop[1].g[sz[1]][i]%mod*shop[2].g[sz[2]][i]%mod<<"\n";
				if(!(i&0xff))res%=mod;
			}
			Q[t][0]=(res+mod)%mod;
		}
	}
}
void MAIN(){
	n=read();q=read();mod=read();
	for(int id=0;id<n;id++){
		int t=read();a[id][0]=1;
		for(int _=1;_<=t;_++){
			int b=read(),c=read();
			for(int i=b;i<=m;i++)(a[id][i]+=a[id][i-b])%=mod;
			b*=(c+1);
			for(int i=m;i>=b;i--)(a[id][i]-=a[id][i-b])%=mod; 
		}
//		for(int i=0;i<=10;i++)cout<<a[id][i]<<' ';cout<<"\n";
	}
	for(int i=1;i<=q;i++){
		scanf("%s",s);Q[i][0]=read();int k=read();
		for(int j=0;j<n;j++)if(s[j]=='1')Q[i][1]|=1<<j;
		if(__builtin_popcount(Q[i][1])*k>Q[i][0])Q[i][0]=0;
		else V[k].pb(i);
	}
	Solve();
	for(int i=1;i<=q;i++)cout<<Q[i][0]<<"\n";
}
signed main(){
//	freopen("a.in","r",stdin);
//	freopen("a.out","w",stdout);
    int T=1;while(T--)MAIN();
//  printf("\nTIME:%lf\n",(double)clock()/CLOCKS_PER_SEC);
	return 0;
}




Details

Tip: Click on the bar to expand more detailed information

Subtask #1:

score: 0
Wrong Answer

Test #1:

score: 0
Wrong Answer
time: 122ms
memory: 18408kb

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
-959431015
98254946
922143620
98254933
922143604
-959431051
-959431065
-76100791
38813247
-959431133
-959431165
-959431203
98254715
38813047
922143273
98254516
38812814
922142999
-899990162
-76101630
38812257
38812058
-899990917
922141847
38811240
92214...

result:

wrong answer 6th lines differ - expected: '38813338', found: '-959431015'

Subtask #2:

score: 0
Skipped

Dependency #1:

0%

Subtask #3:

score: 0
Skipped

Dependency #2:

0%

Subtask #4:

score: 0
Wrong Answer

Test #9:

score: 0
Wrong Answer
time: 1821ms
memory: 19264kb

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
-131875766
797687294
-516999177
-516999177
-959221765
269889529
-445466118
-228421765
331666941
99789529
903956470
112750588
756797435
519045176
-127332353
-636661765
-403963906
494747647
597778941
178845176
-562788353
493445999
-536510471
308912117
-727057412
-501554824
-7873305...

result:

wrong answer 3rd lines differ - expected: '866368587', found: '-131875766'

Subtask #5:

score: 0
Skipped

Dependency #3:

0%

Subtask #6:

score: 0
Skipped

Dependency #4:

0%

Subtask #7:

score: 0
Skipped

Dependency #6:

0%