Spoj 2829 TLE 高维前缀和

title

LUOGU SP2829

给定 \(n\) 和 \(n\) 个数 \(c[i]\) ，要求构造出一串序 列 \(a[i]\) ，要求：

1. \(a[i] \in [0,2^m)\) ；

2. \(a[i] \%c[i]!=0\) ；

3. \(a[i]\&a[i+1]=0\) ；

求方案数（模数 \(10^9\) ）。

analysis

对 \(3\) 的处理：用 \(f[i]\) 表示以 \(i\) 结尾的方案数。

要满足 \(a[i]\&a[i+1]=0\) ，那么 \(f[i]\) 是从上一次结果中 \(i\&j=0\) 的结果转移过来的。

而 \(i\&j=0\) 又可等价于 \(i\&(\sim j)=i\) ，所以每次要对上一次的结果的下标进行取反。

所以，求 \(f[i]\) ，就相当于求以 \(i\) 为子集的上一次求出的 \(f\) 的高维前缀和。

对于 \(2\) 的处理：就把下标为 \(c[i]\) 的倍数的 \(f[i](i\%c[i]=0)\) 的取零即可。

code

#include<bits/stdc++.h>

const int maxn=55,maxm=15,mod=1e9;

namespace IO
{
    char buf[1<<15],*fs,*ft;
    inline char getc() { return (ft==fs&&(ft=(fs=buf)+fread(buf,1,1<<15,stdin),ft==fs))?0:*fs++; }
    template<typename T>inline void read(T &x)
    {
        x=0;
        T f=1, ch=getchar();
        while (!isdigit(ch) && ch^'-') ch=getchar();
        if (ch=='-') f=-1, ch=getchar();
        while (isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48), ch=getchar();
        x*=f;
    }

    char Out[1<<24],*fe=Out;
    inline void flush() { fwrite(Out,1,fe-Out,stdout); fe=Out; }
    template<typename T>inline void write(T x,char str)
    {
        if (!x) *fe++=48;
        if (x<0) *fe++='-', x=-x;
        T num=0, ch[20];
        while (x) ch[++num]=x%10+48, x/=10;
        while (num) *fe++=ch[num--];
        *fe++=str;
    }
}

using IO::read;
using IO::write;

int c[maxn],f[1<<maxm];
int main()
{
    int T;read(T);
    while (T--)
    {
        int n,m;read(n);read(m);
        for (int i=1; i<=n; ++i) read(c[i]);
        for (int i=0; i< (1<<m); ++i) f[i]=(i==0);//赋初值，相当于 f[0]=1, f[1~(1<<m)]=0
        for (int i=1; i<=n; ++i)
        {
            for (int j=0; j< (1<<m); j+=2) std::swap(f[j],f[j^((1<<m)-1)]);//把下标互补的值互换，下标之异或和等于 (1<<m)-1
            for (int j=0; j< m; ++j)
                for (int k=0; k< (1<<m); ++k)
                    if (~k&(1<<j)) f[k]=(f[k]+f[k|(1<<j)])%mod;//高维前缀和
            for (int j=0; j< (1<<m); j+=c[i]) f[j]=0;//把下标为 c[i] 的倍数的清零
        }
        int ans=0;
        for (int i=0; i< (1<<m); ++i) ans=(ans+f[i])%mod;//方案数叠加取模
        write(ans,'\n');
    }
    IO::flush();
    return 0;
}

转载于:https://www.cnblogs.com/G-hsm/p/Spoj_TLE.html