hdu 吉哥系列故事——恨7不成妻 数位DP

最新推荐文章于 2020-09-24 15:44:42 发布
最新推荐文章于 2020-09-24 15:44:42 发布
阅读量714 收藏 1
点赞数
分类专栏： DP
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/ACMmaxx/article/details/19416717
版权

唉，以前就做过了，重做了一遍。居然花了4、5个小时，才做出来。各种惨不忍睹的错误o(╯□╰)o

其实，无论是求和，还是求平方和，都是较裸的数位DP。求和还是比较简单的，这里就谈谈怎么求平方和。

比如，我们要求ABCX*ABCX。由ABCX==ABC0 + X，整理有，ABCX*ABCX == ABC0*ABC0 + X*X + 2*ABC0*X

我们用dp[i][0/1][p][q]表示，从高位开始，已经填好了i个数字，当前是否（1/0）将这前i位置满，每位和除以7余p，前i位除以7余q，同时，前i位中，不能有7，这个状态下，有多少个数字。当然，sum[ ][ ][ ] []则表示这状态下，数字的和，sum2[ ][ ][ ][ ]则是平方和。

当我们要计算sum2[k][ ][ ][ ]的时候，按我们ABCX那个例子知，可以由dp[k-1][ ][ ][ ]和sum[k-1][ ][ ][ ]推出来。


//#pragma comment(linker, "/STACK:102400000,102400000")
#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
#include<cmath>
#include<cctype>
#include<string>
#include<algorithm>
#include<iostream>
#include<ctime>
#include<map>
#include<set>
using namespace std;
#define MP(x,y) make_pair((x),(y))
#define PB(x) push_back(x)
typedef __int64 LL;
//typedef unsigned __int64 ULL;
/* ****************** */
const int INF=1000111222;
const double INFF=1e200;
const double eps=1e-8;
const LL mod=1000000007;
const int NN=1100;
const int MM=10010;
/* ****************** */

int a[22];
LL pow10[22];
LL dp[22][2][8][8];
LL sum[22][2][8][8];
LL sum2[22][2][8][8];

void calc_cnt(int tol)
{
    memset(dp,0,sizeof(dp));
    int i,j,k,k1;
    int nj,nk;
    dp[0][1][0][0]=1;
    for(i=0;i<tol;i++)
    {
        for(j=0;j<7;j++)
            for(k=0;k<7;k++)
            {
                if(dp[i][0][j][k])
                {
                    for(k1=0;k1<10;k1++)
                    {
                        if(k1==7)
                            continue;
                        nj=(j+k1)%7;
                        nk=(k*10+k1)%7;
                        dp[i+1][0][nj][nk]+=dp[i][0][j][k];
                        if(dp[i+1][0][nj][nk]>=mod)
                            dp[i+1][0][nj][nk]-=mod;
                    }
                }
                if(dp[i][1][j][k])
                {
                    for(k1=0;k1<a[i+1];k1++)
                    {
                        if(k1==7)
                            continue;
                        nj=(j+k1)%7;
                        nk=(k*10+k1)%7;
                        dp[i+1][0][nj][nk]+=dp[i][1][j][k];
                        if(dp[i+1][0][nj][nk]>=mod)
                            dp[i+1][0][nj][nk]-=mod;
                    }
                    if(k1!=7)
                    {
                        nj=(j+k1)%7;
                        nk=(k*10+k1)%7;
                        dp[i+1][1][nj][nk]+=dp[i][1][j][k];
                        if(dp[i+1][1][nj][nk]>=mod)
                            dp[i+1][1][nj][nk]-=mod;
                    }
                }
            }
    }
/*
    puts("fuck!");
    for(i=0;i<=tol;i++)
    {
        for(j=0;j<7;j++)
        {
            for(k=0;k<7;k++)
            {
                cout<<dp[i][0][j][k]<<" ";
            }
            cout<<endl;
        }
        puts("-----");
    }
    puts("fuck!!!!!!!!!");
    */
}

