QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#656240 | #5169. 夹娃娃 | Jellyprx | 0 | 132ms | 12156kb | C++14 | 3.5kb | 2024-10-19 11:56:15 | 2024-10-19 11:56:15 |
Judging History
answer
#include<iostream>
#include<cstdio>
#include<vector>
#include<cstring>
#define lld long long
using namespace std;
const lld p=998244353;
const int N=2050,len=2048,lim=520,M=15;
lld qpow(lld a,int y){
lld res=1;
for(;y;y>>=1,a=a*a%p) if(y&1) res=res*a%p;
return res;
}
lld G=114514,invG=qpow(G,p-2);
int rev[N];
void ntt(lld *a,int lim,int typ){
for(int i=0;i<lim;i++) if(i>rev[i]) swap(a[i],a[rev[i]]);
for(int mid=1;mid<lim;mid<<=1){
lld wn=qpow(typ==1?G:invG,(p-1)/(mid<<1));
for(int R=mid<<1,j=0;j<lim;j+=R){
lld w=1;
for(int k=0;k<mid;k++,w=w*wn%p){
lld x=a[j+k],y=w*a[mid+j+k]%p;
a[j+k]=(x+y)%p;
a[mid+j+k]=(x-y+p)%p;
}
}
}
if(typ==1) return ;
const lld inv=qpow(lim,p-2);
for(int i=0;i<lim;i++) a[i]=a[i]*inv%p;
}
int n,Q,trash;
struct Poly{
lld a[N];
}v[M+1],v2[M+1],f[2][1<<8],g[2][1<<8];
lld b[N];
struct Qry{
char s[M+1];int m,k;
}q[52100];int ans[52100];
vector<int>qry[lim+1];
int id[1<<10];
int main(){
// freopen("doll.in","r",stdin);
// freopen("doll.out","w",stdout);
scanf("%d%d%d",&n,&Q,&trash);
for(int i=1;i<len;i++) rev[i]=(rev[i>>1]>>1)|((i&1)<<10);
for(int i=1;i<=n;i++){
v[i].a[0]=1,ntt(v[i].a,len,1);
int m;scanf("%d",&m);
for(int k=1;k<=m;k++){
int val,cnt;scanf("%d%d",&val,&cnt);
for(int j=0;j<=cnt&&val*j<=lim;j++) b[val*j]=1;
// printf("![b]:\n");
// for(int j=0;j<=20;j++) printf("%lld ",b[j]);
ntt(b,len,1);
// putchar('\n');
// for(int j=0;j<=20;j++) printf("%lld ",b[j]);putchar('\n');
// for(int j=0;j<=20;j++) printf("%lld ",v[i].a[j]);putchar('\n');
for(int j=0;j<len;j++) v[i].a[j]=v[i].a[j]*b[j]%p,b[j]=0;
}
for(int j=0;j<len;j++) v2[i].a[j]=v[i].a[j];
ntt(v2[i].a,len,-1);
// printf("\n%d:\n",i);
// for(int j=0;j<=20;j++) printf("%lld ",v2[i].a[j]);
// printf("\n___________\n");
}
int L=1,M=1+n>>1,R=n;
for(int i=1;i<=M;i++) id[1<<i-1]=i;
f[0][0].a[0]=1,ntt(f[0][0].a,len,1);
f[1][0].a[0]=1,ntt(f[1][0].a,len,1);
for(int s=1;s<(1<<M);s++){
int lbit=(s&-s),u=id[lbit],ss=s^lbit;
for(int j=0;j<len;j++) f[0][s].a[j]=f[0][ss].a[j]*v[u].a[j]%p;
}
for(int s=1;s<(1<<R-M);s++){
int lbit=(s&-s),u=id[lbit]+M,ss=s^lbit;
for(int j=0;j<len;j++) f[1][s].a[j]=f[1][ss].a[j]*v[u].a[j]%p;
}
for(int i=1;i<=Q;i++){
scanf("%s%d%d",q[i].s+1,&q[i].m,&q[i].k);
qry[q[i].k].push_back(i);
}
for(int i=lim;i>=1;i--){
if(!qry[i].size()) continue;
memset(g[0][0].a,0,sizeof g[0][0].a),g[0][0].a[0]=1,ntt(g[0][0].a,len,1);
memset(g[1][0].a,0,sizeof g[1][0].a),g[1][0].a[0]=1,ntt(g[1][0].a,len,1);
for(int s=1;s<(1<<M);s++){
int lbit=(s&-s),u=id[lbit],ss=s^lbit;
memset(b,0,sizeof b);
for(int j=i;j<=lim;j++) b[j]=v2[u].a[j];
ntt(b,len,1);
for(int j=0;j<len;j++) g[0][s].a[j]=g[0][ss].a[j]*b[j]%p;
}
for(int s=1;s<(1<<R-M);s++){
int lbit=(s&-s),u=id[lbit]+M,ss=s^lbit;
memset(b,0,sizeof b);
for(int j=i;j<=lim;j++) b[j]=v2[u].a[j];
ntt(b,len,1);
for(int j=0;j<len;j++) g[1][s].a[j]=g[1][ss].a[j]*b[j]%p;
}
for(int x:qry[i]){
int f0=0,g0=0,f1=0,g1=0;
for(int j=1;j<=M;j++){
if(q[x].s[j]=='0') f0|=1<<j-1;
else g0|=1<<j-1;
}
for(int j=1;j<=R-M;j++){
if(q[x].s[j+M]=='0') f1|=1<<j-1;
else g1|=1<<j-1;
}
// printf("%d %d %d %d %d\n",x,f0,g0,f1,g1);
for(int j=0;j<len;j++) b[j]=f[0][f0].a[j]*g[0][g0].a[j]%p*f[1][f1].a[j]%p*g[1][g1].a[j]%p;
ntt(b,len,-1);
for(int j=0;j<=q[x].m;j++) ans[x]=(ans[x]+b[j])%p;
}
}
for(int i=1;i<=Q;i++) printf("%lld\n",ans[i]);
return 0;
}
詳細信息
Subtask #1:
score: 0
Wrong Answer
Test #1:
score: 0
Wrong Answer
time: 132ms
memory: 12156kb
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:
198481100 673328453 68513962 663175913 379942377 455477082 708201926 564752346 225354377 205613639 512248544 513860494 865933396 722470231 514429245 946000542 250894447 78747931 574470029 587410592 696953216 273123636 30209783 323254379 946810402 702936784 988657834 364138161 117209834 242068096 668...
result:
wrong answer 1st lines differ - expected: '38813347', found: '198481100'
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%