[NOI2015]寿司晚宴

题面

题意

给出n，在2~n中选择一些数，分成两组，使不存在两个数a，b，分别属于这两组数,且gcd(a,b)=1,问有几种分法

做法

首先题意可以转化为求两组数，使其两两之间不存在相同的质因数，因此可以考虑选择一些质因数分成两组，但是n<=500，其中有约100个质因数，无法转压，但是我们可以发现所有小于等于500的数的质因子中，大于19的质数至多只有一个，而19以内只有八个质数，因此我们可以对这8个质数进行状压，对大质数单独考虑。
dp时记录f1,f2两个数组，分别表示将此时考虑的大质数不放在B组和不放在A组的贡献（因此之后合并时要去重，即dp=f1+f2-dp，dp指未考虑这个大质数时的贡献），然后每次在考虑某个大质数（或是不包含大质数的数）前，将dp的值赋给f1,f2,考虑完所有包含此大质数的数后，再用上述方法将f1,f2的值计入dp即可。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
#define N 510
#define MN 255
#define ZY 300
using namespace std;

ll n,M,ans,dp[ZY][ZY],f1[ZY][ZY],f2[ZY][ZY],zs[10]= {2,3,5,7,11,13,17,19};
bool flag;
struct Num
{
    ll zt,big;
    void in(ll u)
    {
        ll i,j;
        for(i=0; i<8; i++)
        {
            if(u%zs[i]) continue;
            zt|=(1 << i);
            for(; u%zs[i]==0; u/=zs[i]);
        }
        big=u;
    }
    bool operator < (const Num &u) const
    {
        return big<u.big;
    }
} num[N];

int main()
{
    ll i,j,k;
    cin>>n>>M;
    for(i=2; i<=n; i++) num[i].in(i);
    sort(num+2,num+n+1);
    dp[0][0]=1;
    for(i=2; i<=n; i++)
    {
        if(i==2 || num[i].big==1 || num[i].big!=num[i-1].big)
        {
            memcpy(f1,dp,sizeof(dp));
            memcpy(f2,dp,sizeof(dp));
        }
        for(j=MN;j>=0;j--)
        {
            for(k=MN;k>=0;k--)
            {
                if(j&k) continue;
                if((num[i].zt&k)==0) f1[j|num[i].zt][k]+=f1[j][k],f1[j|num[i].zt][k]%=M;
                if((num[i].zt&j)==0) f2[j][k|num[i].zt]+=f2[j][k],f2[j][k|num[i].zt]%=M;
            }
        }
        if(i==n || num[i].big!=num[i+1].big || num[i].big==1)
        {
            for(j=MN;j>=0;j--)
            {
                for(k=MN;k>=0;k--)
                {
                    if(j&k) continue;
                    dp[j][k]=f1[j][k]+f2[j][k]-dp[j][k]+M;
                    dp[j][k]%=M;
                }
            }
        }
    }
    for(i=MN;i>=0;i--)
    {
        for(j=MN;j>=0;j--)
        {
            if(i&j) continue;
            ans+=dp[i][j];
            ans%=M;
        }
    }
    cout<<ans;
}