spoj 10606 Balanced Numbers 数位dp

题目链接

一个数称为平衡数， 满足他各个数位里面的数， 奇数出现偶数次， 偶数出现奇数次， 求一个范围内的平衡数个数。

用三进制压缩， 一个数没有出现用0表示， 出现奇数次用1表示， 出现偶数次用2表示， 这样只需要开一个20*60000的数组。

#include<bits/stdc++.h>
using namespace std;
#define pb(x) push_back(x)
#define ll long long
#define mk(x, y) make_pair(x, y)
#define lson l, m, rt<<1
#define mem(a) memset(a, 0, sizeof(a))
#define rson m+1, r, rt<<1|1
#define mem1(a) memset(a, -1, sizeof(a))
#define mem2(a) memset(a, 0x3f, sizeof(a))
#define rep(i, a, n) for(int i = a; i<n; i++)
#define ull unsigned long long
typedef pair<int, int> pll;
const double PI = acos(-1.0);
const double eps = 1e-8;
const int mod = 1e9+7;
const int inf = 1061109567;
const int dir[][2] = { {-1, 0}, {1, 0}, {0, -1}, {0, 1} };
int digit[20], a[20];
ll dp[20][60000];
int judge(int num) {
    int cnt = 0;
    for(int i = 0; i<10; i++) {
        a[i] = num%3;
        num/=3;
    }
    for(int i = 0; i<10; i++) {
        if(i%2==0&&a[i]==2)
            return 0;
        if(i%2==1&&a[i]==1)
            return 0;
    }
    return 1;
}
int cal(int num, int tmp) {
    int cnt = 0;
    for(int i = 0; i<10; i++) {
        a[i] = num%3;
        num/=3;
    }
    a[tmp]++;
    if(a[tmp]==3)
        a[tmp]=1;
    for(int i = 9; i>=0; i--) {
        num = num*3+a[i];
    }
    return num;
}
ll dfs(int len, int num, int fp, bool first) {
    if(!len) {
        return judge(num);
    }
    if(!fp&&dp[len][num]!=-1) {
        return dp[len][num];
    }
    ll ret = 0;
    int maxx = fp?digit[len]:9;
    for (int i = 0; i<=maxx; i++) {
        ret += dfs(len-1, (first&&i==0)?0:cal(num, i), fp&&i==maxx, i==0&&first);
    }
    if(!fp)
        return dp[len][num] = ret;
    return ret;
}
ll cal(ll n) {
    int len = 0;
    while(n) {
        digit[++len] = n%10;
        n/=10;
    }
    return dfs(len, 0, 1, true);
}
int main()
{
    mem1(dp);
    int t;
    ll a, b;
    cin>>t;
    while(t--) {
        scanf("%lld%lld", &a, &b);    //I64d会超时......
        printf("%lld
", cal(b)-cal(a-1));
    }
}