首页 > 其他分享> 文章详细

一本通1590恨 7 不成妻

2019-02-11 21:44:25 阅读：398 来源： 互联网

标签：tmp 10 ll 一本 1590 ans Mod 不成 he

1590：恨 7 不成妻

时间限制: 1000 ms 内存限制: 524288 KB

【题目描述】

单身！

依然单身！

吉哥依然单身！

DS 级码农吉哥依然单身！

所以，他平生最恨情人节，不管是 214 还是 77，他都讨厌！

吉哥观察了 214 和 77 这两个数，发现：

2+1+4=7

7+7=7×2

77=7×11

最终，他发现原来这一切归根到底都是因为和 7 有关！所以，他现在甚至讨厌一切和 7 有关的数！

什么样的数和 7 有关呢？如果一个整数符合下面三个条件之一，那么我们就说这个整数和 7 有关：

整数中某一位是 7；

整数的每一位加起来的和是 7 的整数倍；

这个整数是 7 的整数倍。

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

【输入】

输入数据的第一行是测试数据组数 T，然后接下来的 T 行表示 T 组测试数据。

每组数据在一行内包含两个正整数 L,R。

【输出】

对于每组数据，请计算 [L,R] 中和 7 无关的数字的平方和，并将结果对 10^9+7 取模后输出。

【输入样例】

【输出样例】

236
221
0

【提示】

数据范围与提示：

对于全部数据，1≤T≤50,1≤L≤R≤10^18 。

sol：如果只是统计个数就简单多了，然并卵是统计平方和，于是就愉快地看题解去了

（NOOOOOOOOOOOO）

计算与7无关的数的和当前这位的贡献就是（要填的数i）*（10^位数），设这个贡献为A

下一个状态是tmp，这个状态的平方和+=（A+tmp[与7无关的数的和]）², 就是A²+2*A*tmp[与7无关的数的和]+{（tmp[与7无关的数的和]²）=（tmp[与7无关的数的平方和]）}

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
inline ll read()
{
    ll s=0;
    bool f=0;
    char ch=' ';
    while(!isdigit(ch))
    {
        f|=(ch=='-');
        ch=getchar();
    }
    while(isdigit(ch))
    {
        s=(s<<3)+(s<<1)+(ch^48);
        ch=getchar();
    }
    return (f)?(-s):(s);
}
#define R(x) x=read()
inline void write(ll x)
{
    if(x<0)
    {
        putchar('-');
        x=-x;
    }
    if(x<10)
    {
        putchar(x+'0');
        return;
    }
    write(x/10);
    putchar((x%10)+'0');
    return;
}
inline void writeln(ll x)
{
    write(x);
    putchar('\n');
    return;
}
#define W(x) write(x),putchar(' ')
#define Wl(x) writeln(x)
const ll Mod=1000000007;
int T,Num[20];
ll Bin[20];
struct Record
{
    ll Ges,Shuz_he,Pingf_he;
    bool Vis;
}dp[20][10][10][2][2];
inline Record dfs(int Weiz,int Pre1,int Pre2,bool Shangj,bool Qiand0)
{
    if(dp[Weiz][Pre1][Pre2][Shangj][Qiand0].Vis) return dp[Weiz][Pre1][Pre2][Shangj][Qiand0];
    if(Weiz==0)
    {
        return dp[Weiz][Pre1][Pre2][Shangj][Qiand0]=(Record){(Pre1&&Pre2)?1:0,0,0,1};
    }
    Record ans=(Record){0,0,0,1};
    int i,Up=(Shangj)?(Num[Weiz]):(9);
    for(i=0;i<=Up;i++) if(i!=7)
    {
        bool Bo1=(Shangj&&i==Up),Bo2=(Qiand0&&i==0);
        Record tmp=dfs(Weiz-1,(Pre1+i)%7,(Pre2*10+i)%7,Bo1,Bo2);
        ll A=i*Bin[Weiz]%Mod;
        ans.Ges+=tmp.Ges;
        ans.Ges-=(ans.Ges>=Mod)?Mod:0;
        ans.Shuz_he+=(tmp.Shuz_he+A*tmp.Ges%Mod)%Mod;
        ans.Shuz_he-=(ans.Shuz_he>=Mod)?Mod:0;
        ans.Pingf_he+=(tmp.Pingf_he+2*A*tmp.Shuz_he%Mod+A*A%Mod*tmp.Ges%Mod)%Mod;
        ans.Pingf_he-=(ans.Pingf_he>=Mod)?Mod:0;
    }
    return (dp[Weiz][Pre1][Pre2][Shangj][Qiand0]=ans);
}
inline ll Solve(ll n)
{
    *Num=0;
    while(n)
    {
        Num[++*Num]=n%10;
        n/=10;
    }
    int i,j,k,ii,jj;
    for(i=0;i<20;i++) for(j=0;j<10;j++) for(k=0;k<10;k++)
    for(ii=0;ii<2;ii++) for(jj=0;jj<2;jj++) dp[i][j][k][ii][jj].Vis=0;
    return dfs(*Num,0,0,1,1).Pingf_he%Mod;
}
int main()
{
    int i;
    Bin[1]=1;
    for(i=2;i<=20;i++) Bin[i]=(Bin[i-1]*10)%Mod;
    ll l,r;
    R(T);
    while(T--)
    {
        R(l); R(r);
        Wl(((Solve(r)-Solve(l-1))%Mod+Mod)%Mod);
    }
    return 0;
}
/*
input
3
1 9
10 11
17 17
output
236
221
0
*/

View Code

标签：tmp,10,ll,一本,1590,ans,Mod,不成,he
来源： https://www.cnblogs.com/gaojunonly1/p/10363255.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9