QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#67073 | #5169. 夹娃娃 | hjl666 | 0 | 252ms | 22764kb | C++14 | 4.9kb | 2022-12-09 22:05:56 | 2022-12-09 22:05:59 |
Judging History
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 pw[M][N],con[N][M],fac[M+5],inv[M+5];
int Qpow(int b,int p){
int 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++){
int res=0;
for(int k=0;k<=j;k++)(res+=f[x][k]*f[y][j-k]%mod)%=mod;
f[S][j]=res;
}
// for(int j=0;j<=10;j++)cout<<f[S][j]<<' ';cout<<"\n";
}
for(int S=0;S<=tot;S++){
int g[M+5];
for(int j=0;j<M;j++){
int res=0;
for(int k=0;k<=m;k++)(res+=f[S][k]*pw[j][k]%mod)%=mod;
g[j]=res;
}
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]=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]=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();int p[M+5],q[M+5];
p[1]=1;
for(int i=1;i<M;i++){
for(int j=min(i,m);j;j--)(p[j+1]+=p[j])%=mod,p[j]=p[j]*(mod-i)%mod;
}
// for(int i=0;i<=10;i++)cout<<p[i]<<' ';cout<<"\n";
for(int i=0;i<M;i++){
int down=((M-1-i)&1?-1:1)*inv[i]*inv[M-1-i]%mod;
if(i==0)memcpy(q,p+1,sizeof p);
else {
memcpy(q,p,sizeof p);
int iv=-inv[i]*fac[i-1]%mod;
for(int j=0;j<=m;j++)q[j+1]=(q[j+1]-((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 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]){
int res=0,m=Q[t][0];int sz[B];
for(int i=0;i<B;i++)sz[i]=(Q[t][1]>>shop[i].l)&((1<<shop[i].n)-1);
// cout<<sz[0]<<' '<<sz[1]<<' '<<sz[2]<<"\n";
for(int i=0;i<M;i++){
res+=con[m][i]*shop[0].g[sz[0]][i]%mod*shop[1].g[sz[1]][i]%mod*shop[2].g[sz[2]][i]%mod;
// 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";
res%=mod;
}
Q[t][1]=(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][1]<<"\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;
}
詳細信息
Subtask #1:
score: 0
Wrong Answer
Test #1:
score: 0
Wrong Answer
time: 252ms
memory: 22764kb
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:
26288071 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:
wrong answer 1st lines differ - expected: '38813347', found: '26288071'
Subtask #2:
score: 0
Skipped
Dependency #1:
0%
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:
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%