void calc_sum(int tol)
{
    memset(sum,0,sizeof(sum));
    int i,j,k,k1;
    int nj,nk;
    LL t;
    for(i=0;i<tol;i++)
    {
        for(j=0;j<7;j++)
            for(k=0;k<7;k++)
            {
                for(k1=0;k1<10;k1++)
                {
                    if(k1==7)
                        continue;
                    nj=(j+k1)%7;
                    nk=(k*10+k1)%7;

                    t=(pow10[tol-1-i]*k1%mod)*dp[i][0][j][k];
                    t%=mod;
                    sum[i+1][0][nj][nk]+=t;
                    if(sum[i+1][0][nj][nk]>=mod)
                        sum[i+1][0][nj][nk]-=mod;

                    t=sum[i][0][j][k];
                    sum[i+1][0][nj][nk]+=t;
                    if(sum[i+1][0][nj][nk]>=mod)
                        sum[i+1][0][nj][nk]-=mod;

                }
                for(k1=0;k1<a[i+1];k1++)
                {
                    if(k1==7)
                        continue;
                    nj=(j+k1)%7;
                    nk=(k*10+k1)%7;

                    t=(pow10[tol-1-i]*k1%mod)*dp[i][1][j][k];
                    t%=mod;
                    sum[i+1][0][nj][nk]+=t;
                    if(sum[i+1][0][nj][nk]>=mod)
                        sum[i+1][0][nj][nk]-=mod;

                    t=sum[i][1][j][k];
                    sum[i+1][0][nj][nk]+=t;
                    if(sum[i+1][0][nj][nk]>=mod)
                        sum[i+1][0][nj][nk]-=mod;
                }
                if(k1!=7)
                {
                    nj=(j+k1)%7;
                    nk=(k*10+k1)%7;

                    t=(pow10[tol-1-i]*k1%mod)*dp[i][1][j][k];
                    t%=mod;
                    sum[i+1][1][nj][nk]+=t;
                    if(sum[i+1][1][nj][nk]>=mod)
                        sum[i+1][1][nj][nk]-=mod;

                    t=sum[i][1][j][k];
                    sum[i+1][1][nj][nk]+=t;
                    if(sum[i+1][1][nj][nk]>=mod)
                        sum[i+1][1][nj][nk]-=mod;
                }
            }
    }
/*
    puts("sum!");
    for(i=0;i<=tol;i++)
    {
        for(j=0;j<7;j++)
        {
            for(k=0;k<7;k++)
            {
                cout<<sum[i][1][j][k]<<" ";
            }
            cout<<endl;
        }
        puts("-----");
    }
    puts("sum!!!!!!!!!");
*/
}

void calc_sum2(int tol)
{
    memset(sum2,0,sizeof(sum2));
    int i,j,k,k1;
    int nj,nk;
    LL t;
    for(i=0;i<tol;i++)
    {
        for(j=0;j<7;j++)
            for(k=0;k<7;k++)
            {
                for(k1=0;k1<10;k1++)
                {
                    if(k1==7)
                        continue;
                    nj=(j+k1)%7;
                    nk=(k*10+k1)%7;

                    sum2[i+1][0][nj][nk]+=sum2[i][0][j][k];
                    if(sum2[i+1][0][nj][nk]>=mod)
                        sum2[i+1][0][nj][nk]-=mod;

                    t=pow10[tol-1-i]*k1%mod;
                    t*=sum[i][0][j][k]*2;
                    t%=mod;
                    sum2[i+1][0][nj][nk]+=t;
                    if(sum2[i+1][0][nj][nk]>=mod)
                        sum2[i+1][0][nj][nk]-=mod;

                    t=pow10[tol-1-i]*k1%mod;
                    t=t*t%mod;
                    t=t*dp[i][0][j][k];
                    t%=mod;
                    sum2[i+1][0][nj][nk]+=t;
                    if(sum2[i+1][0][nj][nk]>=mod)
                        sum2[i+1][0][nj][nk]-=mod;
                }
                for(k1=0;k1<a[i+1];k1++)
                {
                    if(k1==7)
                        continue;
                    nj=(j+k1)%7;
                    nk=(k*10+k1)%7;

                    sum2[i+1][0][nj][nk]+=sum2[i][1][j][k];
                    if(sum2[i+1][0][nj][nk]>=mod)
                        sum2[i+1][0][nj][nk]-=mod;

                    t=pow10[tol-1-i]*k1%mod;
                    t*=sum[i][1][j][k]*2;
                    t%=mod;
                    sum2[i+1][0][nj][nk]+=t;
                    if(sum2[i+1][0][nj][nk]>=mod)
                        sum2[i+1][0][nj][nk]-=mod;

                    t=pow10[tol-1-i]*k1%mod;
                    t=t*t%mod;
                    t=t*dp[i][1][j][k];
                    t%=mod;
                    sum2[i+1][0][nj][nk]+=t;
                    if(sum2[i+1][0][nj][nk]>=mod)
                        sum2[i+1][0][nj][nk]-=mod;
                }
                if(k1!=7)
                {
                    nj=(j+k1)%7;
                    nk=(k*10+k1)%7;

                    sum2[i+1][1][nj][nk]+=sum2[i][1][j][k];
                    if(sum2[i+1][1][nj][nk]>=mod)
                        sum2[i+1][1][nj][nk]-=mod;

                    t=pow10[tol-1-i]*k1%mod;
                    t*=sum[i][1][j][k]*2;
                    t%=mod;
                    sum2[i+1][1][nj][nk]+=t;
                    if(sum2[i+1][1][nj][nk]>=mod)
                        sum2[i+1][1][nj][nk]-=mod;

                    t=pow10[tol-1-i]*k1%mod;
                    t=t*t%mod;
                    t=t*dp[i][1][j][k];
                    t%=mod;
                    sum2[i+1][1][nj][nk]+=t;
                    if(sum2[i+1][1][nj][nk]>=mod)
                        sum2[i+1][1][nj][nk]-=mod;
                }
            }
    }
}

