hdu 5106
真的好累,就不解释代码了。
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
long long int dp[1002][1002],dpl[1002][1002];//dpl[i][j]用来计算前i位j个1的个数,dp[i][j]用来计算前i位,有j个1的数的个数。
int bit[1003],n,c[1001];
char r[1003];
#define mod 1000000007
long long int solve()
{
long long int answer=0,s=0;
int ss=0;
memset(dp,0,sizeof(dp));
memset(bit,0,sizeof(bit));
memset(dpl,0,sizeof(dpl));
int len,k;
len=strlen(r);
int i,j;
for(i=1;i<=len;i++)
{
bit[i]=r[i-1]-'0';
}
for(i=1;i<=len;i++)
{
for(j=0;j<bit[i];j++)
{
dp[i][s+j]+=(ss*2)%mod;
dpl[i][s+j]++;
}
s+=bit[i];
ss=ss*2+bit[i];
ss%=mod;
for(j=0;j<=n;j++)
{
dp[i+1][j]+=dp[i][j]*2%mod;
dpl[i+1][j]+=dpl[i][j];
dpl[i+1][j]%=mod;
dpl[i+1][j+1]+=dpl[i][j];
dpl[i+1][j+1]%=mod;
if(dp[i][j]==0&&j!=0)
continue;
dp[i+1][j+1]+=(dp[i][j]*2+dpl[i][j])%mod;
dp[i+1][j+1]%=mod; //状态转移,其中相当于将数向前移动一位,再加上dpl数组中相应的个数。
//printf("dp[%d][%d]=%d\n",i,j,dp[i][j]);
}
}
return dp[len][n];
}
int main () {
//freopen("d:\\in.txt","r",stdin);
while(~scanf("%d%s",&n,r))
{
int ans=solve();
printf("%d\n",ans);
}
return 0;
}
hdu 5435
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
long long int dp[100005][16],bit[100005];
char a[100005],b[100005];
#define mod 1000000007
long long int solve(char *str,int c)
{
long long int ans=0,s=0;
int len;
len=strlen(str);
memset(dp,0,sizeof(dp));
memset(bit,0,sizeof(bit));
for(int i=1;i<=len;i++)
{
bit[i]=str[i-1]-'0';
//printf("%d\n",bit[i]);
}
int j,k,f;
for(j=1;j<=len;j++)
{
for(k=0;k<bit[j];k++)
{
dp[j][s^k]++;
}
s^=bit[j];
for(int i=0;i<=9;i++)
for(f=0;f<=15;f++)
{
dp[j+1][i^f]+=dp[j][f];
dp[j+1][i^f]%=mod;
//printf("dp[%d][%d]=%d\n",j,f,dp[j][f]);
}
}
if(c!=0)
dp[len][s]++;
dp[len][0]--;
for(j=0;j<=15;j++)
{
ans+=dp[len][j]*j;
//printf("dp[%d][%d]=%d\n",len,j,dp[len][j]);
ans%=mod;
}
return ans;
}
int main () {
int cas;
//freopen("d:\\in.txt","r",stdin);
//freopen("d:\\inn.txt","w",stdout);
scanf("%d",&cas);
int y=1;
while(y<=cas)
{
scanf("%s%s",a,b);
long long int ss=solve(b,1)-solve(a,0);
printf("Case #%d: %d\n",y,(ss+mod)%mod);
//printf("%lld\n",solve(b,1));
y++;
}
return 0;
}
这题又让我想起了cf上类似于数位dp的乱搞题,貌似难度还是递增的,数据范围不断变大。最后有10^100左右
cf 146A (n=50)
cf(95B)(n=10^100)
cf(145C) n=100000
题意大概都一样:
我第二种卡了好久,算了别人代码送上吧:
#include<bits/stdc++.h>
using namespace std;
#define inf 1000000007
#define maxn 100010
char s[maxn];
int main()
{
//freopen("d:\\in.txt", "r", stdin);
int i,j,k,n,m,a,b;
scanf("%s",s);
n=strlen(s);
if(n&1||s[0]-'0'>7)
{
n+=2;if(n&1)n--;
for(i=0;i<n/2;i++)printf("4");
for(i=0;i<n/2;i++)printf("7");
puts("");
return 0;
}
bool flag=0,ok=0;m=n/2;int c4=0,c7=0,last=-1,t=-1;
for(i=0;i<n;i++)
{
if(flag)break;
if(s[i]<'4')
{
if(c4<m)
{
c4++;s[i]='4';flag=1;
}
else {s[i]='7';c7++;flag=1;}
}
else if(s[i]=='4')
{
if(c4<m)
{
c4++;last=i;
}
else {s[i]='7',c7++,flag=1;}
}
else if(s[i]<'7')
{
if(c7<m)
{
c7++;s[i]='7';flag=1;
}
else
{
if(t==-1){ok=1;break;}
s[t]='7',s[t+1]='4',i=t+1,flag=1;
}
}
else if(s[i]=='7')
{
if(c7<m)
{
c7++;if(i&&s[i-1]=='4')t=i-1;
}
else
{
if(t!=-1)
{
s[t]='7',s[t+1]='4',i=t+1,flag=1;
}
else{ok=1;break;}
}
}
else
{
if(c7>=m)
{
if(t==-1){ok=1;break;}
s[t]='7',s[t+1]='4',i=t+1,flag=1;
}
else
{
if(last!=-1)
{
if(last<i-1)
{j=last;s[j]='7',s[j+1]='4';i=j+1;flag=1;}
else{j=last;s[j]='7';s[i]='4';c7++;flag=1;}
}
else{ok=1;break;}
}
}
}
if(ok)
{
n+=2;
for(i=0;i<n/2;i++)
printf("4");
for(i=0;i<n/2;i++)
printf("7");
puts("");
return 0;
}
c4=0,c7=0;
for(j=0;j<i;j++)
printf("%c",s[j]),c4+=(s[j]=='4'),c7+=(s[j]=='7');
for(i=c4+1;i<=m;i++)
printf("4");
for(i=c7+1;i<=m;i++)
printf("7");
puts("");
return 0;
}
其实也算不上严格意义的数位dp,只是都要一位一位来进行思考。