UOJ 129/BZOJ 4197 寿司晚宴 状压DP

//By SiriusRen
#include <cstdio>
#include <algorithm>
using namespace std;
const int N=;
struct Node{int a,p;}node[N];
bool operator<(Node a,Node b){return a.p<b.p;}
char pr[]={,,,,,,,};
int n,p,f[N][N],dp[N][N][],ans;
void div(int x){
    int cpy=x;
    for(int i=;i<;i++)
        if(!(x%pr[i])){
            while(x%pr[i]==)x/=pr[i];
            node[cpy].a|=<<i;
        }
    node[cpy].p=x;
}
int main(){
    scanf("%d%d",&n,&p);
    for(int i=;i<=n;i++)div(i);
    sort(node+,node++n),f[][]=;
    for(int i=;i<=n;i++){
        if(node[i].p==||node[i].p!=node[i-].p)
            for(int j=;j<;j++)
                for(int k=;k<;k++)
                    dp[j][k][]=dp[j][k][]=f[j][k];
        for(int j=;~j;j--)
            for(int k=;~k;k--){
                if(!(node[i].a&k))(dp[j|node[i].a][k][]+=dp[j][k][])%=p;
                if(!(node[i].a&j))(dp[j][k|node[i].a][]+=dp[j][k][])%=p;
            }
        if(node[i].p==||node[i].p!=node[i+].p)
            for(int j=;~j;j--)
                for(int k=;~k;k--)
                    f[j][k]=((dp[j][k][]+dp[j][k][]-f[j][k])%p+p)%p;
    }
    for(int i=;i<;i++)for(int j=;j<;j++)if(!(i&j))(ans+=f[i][j])%=p;
    printf("%d\n",ans);
}