bzoj 3209: 花神的数论题 数位dp

        网上题解都是高大上的组合数学。。然而我等蒟蒻显然只会简单粗暴的数位dp，即：

       令f[i][j][k]表示在前i位中，状态为j，有k个1的方案数。其中j=0表示前i为<n，j=1表示>=n。方便起见可以令n++，然后转移就很方便了，枚举一下01即可。

AC代码如下：

#include<iostream>#include<cstdio>#include<cstring>#define mod 10000007#define ll long longusing namespace std;ll n,dp[105][2][105];int ksm(int x,ll y){int t=1; for (; y; y>>=1,x=(ll)x*x%mod) if (y&1) t=(ll)t*x%mod;return t;}int main(){scanf("%lld",&n); n++; int i,j,k,x;dp[1][1][0]=1;for (i=2; n; n>>=1,i++)for (j=0; j<2; j++)for (k=0; k<=50; k++)for (x=0; x<=max(j,(int)(n&1)) && x<=k; x++)dp[i][j][k]+=dp[i-1][j|(x<(n&1))][k-x];k=i-1; int ans=1;for (i=1; i<=50; i++) ans=(ll)ans*ksm(i,dp[k][0][i])%mod;printf("%d\n",ans);return 0;}

by lych

2016.4.21