hdu 4507 数位dp(求和,求平方和)[通俗易懂]

hdu 4507 数位dp(求和,求平方和)

大家好,又见面了,我是全栈君。

http://acm.hdu.edu.cn/showproblem.php?pid=4507

Problem Description
  单身!

  依旧单身!

  吉哥依旧单身!

  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
   
   
3 1 9 10 11 17 17

 


Sample Output
   
   
236 221 0

 

/***
hdu 4507  数位dp(求和,求平方和)
解题思路:dp[len][sum1][sum2] 表示长度为len对7取模为sum1。各位上的数字和为sum2有多少个满足的数
          一个是与7无关的数的个数。就是简单的数位DP了,非经常规。

第二个与7无关的数的和的维护须要用到第一个个数。 处理到第pos个数位时,加上i*10^pos * 后面的个数 第三个的维护须要用到前面两个 (pre*10^pos + next)^2= (pre*10^pos)^2+2*pre*10^pos*next +next^2*/#include <stdio.h>#include <algorithm>#include <string.h>#include <iostream>using namespace std;typedef long long LL;const LL mod=1e9+7;LL l,r,p[25];int bit[25];struct node{ LL cnt,sum,sqsum;} dp[25][10][10];node dfs(int len,int sum1,int sum2,int flag){ if(len<0) { node tmp; tmp.cnt=(sum1!=0&&sum2!=0); tmp.sum=tmp.sqsum=0; return tmp; } if(flag==0&&dp[len][sum1][sum2].cnt!=-1)return dp[len][sum1][sum2]; node ans,tmp; int end=flag?bit[len]:9; ans.cnt=ans.sqsum=ans.sum=0; for(int i=0; i<=end; i++) { if(i==7)continue; tmp=dfs(len-1,(sum1+i)%7,(sum2*10+i)%7,flag&&i==end); ans.cnt+=tmp.cnt; ans.cnt%=mod; ans.sum+=(tmp.sum+i*p[len]%mod*tmp.cnt%mod)%mod; ans.sum%=mod; ans.sqsum+=(tmp.sqsum+2*p[len]*i%mod*tmp.sum%mod)%mod; ans.sqsum%=mod; ans.sqsum+=(tmp.cnt*p[len]%mod*p[len]%mod*i*i%mod); ans.sqsum%=mod; } if(flag==0)dp[len][sum1][sum2]=ans; return ans;}LL solve(LL n){ int len=0; while(n) { bit[len++]=n%10; n/=10; } return dfs(len-1,0,0,1).sqsum;}int main(){ p[0]=1; for(int i=1; i<20; i++) p[i]=(p[i-1]*10)%mod; for(int i=0; i<25; i++) { for(int j=0; j<10; j++) { for(int k=0; k<10; k++) { dp[i][j][k].cnt=-1; } } } int T; scanf("%d",&T); while(T--) { scanf("%I64d%I64d",&l,&r); printf("%I64d\n",((solve(r)-solve(l-1))%mod+mod)%mod); } return 0;}

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/116251.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)


相关推荐

  • 使用 SQL Server 的 uniqueidentifier 字段类型[通俗易懂]

    使用 SQL Server 的 uniqueidentifier 字段类型[通俗易懂]SQLServer自2008版起引入了 uniqueidentifier 字段,它存储的是一个UUID,或者叫GUID,内部存储为16个字节。SQLServer可用两个函数来生成 uniqueidentifier,分别是 NEWID() 和 NEWSEQUENTIALID(),后者只能用作字段的默认值。Java也有一个UUID工具类 java.uti.UUID, U…

  • idea2021.2.8激活码-激活码分享

    (idea2021.2.8激活码)这是一篇idea技术相关文章,由全栈君为大家提供,主要知识点是关于2021JetBrains全家桶永久激活码的内容IntelliJ2021最新激活注册码,破解教程可免费永久激活,亲测有效,下面是详细链接哦~https://javaforall.cn/100143.html83PVI25FMO-eyJsa…

  • genre-based_deepsort特征判断

    genre-based_deepsort特征判断  https://arxiv.org/abs/1804.01438https://blog.csdn.net/gavinmiaoc/article/details/80648754https://zhuanlan.zhihu.com/p/35296881 https://github.com/seathiefwang/MGN-pytorch

  • 【金融市场基础知识】——中国的金融体系(一)[通俗易懂]

    【金融市场基础知识】——中国的金融体系(一)[通俗易懂]阅读之前看这里????:博主是一名正在学习证券知识的学生,在每个领域我们都应当是学生的心态,也不应该拥有身份标签来限制自己学习的范围,所以博客记录的是在学习过程中一些总结,也希望和大家一起进步,在记录之时,未免存在很多疏漏和不全,如有问题,还请私聊博主指正。博客地址:天阑之蓝的博客,学习过程中不免有困难和迷茫,希望大家都能在这学习的过程中肯定自己,超越自己,最终创造自己。目录中国的金融体系(一)一、中国金融市场的历史、现状及影响因素1、新中国成立以来我国金融市场的发展历史★2、我国金融市场的发展现状

  • 计算机前沿技术介绍

    计算机前沿技术介绍简要的说明了一些当今,热门、前沿技术,只是大致介绍了一下,想详细了解的,请自行查阅。

    2022年10月19日
  • R语言plot参数_plot函数参数

    R语言plot参数_plot函数参数最近用R语言画图,plot函数是用的最多的函数,而他的参数非常繁多,由此总结一下,以供后续方便查阅。plot(x,y=NULL,type="p",xlim=NULL,ylim=NULL,log="",main=NULL,sub=NULL,xlab=NULL,ylab=NULL,ann=par("ann"),axes=TRUE,fra…

    2022年10月15日

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号