抽球游戏(fwt)

地址：https://nanti.jisuanke.com/t/26017

分析：

现在是给定p，求是否存在这样的数列c，我们可以让p进行fwt变换，然后把点值都三次方根，然后再把得到的点值ufwt成系数

这题主要是判断无解的情况：

1、开三次方根后不是整数

2、最后得到的系数中有负数或者和不为给定的n

3、最后ufwt的过程中出现了非整数

#include<bits/stdc++.h>
using namespace std;
const int maxn=64;
int a[maxn+5];
int n;
void fwt(int *a,int n)
{
    for(int d=1;d<n;d<<=1)
        for(int m=d<<1,i=0;i<n;i+=m)
            for(int j=0;j<d;j++)
            {
                int x=a[i+j],y=a[i+j+d];
                a[i+j]=x+y,a[i+j+d]=x-y;
            }
}
bool ufwt(int *a,int n)
{
    for(int d=1;d<n;d<<=1)
        for(int m=d<<1,i=0;i<n;i+=m)
            for(int j=0;j<d;j++)
            {
                int x=a[i+j],y=a[i+j+d];
               // printf("%d %d
",i+j,i+j+d);
                if((x+y)%2!=0) return 0;
                if((x-y)%2!=0) return 0;
                a[i+j]=(x+y)/2,a[i+j+d]=(x-y)/2;
            }
    return 1;
}
bool solve()
{
    fwt(a,64);
    //for(int i=0;i<64;++i) printf("%d ",a[i]);printf("
");
    for(int i=0;i<64;++i)
    {

        int type=a[i]<0;
        int num=abs(a[i]);

        //printf("%.9f
",pow(-1,1.0/3));
        a[i]=round(pow(num,1.0/3));

        //printf("%d %d %d
",i,a[i],num);
        if(a[i]*a[i]*a[i]!=num) return 0;
        if(type) a[i]=-a[i];
    }
    //printf("ok
");
    if(!ufwt(a,64)) return 0;
    for(int i=0;i<64;++i)
    {
        if(a[i]<0) return 0;
        n-=a[i];
    }
    if(n!=0) return 0;
    for(int i=0;i<64;++i)
        for(int j=1;j<=a[i];++j) printf(" %d",i);
    printf("
");
    return 1;
}
int main()
{
    int T;
    scanf("%d",&T);
    for(int cas=1;cas<=T;++cas)
    {
        printf("Case #%d:",cas);
        scanf("%d",&n);
        for(int i=0;i<64;++i) scanf("%d",&a[i]);
        if(!solve())printf(" -1
");
    }
    return 0;
}