hdu 4906 状压dp

题意：给出n,k<=20,L<=10^9 问能不能通过不大于L的数组成一个含n个数的序列使得能从这个序列中挑出几个数的和为k,求满足条件序列的个数

题解：看到k<20可以用sum的状态做状压dp的状态

状态转移方程next=(1<<(p-1))|j|((j<<p)&SIZE); dp[next]=(dp[next]+v)%mod;

#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll __int64
using namespace std;
const ll mod=1e9+7;
ll dp[1<<21],v;
int main()
{
    int T,n,k,l,MIN,SIZE,i,j,p;
    scanf("%d",&T);
    while(T--){
        scanf("%d%d%d",&n,&k,&l);
        memset(dp,0,sizeof(dp));
        dp[0]=1;//赋初值不然没法启动
        SIZE=(1<<k)-1;//最大状态
        ll MIN=min(l,k);//因为大于k的值是没用的可以在后面直接加上去
        for(i=1;i<=n;i++){//因为要加n个数这里枚举次数
            for(j=SIZE;j>=0;j--){//倒着枚举状态是类似背包一样防止重复加
                v=dp[j];
                if(v==0)continue;
                for(p=1;p<=MIN;p++){
                    int next=(1<<(p-1))|j|((j<<p)&SIZE);//加上数p后的状态1.p本身2.j本身3.j的每一位加上p后的状态
                    dp[next]=(dp[next]+v)%mod;
                }
                if(k<l)dp[j]=(dp[j]+(ll)v*(l-k))%mod;//剩下部分直接加上去不知道为什么看有些blog没限制条件答案也是对的
            }//是数据太弱吗
        }
        ll ans=0;
        for(i=0;i<=SIZE;i++)
            if(i&(1<<(k-1))) (ans+=dp[i])%=mod;
        printf("%I64d\n",ans);
    }
    return 0;
}