[NOI2015][BZOJ4197][状压DP]寿司晚宴

Description
为了庆祝 NOI 的成功开幕，主办方为大家准备了一场寿司晚宴。小 G 和小 W 作为参加 NOI 的选手，也被邀请参加了寿司晚宴。
在晚宴上,主办方为大家提供了 n−1 种不同的寿司，编号 1,2,3,…,n−1，其中第 i 种寿司的美味度为 i+1 （即寿司的美味度为从 2 到 n）。
现在小 G 和小 W 希望每人选一些寿司种类来品尝，他们规定一种品尝方案为不和谐的当且仅当：小 G 品尝的寿司种类中存在一种美味度为 x 的寿司，小 W 品尝的寿司中存在一种美味度为 y 的寿司，而 x 与 y 不互质。
现在小 G 和小 W 希望统计一共有多少种和谐的品尝寿司的方案（对给定的正整数 p 取模）。注意一个人可以不吃任何寿司。
Input
输入文件的第 1 行包含 2 个正整数 n,p，中间用单个空格隔开，表示共有 n 种寿司，最终和谐的方案数要对 p 取模。
Output
输出一行包含 1 个整数，表示所求的方案模 p 的结果。
Sample Input
3 10000
Sample Output
9
HINT
2≤n≤500
0

第
一眼状压无疑。但是500以内的质数有80+，2^80是无法完成的。考虑到一个性质，n以内的数最多只含有一个大于根号n的质因子。所以考虑将质数分类，其中小于根号500的质数只有八个，将其压成2^8二进制数即可，其余的特殊保存在结构体中。接下来就可以根据speical分类来dp了,f[i][j]表示第一个人选的i状态第二个人选的j状态的方案数。因为同一个special的素数之间不能互相转移，所以对于每一类spe开始dp时，我们可以生成一个g[1/2][i][j]代表当前1/2做决策，1状态i,2状态j的方案数，预处理就应memcpy(g[1],f,sizeof(f)); memcpy(g[2],f,sizeof(f));
值得注意的是，类似于01背包，状态应从255倒序枚举，表示不能从当前的spe转移过来。
最后一类Dp完之后，更新f[i][j]=g[1][i][j]+g[2][i][j]，但应减去一个f[i][j]，应为加的时候加多了

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=500+5;
int n,mod,f[301][301],g[3][301][301],ans;
struct data
{
    int spe,nor;
    bool operator<(const data&x)const
    {
        if (spe^x.spe) return spe<x.spe;
        return nor<x.nor;
    }
}Prime[maxn];
#define spe(x) Prime[x].spe
#define nor(x) Prime[x].nor
const int prime[8]={2,3,5,7,11,13,17,19};
int main()
{
//  freopen("dinner.in","r",stdin);
    scanf("%d%d",&n,&mod);
    for (int i=1;i<=n;i++)
    {
        int temp=i;
        for (int j=0;j<8;j++) if (!(temp%prime[j]))
        {
            nor(i)|=1<<j;
            while (!(temp%prime[j])) temp/=prime[j];
        }
        spe(i)=temp;
    }
    sort(Prime+2,Prime+n+1);
    f[0][0]=1;
    for (int i=2;i<=n;i++)
    {
        if (i==2|spe(i)==1||spe(i)!=spe(i-1))
        {
            memcpy(g[1],f,sizeof(f));
            memcpy(g[2],f,sizeof(f));
        }
        for (int j=255;j>=0;j--)
          for (int k=255;k>=0;k--)
          {
            if ((k&nor(i))==0) g[1][j|nor(i)][k]=(g[1][j|nor(i)][k]+g[1][j][k])%mod;
            if ((j&nor(i))==0) g[2][j][k|nor(i)]=(g[2][j][k|nor(i)]+g[2][j][k])%mod;
          }
        if (i==n||spe(i)==1||spe(i)!=spe(i+1))
        {
          for (int j=0;j<=255;j++)
            for (int k=0;k<=255;k++)
              f[j][k]=((g[1][j][k]+g[2][j][k]-f[j][k])%mod+mod)%mod;
        }
    }
    for (int i=0;i<=255;i++)
      for (int j=0;j<=255;j++)
        if (!(i&j)) ans=(ans+f[i][j])%mod;
    printf("%d",ans);
    return 0;
}