LL solve(LL x)
{
    int i,j,tol=0;
    LL ans=0;
    do
    {
        a[++tol]=x%10;
        x=x/10;
    }while(x);
    for(i=1;i<=tol/2;i++)
        swap(a[i],a[tol+1-i]);

  //  puts("fen");
  //  for(i=1;i<=tol;i++)
  //      printf("%d%c",a[i],i==tol?'\n':' ');

    pow10[0]=1;
    for(i=1;i<22;i++)
    {
        pow10[i]=pow10[i-1]*10;
        pow10[i]%=mod;
    }

    calc_cnt(tol);
    calc_sum(tol);
    calc_sum2(tol);

    for(i=1;i<7;i++)
        for(j=1;j<7;j++)
        {
            ans+=sum2[tol][0][i][j];
            if(ans>=mod)
                ans-=mod;
            ans+=sum2[tol][1][i][j];
            if(ans>=mod)
                ans-=mod;
        }
    return ans;
}

int main()
{
   // freopen("E:\\in.txt","r",stdin);
  //  freopen("E:\\no.txt","w",stdout);

    int cas;
    LL s,e,ans;
    scanf("%d",&cas);
    while(cas--)
    {
        cin>>s>>e;

     //   cout<<"s=="<<solve(s-1)<<endl;
    //    cout<<"e=="<<solve(e)<<endl;

        ans=solve(e)-solve(s-1);
        if(ans<0)
            ans+=mod;
        cout<<ans<<endl;
    }
    return 0;
}



