【Luogu】P2150寿司晚宴（状压DP）

　　题目链接

　　反正……我是没什么想法了，全程看题解

　　（或者说自己想了半天错解）

　　因为大于根n的质数最多只会在一个数里出现一种，所以可以把数拆成两部分：小数的二进制集合和大数。

　　然后把大数一样的放到一起DP，设s[i][j]表示第一个集合为i第二个为j的方案数，f[i][j][k]是第一个集合为i第二个集合为j，当前数放进第k个集合里……

　　转移当然比我那个sb的二维状态好转移啦

　　

#include<cstdio>
#include<cctype>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#include<cmath>
#define maxs 11
#define maxn 600
using namespace std;

inline long long read(){
    long long num=,f=;
    char ch=getchar();
    while(!isdigit(ch)){
        if(ch=='-')    f=-;
        ch=getchar();
    }
    while(isdigit(ch)){
        num=num*+ch-'';
        ch=getchar();
    }
    return num*f;
}

int prime[]={,,,,,,,,,,,,,,,};
int s[maxn][maxn],f[maxn][maxn][];

struct Number{
    int state,p;
    bool operator <(const Number a)const{
        return p<a.p;
    }
}q[maxn];

int main(){
    int n=read(),mod=read();
    for(int i=;i<=n;++i){
        int ret=i;
        q[i-].state=;
        for(int j=;j<;++j){
            int now=prime[j];
            if(ret%now)    continue;
            q[i-].state|=(<<(j-));
            while(ret%now==)    ret/=now;
        }
        q[i-].p=ret;
    }
    sort(q+,q+n);
    s[][]=;
    for(int i=;i<n;++i){
        Number now=q[i];
        if(now.p==||now.p!=q[i-].p)
            for(int S1=;S1<=;++S1)
                for(int S2=;S2<=;++S2)    f[S1][S2][]=f[S1][S2][]=s[S1][S2];

        for(int S1=;S1>=;--S1)
            for(int S2=;S2>=;--S2){
                if((now.state&S2)==)    f[now.state|S1][S2][]=(1LL*f[now.state|S1][S2][]+f[S1][S2][])%mod;
                if((now.state&S1)==)    f[S1][now.state|S2][]=(1LL*f[S1][now.state|S2][]+f[S1][S2][])%mod;
            }

        if(now.p==||now.p!=q[i+].p)
            for(int S1=;S1<=;++S1)
                for(int S2=;S2<=;++S2)    s[S1][S2]=((1LL*f[S1][S2][]+f[S1][S2][])%mod-s[S1][S2]+mod)%mod;
    }
    int ans=;
    for(int i=;i<=;++i)
        for(int j=;j<=;++j)
            if((i&j)==)    ans=(ans+s[i][j])%mod;
    printf("%d\n",ans);
    return ;
}