NOI 2015 寿司晚宴 (状压DP+分组背包)

题目大意：两个人从2~n中随意取几个数(不取也算作一种方案)，被一个人取过的数不能被另一个人再取。两个人合法的取法是，其中一个人取的任何数必须与另一个人取的每一个数都互质，求所有合法的方案数

(数据范围毕竟很小，乍一看也不是啥打表找规律的题)

和我之前做过的一道题很类似hdu 6125，但这道题由于题面看起来很玄学，所以正解更难想

但还是 状压DP+分组背包 的套路

因为500以内的任何一个数，只会有一个大于19的质因子，所以对2 3 5 7 11 13 17 19这8个质数进行状压，然后每个数都质因数分解，把小于等于19的质因子存入状态，剩下的因子分组背包搞搞就行了，注意如果剩下的因子是1要单独算一组，否则会出大事情，比如2和3并不是同一组的，如果再来一个4，和2是同一组的，转移就会出错

具体DP的实现呢，定义是第一个人取了状态为s1的数，第二个人取了状态为s2的数

分组背包要把同一组的东西放到连续的一段序列上

对于这道题而言，如果某个人取了某一组的任何一个，那么这一组的其它物品也只能被这个人取/不取

所以额外定义两个状态f1,f2，含义和dp的意义是一样的，只不过在同一组内是f1和f2这两个状态自己和自己转移，然后把答案贡献给dp，即这一组对整体的贡献，然后把dp重新赋给f1,f2，再进行下一组背包

方程   

由于f1,f2为了下一层转移，都被加了一次dp值，所以最后要减掉一个dp

即

而f1,f2转移也有技巧，常规的自己和自己转移为了避免传递性，要另外开一个数组进行转移。但因为这道题的转移方程都是位与|操作，具有递增性，所以倒序枚举就可以减少一些常数

 #include <queue>
 #include <cstdio>
 #include <cstring>
 #include <algorithm>
 #define ui unsigned int
 #define ll long long
 #define il inline
 #define N (1<<8)+3
 #define inf 0x3f3f3f3f
 using namespace std;

 int n;
 ll p;
 ll f1[N][N],f2[N][N],dp[N][N];
 int pr[]={,,,,,,,};
 struct node{
     int w,f;
     friend bool operator < (const node &a,const node &b){
         if(a.w!=b.w) return a.w<b.w;
         else return a.f<b.f;
     }
 }s[N];
 void get_son()
 {
     for(int i=;i<=n;i++)
     {
         int x=i;
         for(int j=;j<;j++)
         {
             if(x%pr[j]==) x/=pr[j],s[i].f|=(<<j);
             while(x%pr[j]==) x/=pr[j];
         }s[i].w=x;
     }
 }

 int main()
 {
     scanf("%d%lld",&n,&p);
     get_son();
     sort(s+,s+n+);
     f1[][]=f2[][]=dp[][]=;
     int m=(<<)-;
     for(int i=;i<=n;i++)
     {
         for(int s1=m;s1>=;s1--)
             for(int s2=m;s2>=;s2--){
                 if(!((s1|s[i].f)&s2)) f1[s1|s[i].f][s2]=(f1[s1|s[i].f][s2]+f1[s1][s2])%p;
                 if(!(s1&(s2|s[i].f))) f2[s1][s2|s[i].f]=(f2[s1][s2|s[i].f]+f2[s1][s2])%p;}
         if(s[i].w==||s[i+].w!=s[i].w)
             for(int s1=m;s1>=;s1--)
                 for(int s2=m;s2>=;s2--)
                     f1[s1][s2]=f2[s1][s2]=dp[s1][s2]=(f1[s1][s2]+f2[s1][s2]-dp[s1][s2]+p)%p;
     }
     ll ans=;
     for(int s1=m;s1>=;s1--)
         for(int s2=m;s2>=;s2--)
             ans+=dp[s1][s2],ans%=p;
     printf("%lld\n",ans);
     return ;
 }