QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#806027#9870. ItemswyhaoWA 242ms34524kbC++142.4kb2024-12-08 20:41:452024-12-08 20:42:00

Judging History

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

  • [2024-12-08 20:42:00]
  • 评测
  • 测评结果:WA
  • 用时:242ms
  • 内存:34524kb
  • [2024-12-08 20:41:45]
  • 提交

answer

#include<bits/stdc++.h>
using namespace std;
const int N=500005,P=998244353;
typedef long long ll;
namespace NTT{
    int r[N],n,inv;
    int pw(int x,int k){
        int ans=1;
        for(;k;k>>=1,x=1ll*x*x%P){
            if(k&1) ans=1ll*ans*x%P;
        }
        return ans;
    }
    void init(int lim){
        n=lim;inv=pw(n,P-2);
        for(int i=0;i<lim;i++){
            r[i]=(r[i>>1]>>1)|((i&1)?(n>>1):0);
        }
    }
    void fft(int *f,int v){
        for(int i=0;i<n;i++){
            if(i<r[i]){
                int z=f[i];f[i]=f[r[i]];f[r[i]]=z;
            }
        }
        for(int lim=2;lim<=n;lim<<=1){
            int k=lim>>1,dw=pw(3,(P-1)/lim),w;
            if(v==-1) dw=pw(dw,P-2);
            for(int i=0;i<n;i+=lim){
                w=1;
                for(int j=i;j<i+k;j++){
                    int x=f[j],y=1ll*f[j+k]*w%P;
                    f[j]=(x+y)%P;
                    f[j+k]=(x-y+P)%P;
                    w=1ll*w*dw%P;
                }
            }
        }
        if(v==-1){
            for(int i=0;i<n;i++){
                f[i]=1ll*f[i]*inv%P;
            }
        }
    }
    int X[N],Y[N];
    void tms(int *f,int *g,int *h,int lim){
        // init(lim*2);
        for(int i=0;i<n;i++) X[i]=f[i],Y[i]=g[i];
        fft(X,1);fft(Y,1);
        for(int i=0;i<n;i++) X[i]=1ll*X[i]*Y[i]%P;
        fft(X,-1);
        for(int i=0;i<n;i++) h[i]=X[i];
    }
}
int n;
ll m;
int w[N];
int f[20][N>>1];
int h[N],g[N];
void solve(){
    scanf("%d%lld",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",&w[i]),w[i]-=m/n;
    m%=n;
    int lim =1;
    while(lim<=(2*n+1)) lim<<=1;
    for(int i=0;i<=18;i++){
        for(int j=0;j<lim;j++){
            f[i][j]=0;
        }
    }
    memset(h,0,sizeof h);
    memset(g,0,sizeof g);
    NTT::init(lim);
    for(int i=1;i<=n;i++) f[0][w[i]+n] = 1;
    for(int i=1;i<=18;i++){
        NTT::tms(f[i-1],f[i-1],h,lim);
        for(int j=0;j<2*n+1;j++) f[i][j] = h[j+n]?1:0;
    }
    g[n] = 1;
    for(int i=0;i<=18;i++){
        if((n>>i)&1){
            NTT::tms(f[i],g,h,lim);
            for(int j=0;j<2*n+1;j++) g[j] = h[j+n]?1:0;
        }
        // for(int j=0;j<lim;j++) printf("%d ",f[i][j]);puts("");
    }
    if(g[n+m]) puts("Yes");
    else puts("No");
}
int main(){
    // freopen("ex.in","r",stdin);
    int _;
    scanf("%d",&_);
    while(_--){
        solve();
    }
    return 0;
}

詳細信息

Test #1:

score: 100
Accepted
time: 0ms
memory: 33544kb

input:

4
5 25
0 0 0 0 5
5 11
4 4 4 5 5
5 0
1 2 3 4 5
5 25
0 1 2 3 4

output:

Yes
No
No
No

result:

ok 4 token(s): yes count is 1, no count is 3

Test #2:

score: -100
Wrong Answer
time: 242ms
memory: 34524kb

input:

1314
1 0
0
1 0
1
1 1
0
1 1
1
2 0
0 0
2 0
0 1
2 0
0 2
2 0
1 0
2 0
1 1
2 0
1 2
2 0
2 0
2 0
2 1
2 0
2 2
2 1
0 0
2 1
0 1
2 1
0 2
2 1
1 0
2 1
1 1
2 1
1 2
2 1
2 0
2 1
2 1
2 1
2 2
2 2
0 0
2 2
0 1
2 2
0 2
2 2
1 0
2 2
1 1
2 2
1 2
2 2
2 0
2 2
2 1
2 2
2 2
2 3
0 0
2 3
0 1
2 3
0 2
2 3
1 0
2 3
1 1
2 3
1 2
2 3
2 0...

output:

Yes
No
No
Yes
Yes
Yes
Yes
Yes
No
No
Yes
No
No
No
Yes
No
Yes
No
No
No
No
No
No
Yes
Yes
Yes
Yes
Yes
Yes
Yes
No
No
No
No
No
No
Yes
No
Yes
No
No
No
Yes
No
No
Yes
Yes
Yes
Yes
Yes
Yes
Yes
Yes
Yes
Yes
Yes
Yes
Yes
Yes
Yes
Yes
Yes
Yes
Yes
Yes
Yes
Yes
Yes
Yes
Yes
No
No
No
Yes
No
No
No
Yes
No
No
No
Yes
Yes
Yes...

result:

wrong answer expected YES, found NO [179th token]