hdu6125 Free from square 状压+分组

http://acm.hdu.edu.cn/showproblem.php?pid=6125

题意：从1～n里选不超过k个数使得选的数的成绩不能被平方数（不包括1）整除

想到了从因子角度考虑，不能被非1平方数整数，也就是选的数里没有相同的质因子，会想到按质因子分组背包，这样能够保证每一个质因子只能选一个，但有些数不单单只有一个质因子，也就是可以被分到多个组，这样显然是不行的，因为每个数只能被选一次

其实能被分到多组的质因子并没有多少个:{2,3,5,7,11,13,17,19} 去掉这八个小于根号n的质因子，剩下的每个数都可以被分到唯一的一组（大于根号n的两个质因子相乘肯定大于n，所以不会有分到多组的情况）

这样对这八个质因子用状压来实现不重复，而分组实现了像23，29这些因子的不重复（比如23和46是不能同时选的）

#include<bits/stdc++.h>
#define eps 1e-9
#define PI 3.141592653589793
#define bsize 256
#define MEM(a) memset(a,0,sizeof(a))
typedef long long ll;
using namespace std;
long long dp[2][(1<<8)+10][505];
int state[505],f[505];
vector<int>w[505];
const int mo=1e9+7;
int a[10]={2,3,5,7,11,13,17,19};
void init()
{
    int i,j;
    for(i=1;i<=500;i++)
    {
        f[i]=i;
        for(j=0;j<8;j++)
        {
            if(i%(a[j]*a[j])==0)
            {
                state[i]=-1;
                break;
            }
            else if(i%a[j]==0)
            f[i]/=a[j],state[i]|=(1<<j);
        }
    }
    for(i=1;i<=500;i++)
    {
        if(state[i]==-1)
            continue;
        if(f[i]==1)
            w[i].push_back(i);
        else
            w[f[i]].push_back(i);
    }
    return ;
}
int main()
{
    int T,i,j,n,k,z,x;
    cin>>T;
    init();
    while(T--)
    {
        memset(dp,0,sizeof(dp));
        dp[0][0][0]=1;
        cin>>n>>k;
        int now=0,pre=1;
        for(i=1;i<=n;i++)
        {
            if(w[i].size()==0)
            continue;
            swap(now,pre);
            for(j=0;j<(1<<8);j++)
            {
                for(z=0;z<=k;z++)
                dp[now][j][z]=dp[pre][j][z];
            }
            for(x=0;x<w[i].size();x++)
            {
                int nn=w[i][x];
                if(nn>n)
                    continue;
                for(j=0;j<(1<<8);j++)
                {
                    if(state[nn]&j)
                    continue;
                    for(z=1;z<=k;z++)
                    {
                        dp[now][state[nn]|j][z]=(dp[now][state[nn]|j][z]+dp[pre][j][z-1])%mo;
                    }
                }
            }
        }

        long long ans=0;
        for(i=0;i<(1<<8);i++)
        {
            for(j=1;j<=k;j++)
            ans=(ans+dp[now][i][j])%mo;
        }
        cout<<ans<<endl;
    }
    return 0;
 }