Day9 - J - 吉哥系列故事——恨7不成妻 HDU - 4507

单身!
　　依然单身！
　　吉哥依然单身！
　　DS级码农吉哥依然单身！
　　所以，他生平最恨情人节，不管是214还是77，他都讨厌！
　　
　　吉哥观察了214和77这两个数，发现：
　　2+1+4=7
　　7+7=7*2
　　77=7*11
　　最终，他发现原来这一切归根到底都是因为和7有关！所以，他现在甚至讨厌一切和7有关的数！

　　什么样的数和7有关呢？

　　如果一个整数符合下面3个条件之一，那么我们就说这个整数和7有关——
　　1、整数中某一位是7；
　　2、整数的每一位加起来的和是7的整数倍；
　　3、这个整数是7的整数倍；

　　现在问题来了：吉哥想知道在一定区间内和7无关的数字的平方和。

Input输入数据的第一行是case数T(1 <= T <= 50)，然后接下来的T行表示T个case;每个case在一行内包含两个正整数L, R(1 <= L <= R <= 10^18)。
Output请计算[L,R]中和7无关的数字的平方和，并将结果对10^9 + 7 求模后输出。Sample Input

Sample Output

236

221

0

思路：还是数位DP，套用模板，求的是平方和，记录其后数位满足条件的个数即可，比如233, 234, 235满足，则平方和为(200+33)^2+(200+34)^2+(200+35)^2 = 3*200^2 + 2*200*(33+34+35)+33^2+34^2+35^2,记录下每次的和与平方和维护即可

typedef long long LL;

typedef pair<LL, LL> PLL;

const int MOD = 1e9+;

struct Node {

    LL num, sum, ssum;

} dp[][][];

LL a[], p[];

Node dfs(int pos, int pre, int presum, bool limit) {

    if(pos == -) {

        Node ret;

        ret.num = (pre !=  && presum != );

        ret.sum = ret.ssum = ;

        return ret;

    }

    if(!limit && dp[pos][pre][presum].num != -) return dp[pos][pre][presum];

    int up = limit?a[pos]:;

    Node ans;

    ans.num = ans.sum = ans.ssum = ;

    for(int i = ; i <= up; ++i) {

        if(i == ) continue;

        Node nex = dfs(pos-, (pre+i)%, (presum*+i)%, limit&&i==a[pos]);

        if(nex.num == ) continue;

        ans.num = ((ans.num + nex.num) % MOD + MOD) % MOD;

        ans.sum = ((ans.sum + nex.sum + (p[pos]*i)%MOD*nex.num) % MOD + MOD) % MOD;

        ans.ssum = ((ans.ssum + nex.ssum + ((*p[pos]*i) % MOD)*nex.sum) % MOD + MOD) % MOD;

        ans.ssum = ((ans.ssum + (p[pos]*nex.num)%MOD * p[pos] % MOD *i*i) % MOD + MOD) % MOD;

    }

    if(!limit) dp[pos][pre][presum] = ans;

    return ans;

}

LL solve(LL n, LL m) {

    int pos = ;

    while(n) {

        a[pos++] = n % ;

        n /= ;

    }

    LL t1 = dfs(pos-, , , true).ssum;

    pos = ;

    while(m) {

        a[pos++] = m % ;

        m /= ;

    }

    t1 = dfs(pos-, , , true).ssum - t1;

    return (t1%MOD+MOD)%MOD;

}

void run_case() {

    LL n, m;

    cin >> n >> m;

    cout << solve(n-, m) << "\n";

}

int main() {

    ios::sync_with_stdio(false), cin.tie();

    int t;

    cin >> t;

    for(int i = ; i < ; ++i)

        for(int j = ; j < ; ++j)

            for(int k = ; k < ; ++k)

                dp[i][j][k].num = -;

    p[] = ;

    for(int i = ; i < ; ++i)

        p[i] = (p[i-]*)%MOD;

    while(t--)

        run_case();

    return ;

}