HDU3555数位dp

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#define N 100010
#define LL long long
#define debug(x) cout<<#x<<":"<<(x)<<endl;
using namespace std;
int a[25];
LL dp[25][3];
int main()
{
    // input
    memset(dp, 0, sizeof(dp));
    // 分情况递推, dp[0][0] = 1 需要自己推导
    dp[0][0] = 1;
    for(int i=1; i<=22; i++)
    {
        dp[i][0] = dp[i-1][0] * 10 - dp[i-1][1];
        dp[i][1] = dp[i-1][0];
        dp[i][2] = dp[i-1][2] * 10 + dp[i-1][1];
    }
    LL n;
    int t;
    scanf("%d", &t);
    while( t-- )
    {
        scanf("%lld", &n);
        n++;  // 整个算法是没有计上n自身的, 故49 = 1
        LL len = 0, ans = 0;
        bool flag = false;
        while ( n )
        {
            a[++len] = n % 10;
            n /= 10;
        }
        a[len + 1] = 0;
        // 下面的for循环可以看看
        // http://wenku.baidu.com/link?url=gWlPRquCGXVsRbD8Modjq0iFkfTjwz6ZUm6zXmFYwOmzG5nNk7uqNULB9RNKg5s3UF8VAUuckeE5lZgZQnDOhnBrCzncFczX0J5maFKxb-3
        // 第10页开始
        for (int i=len; i; i--)
        {
            ans += dp[i-1][2] * a[i];
            if (flag)
                ans += dp[i-1][0] * a[i];
            if(!flag && a[i] > 4)
                ans += dp[i-1][1];
            if (a[i+1] == 4 && a[i] == 9) flag = true;
        }
        printf("%lld\n", ans);
    }
    return 0;
}