hdu 5106 hdu 5435 数位dp

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,只是都要一位一位来进行思考。