A serious math problem HDU - 5435 （数位dp+大数+计数问题）

思路：每一位取值只能是0到9,可以证明0到9中的数进行异或操作，不管多少个数进行异或，最后的值都不可能超过16，因为16转化为二进制是10000,已经是5位二进制数了，最坏情况是4位都是1,1111==15。所以我们可以利用这个关系，解决和太大，无法利用dp记忆化的问题。

也就是说求异或和的问题，转化成了一个，先求出最后答案为0-15的异或值有多少个的问题，只要计算出每种异或值有多少个，求总和只要乘起来就行了。

转化的思路已经完成了，下面就是考虑如何求出每种异或值有多少个的问题，这就转化成了数位dp计数问题，数位dp计数问题应该是比较常见的一种问题。

dp[i][j][k]表示前i位，异或值为j，k==0表示枚举到了第i位最高位，k==1表示还没有枚举到第i位的最高位。

dp的方式，选择按位枚举+递推。也就是说，当我们知道前i为的所有枚举情况的时候，去推前i+1为的情况。

dp[i+1][j^k][1]+=dp[i][j][0]  k>=0&&k<num[i+1]，表示第i为已经枚举到了最高位，而i+1位还没有到最高位，所以取值只能是k>=0&&k<num[i]+1,当k==num[i]时，表明i+1位已经枚举到了最高位，这个需要记录在dp[i+1][j^k][0]+=dp[i][j][0]上。

 

dp[i+1][j^k][1]+=dp[i][j][1] 表示第i位还没有枚举到最高位，那么i+1位可以随便取0到9的任意数，直接遍历一遍就ok。

这也是我第一次写非递归的数位dp，写起来比较痛苦，看着别人博客想了很久，最后想通了，也让我对数位dp有了新的认识。对于数位dp，但你发现需要记录的数太大，没办法开那么大的dp数组时，你需要做的就是去找关系，将其问题分解成小数据，可以开的dp数组，以前做过将模数传进来的题目，今天又遇见了这样记录每种数有多少个的问题，他们的大致思想都一样。一开始没想到，希望以后遇到困难，都想着先转化一下。

还有要注意的是：这个需要字符串读入，因为长度为1e5

代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const int maxn=1e5+7;
const int mod=1e9+7;
char s1[maxn],s2[maxn];
ll dp[maxn][16][2];
int bit[maxn];
ll solve(char *s)
{
    int n=strlen(s);
    for(int i=1;i<=n;i++) bit[i]=s[i-1]-'0';
    memset(dp,0,sizeof(dp));
    dp[0][0][0]=1;
    for(int i=0;i<n;i++)
    {
        int num=bit[i+1];
        for(int j=0;j<=15;j++)
        {
            dp[i][j][0]%=mod;
            dp[i][j][1]%=mod;
            if(dp[i][j][0]>0)//如果前i位已经枚举完了
            {
                dp[i+1][j^num][0]+=dp[i][j][0];//最后一位
                for(int k=0;k<num;k++)//当第i位取到最大值是，i+1位只能取0-num[i]
                {
                    dp[i+1][j^k][1]+=dp[i][j][0];
                }
            }
            if(dp[i][j][1]>0)//前i位还没有枚举到最高位的情况
            {
                for(int k=0;k<=9;k++)
                {
                    dp[i+1][k^j][1]+=dp[i][j][1];
                }
            }
        }
    }
    ll ans=0;
    for(int i=1;i<=15;i++)
    {
        ans+=(dp[n][i][0]+dp[n][i][1])*1ll*i%mod;
        ans%=mod;
    }
    ans%=mod;
    return ans;
}
int main()
{
    #ifndef ONLINE_JUDGE
        freopen("in.txt","r",stdin);
        freopen("out.txt","w",stdout);
    #endif
    int T;
    cin>>T;
    int Case=0;
    while(T--)
    {
        scanf("%s%s",s1,s2);
        ll ans=0;
        int n=strlen(s1);
        ans=s1[0]-'0';
        for(int i=1;i<n;i++)
        {
            ans^=s1[i]-'0';
        }
        printf("Case #%d: %lld\n",++Case,(solve(s2)-solve(s1)+ans+mod)%mod);
    }
    return 0;
}