ACMmaxx
关注
专栏目录
HDU 4507 吉哥系列故事——恨7不成妻 数位dp
每天都要进步哦~
06-04 364
题目链接：吉哥系列故事——恨7不成妻。此题和HDU 3652 B-number这道题很相似，如果题解看不懂，可以先看3652这道题。 代码如下： #include #include #include #include #include #define mod1 (m+7-(j*(bit[i]%7))%7)%7 #define mod2 (n+7-j%7)%7 using
hdu 4507 吉哥系列故事——恨7不成妻 (数位DP求和)
咸鱼有梦想
11-07 1597
吉哥系列故事——恨7不成妻 Time Limit: 1000/500 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others) Total Submission(s): 4951    Accepted Submission(s): 1620 Problem Description 　　单身! 　　依然单身！
HDU 4507 吉哥系列故事——恨7不成妻(数位DP)
smilestruggler_Fly To Higher
09-30 209
题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4507 题目大意：找出区间内不满足整数中某一位是7或整数的每一位加起来的和是7的整数倍或这个整数是7的整数倍的数的平方和 题目思路：如果只是单纯的求解这种数的个数，相信都嗤之以鼻，就是一个模板题，但是这题是求平方和，应该怎么办？首先我们可以发现，如果知道前面的数字是若干数的平方和，如1^...
Hdu-4507 吉哥系列故事——恨7不成妻(数位DP)
WZJRJ28
08-19 368
Description 　　单身! 　　依然单身！ 　　吉哥依然单身！ 　　DS级码农吉哥依然单身！ 　　所以，他生平最恨情人节，不管是214还是77，他都讨厌！ 　　 　　吉哥观察了214和77这两个数，发现： 　　2+1+4=7 　　7+7=7*2 　　77=7*11 　　最终，他发现原来这一切归根到底都是因为和7有关！所以，他现在甚至讨厌一切和7
HDU4507 吉哥系列故事——恨7不成妻（数位DP）
weixin_33962621的博客
01-08 125
题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4507 题意：如果一个整数符合下面3个条件之一，那么我们就说这个整数和7有关： 　　1、整数中某一位是7；　　2、整数的每一位加起来的和是7的整数倍；　　3、这个整数是7的整数倍；　　 现在吉哥想知道在一定区间内[L,R]和7无关的数字的平方和。 思路： #include &lt;io...
HDU 4507 吉哥系列故事——恨7不成妻 数位DP
_紫香妃的博客
09-24 194
文章目录题目描述 题目链接 题目描述 给出
HDU - 4507 吉哥系列故事——恨7不成妻(数位dp)
Falcon的博客
08-16 531
题目链接：点击查看 题目大意：给出闭区间[l,r]，求出区间内与7不沾边的数的平方和 与7有关的条件： 整数中某一位是7； 整数的每一位加起来的和是7的整数倍； 这个整数是7的整数倍； 题目分析：抛去平方和，如果改成统计有多少个符合条件的数字，那就是传统的数位dp了，只需维护一个dp[pos][sum1][sum2]即可，pos表示位数，sum1表示第二个条件的余数，sum2表示第三...
HDU - 4507 吉哥系列故事——恨7不成妻 数位dp
还我笑颜
08-11 286
题目链接：https://cn.vjudge.net/problem/HDU-4507 　　单身! 　　依然单身！ 　　吉哥依然单身！ 　　DS级码农吉哥依然单身！ 　　所以，他生平最恨情人节，不管是214还是77，他都讨厌！ 　　 　　吉哥观察了214和77这两个数，发现： 　　2+1+4=7 　　7+7=7*2 　　77=7*11 　　最终，他发现原来这一切归根到底都是因...
hdu 4507 吉哥系列故事——恨7不成妻
overload1997
10-14 348
吉哥系列故事——恨7不成妻 Time Limit: 1000/500 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others) Total Submission(s): 3493    Accepted Submission(s): 1091 Problem Description 　　单身! 　　依然单身
数位dp——某次队内分享任务
07-21
在介绍数位DP之前，我们先来看一个具体的例子，即题目【不要62 HDU-2089】。这个问题要求计算在给定区间[n, m]中不包含数字4或连续数字62的吉利数的数量。朴素的解决方案是预处理1到1e6的所有数，判断它们是否吉利，...
HDU - 4507 吉哥系列故事——恨7不成妻（数位DP求平方和）
Plus Ultra！
11-22 233
链接：HDU - 4507 吉哥系列故事——恨7不成妻 题意： 如果一个整数符合下面3个条件之一，那么我们就说这个整数和7有关 　　1、整数中某一位是7； 　　2、整数的每一位加起来的和是7的整数倍； 　　3、这个整数是7的整数倍； 现在问题来了：吉哥想知道在一定区间内和7无关的数字的平方和。 分析： 设数位为：dlen−1dlen−2⋯d2d1d0d_{len-1}d_{len-2}\cdo...
hdu 4507 吉哥系列故事——恨7不成妻 题解（数位DP（超难））
LightningUZ的博客
05-28 415
原题链接： HDU：点我QωQ 题意简述 定义一个正整数和777有关，满足以下三个之一： 某一位是777 数位和是777的倍数 原数是777的倍数 那么就和777有关。给定QQQ和QQQ个区间[l,r][l,r][l,r]，请求出l,rl,rl,r里面和777无关的数的平方和。 数据 输入 第一行是一个Q(Q&lt;=50)Q(Q&lt;=50)Q(Q<=50) 接下来Q...
杭电多校第三场X Number——数位dp+指数型母函数
mumei的博客
09-19 227
题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=6796 题目大意：给你一个区间[l,r]，让你求出这个区间内满足数位中众数为d的数的个数。 题解：一道披着数位dp皮的组合计数题。 如果用常规数位dp的思路去维护，不难发现状态很难维护，因此我们要考虑换一种思路。 考虑数位dp枚举当前位时，当前位不是数位的上限的情况。 eg:假设数位上限为123456789，当前枚举到第七位，之前的数位为123455（即123455XXX）， 即前面数位有1个1.
hihoCoder 1035 自驾旅行 树形DP
ACMmaxx的专栏
07-21 3293
题目链接：http://hihocoder.com/problemset/problem/1035
codeforces 449D Jzzhu and Numbers 容斥+DP
ACMmaxx的专栏
07-24 1914
很有意思的一题，题目大意
BZOJ 2757: [SCOI2012]Blinker的仰慕者 数位DP
ACMmaxx的专栏
07-11 1856
好繁琐，因为难处理，分为k==0和k！=0两类处理。 代码很丑。。。手残shang
UVALive 6669 Hidden Tree DP
ACMmaxx的专栏
07-17 1609
题目大意：给一个序列a[]，
FZU 2133 kk排座位 轮廓线状态压缩
ACMmaxx的专栏
07-21 1125
题目链接：http://acm.fzu.edu.cn/problem.php?pid=2133
数位类统计问题与数位DP解析
"这篇资源主要介绍了数位动态规划（数位DP）在解决数位类统计问题中的应用，包括基本概念、思想与方法，并通过具体的HDU2089题目进行了实例解析。" 数位动态规划（数位DP）是一种在计算机科学，特别是算法竞赛中...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值