1590 恨 7 不成妻(HDU4507 LOJ10168) 暴力0分 能处理平方和数位DP(包含数本身 数位和 两种状态) 到处要模否则溢出

最新推荐文章于 2024-09-22 11:42:19 发布
最新推荐文章于 2024-09-22 11:42:19 发布
阅读量432 收藏
点赞数
分类专栏： 信息学奥赛一本通 提高篇 文章标签： 动态规划 算法 c++
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/mrcrack/article/details/123476097
版权

总目录

在线测评地址(ybt)

在线测评地址(LOJ)

1.暴力0分

ybt

未通过

测试点结果内存时间
测试点1运行超时576KB997MS
测试点2运行超时592KB998MS

LOJ

 暴力0分代码如下：

#include <bits/stdc++.h>
#define LL long long
#define mod 1000000007
using namespace std;
int main(){
	int t,flag,sum;
	LL lt,rt,i,x,ans;
	scanf("%d",&t);
	while(t--){
		scanf("%lld%lld",&lt,&rt);
		ans=0;
		for(i=lt;i<=rt;i++){
			x=i;
			if(x%7==0)continue;
			flag=0,sum=0;
			while(x){
				if(x%10==7){
					flag=1;
					break;
				}else sum+=x%10;
				x/=10;
			}
			if(flag)continue;
			if(sum%7==0)continue;
			ans=(ans+i*i%mod)%mod;
		}
		printf("%lld\n",ans);
	}
	return 0;
}

2.数位DP

(1)整数中某一位是 7；(数位DP很容易编写)

(2)整数的每一位加起来的和是 7的整数倍；(作为中间值进行传递，数位DP也可编写)

(3)这个整数是 7的整数倍。(作为中间值进行传递，数位DP也可编写，不过编写没把握)

困难:将(2)(3)综合起来比较困难，还有dp数组应该如何设计.dp[pos][]另一个维度应该如何设计？

我们先定义结构体，其中存储三个变量。

  1. 个数 cnt

  2. 和 sum

  3. 平方和 sqsum

定义dp[i][j][k] 表示前i位数,各位数之和mod7（j）,这个数mod7中还能合法的数的个数(k)。

然后涮刷刷一边数位dp dfs的模板

之后就到了最难的一点，如何统计平方和

以两位数2X为例：
初始化
对应的个数 cnt=0
对应的和 sum=0
对应的平方和 sqsum=0

第0位符合题意的数据是1,2,3,4,5,6,8,9
对应的个数 cnt=8
对应的和 sum=1+2+3+4+5+6+8+9=38
对应的平方和 sqsum=1*1+2*2+3*3+4*4+5*5+6*6+8*8+9*9=236

到达第1位0的计算，符合题意的数据是01,02,03,04,05,06,08,09
对应的个数 cnt=8
对应的和 sum=1+2+3+4+5+6+8+9=38
对应的平方和 sqsum=1*1+2*2+3*3+4*4+5*5+6*6+8*8+9*9=236

到达第1位1的计算，符合题意的数据是01,02,03,04,05,06,08,09,10,11,12,13,15,18,19
对应的个数 cnt=15
对应的和 sum=1+2+3+4+5+6+8+9+10+11+12+13+15+18+19=136
对应的平方和 sqsum=1*1+2*2+3*3+4*4+5*5+6*6+8*8+9*9+10*10+11*11+12*12+13*13+15*15+18*18+19*19
=1680

到达第1位2的计算，符合题意的数据是01,02,03,04,05,06,08,09,10,11,12,13,15,18,19,20,22,23,24,26,29
对应的个数 cnt=21
对应的和 sum=1+2+3+4+5+6+8+9+10+11+12+13+15+18+19+20+22+23+24+26+29=280
对应的平方和 sqsum=1*1+2*2+3*3+4*4+5*5+6*6+8*8+9*9+10*10+11*11+12*12+13*13+15*15+18*18+19*19+20*20+22*22+23*23+24*24+26*26+29*29=5186

我们的转移是这样滴

首先关于个数的转移是就直接把后面的叠加不就完了，标准转移。

以两位数2X为例：
初始化
对应的个数 cnt=0


第0位符合题意的数据是1,2,3,4,5,6,8,9
对应的个数 cnt=8


到达第1位0的计算，符合题意的数据是01,02,03,04,05,06,08,09
对应的个数 cnt=8


到达第1位1的计算，符合题意的数据是01,02,03,04,05,06,08,09,10,11,12,13,15,18,19
01,02,03,04,05,06,08,09有8个数
10,11,12,13,15,18,19有7个数
对应的个数 cnt=8+7=15


到达第1位2的计算，符合题意的数据是01,02,03,04,05,06,08,09,10,11,12,13,15,18,19,20,22,23,24,26,29
01,02,03,04,05,06,08,09,10,11,12,13,15,18,19有15个数
20,22,23,24,26,29有6个数
对应的个数 cnt=15+6=21

和的转移则是稍微动下脑子就可以轻松写出来，sum_now += cnt_son*(10^pos*i)。

以两位数2X为例：
初始化
对应的个数 cnt=0
对应的和 sum=0


第0位符合题意的数据是1,2,3,4,5,6,8,9
对应的个数 cnt=8
对应的和 sum=1+2+3+4+5+6+8+9=38


到达第1位0的计算，符合题意的数据是01,02,03,04,05,06,08,09
对应的个数 cnt=8
对应的和 sum=1+2+3+4+5+6+8+9=38


到达第1位1的计算，符合题意的数据是01,02,03,04,05,06,08,09,10,11,12,13,15,18,19
对应的个数 cnt=8+7=15
对应的和 
sum=1+2+3+4+5+6+8+9+10+11+12+13+15+18+19
=38+(10+11+12+13+15+18+19)
=38+(10+10+1+10+2+10+3+10+5+10+8+10+9)
=38+(10*7+1+2+3+5+8+9)
=38+(10*7+28)
=136


到达第1位2的计算，符合题意的数据是01,02,03,04,05,06,08,09,10,11,12,13,15,18,19,20,22,23,24,26,29
对应的个数 cnt=15+6=21
对应的和 
sum=1+2+3+4+5+6+8+9+10+11+12+13+15+18+19+20+22+23+24+26+29
=136+(20+22+23+24+26+29)
=136+(20*6+2+3+4+6+9)
=136+(20*6+24)
=280

当前位数对于和的贡献是存在于之后的都合法，那么这一位对于和的叠加就应该是它本身的值乘上之后合法的数目。（就是有多少加多少。。实在看不懂琢磨琢磨代码应该就懂了）在基于这两点后，我们对于sqsum的求解就可以借由一个初中大家都学过的公式了

(a+b)^2 = a^2+2ab+b^2;

显然之前的sqsum已经将之前所有数的平方和都求完了，实际上我们求出的是b1^2+b2^2+b3^2……bn^2。

那么我们发现在这一位的a就是i*10^pos，并且对于每一个b都相同，所以转移可以这么写。

sqsum_now = sqsum_son + 2*sum_son*(10^pos*i) +(10^pos*i)*(10^pos*i)*cnt_son;

以两位数2X为例：
初始化
对应的个数 cnt=0
对应的和 sum=0
对应的平方和 sqsum=0

第0位符合题意的数据是1,2,3,4,5,6,8,9
对应的个数 cnt=8
对应的和 sum=38
对应的平方和 sqsum=1*1+2*2+3*3+4*4+5*5+6*6+8*8+9*9=236

到达第1位0的计算，符合题意的数据是01,02,03,04,05,06,08,09
对应的个数 cnt=8
对应的和 sum=38
对应的平方和 sqsum=1*1+2*2+3*3+4*4+5*5+6*6+8*8+9*9=236

到达第1位1的计算，符合题意的数据是01,02,03,04,05,06,08,09,10,11,12,13,15,18,19
对应的个数 cnt=8+7=15
对应的和 sum==38+(10*7+1+2+3+5+8+9)
=38+(10*7+28)
=136
对应的平方和 sqsum=1*1+2*2+3*3+4*4+5*5+6*6+8*8+9*9+10*10+11*11+12*12+13*13+15*15+18*18+19*19
=236+(10*10+11*11+12*12+13*13+15*15+18*18+19*19)
=236+(10*10+(10+1)*(10+1)+(10+2)*(10+2)+(10+3)*(10+3)+(10+5)*(10+5)+(10+8)*(10+8)+(10+9)*(10+9))
=236+(10*10*7+2*10*(1+2+3+5+8+9)+(1*1+2*2+3*3+5*5+8*8+9*9))
=236+(10*10*7+2*10*28+184)
=1680

到达第1位2的计算，符合题意的数据是01,02,03,04,05,06,08,09,10,11,12,13,15,18,19,20,22,23,24,26,29
对应的个数 cnt=15+6=21
对应的和 sum=136+(20*6+2+3+4+6+9)
=136+(20*6+24)
=280

对应的平方和 sqsum=1*1+2*2+3*3+4*4+5*5+6*6+8*8+9*9+10*10+11*11+12*12+13*13+15*15+18*18+19*19
=1680+(20*20+22*22+23*23+24*24+26*26+29*29)
=1680+(20*20+(20+2)*(20+2)+(20+3)*(20+3)+(20+4)*(20+4)+(20+6)*(20+6)+(20+9)*(20+9))
=1680+(20*20*6+2*20*(2+3+4+6+9)+(2*2+3*3+4*4+6*6+9*9))
=1680+(20*20*6+2*20*24+146)
=5186

至此，所有状态转移完毕.

编码过程中，还有一个难点，到处要模否则溢出。

ybt

通过

测试点结果内存时间
测试点1答案正确660KB5MS
测试点2答案正确652KB7MS

LOJ

 能处理平方和数位DP代码如下：

#include <bits/stdc++.h>
#define LL long long
#define mod 1000000007
using namespace std;
struct node{
	LL cnt,sum,sqsum;//数量,和,平方和 
}dp[20][10][10];//dp[pos][num%7][sum%7]//num数本身,sum数位和 
int bit[20];
LL p[20];
node dfs(int pos,int num,int sum,int limit){
	node ans,tmp;
	int i,up;
	if(pos==-1){
		tmp.cnt=(num!=0&&sum!=0);
		tmp.sum=0;
		tmp.sqsum=0;
		return tmp;
	}
	if(!limit&&dp[pos][num][sum].cnt!=-1)return dp[pos][num][sum];
	up=limit?bit[pos]:9;
	ans.cnt=ans.sum=ans.sqsum=0;
	for(i=0;i<=up;i++){
		if(i==7)continue;
		tmp=dfs(pos-1,(num*10+i)%7,(sum+i)%7,limit&&i==up);
		ans.cnt+=tmp.cnt;
		ans.cnt%=mod;
		ans.sum+=((p[pos]*i)%mod*tmp.cnt%mod+tmp.sum)%mod;
		ans.sum%=mod;
		ans.sqsum+=(p[pos]*i)%mod*(p[pos]*i)%mod*tmp.cnt%mod;
		ans.sqsum+=(2*(p[pos]*i)%mod*tmp.sum%mod+tmp.sqsum)%mod;
		ans.sqsum%=mod;
	}
	if(!limit)dp[pos][num][sum]=ans;
	return ans;
}
LL solve(LL x){
	int pos=0;
	while(x){
		bit[pos++]=x%10;
		x/=10;
	}
	return dfs(pos-1,0,0,1).sqsum;
}
int main(){
	int T,i;
	LL lt,rt;
	scanf("%d",&T);
	p[0]=1;
	for(i=1;i<=18;i++)p[i]=p[i-1]*10%mod;
	while(T--){
		LL ret;
		scanf("%lld%lld",&lt,&rt);
		memset(dp,-1,sizeof(dp));
		ret=((solve(rt)-solve(lt-1))%mod+mod)%mod;//避免出现负数的常见手法 
		printf("%lld\n",ret);
	}
	return 0;
}

mrcrack
关注
专栏目录
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...
1590：恨 7 不成妻(数位DP）
RunningBeef的博客
08-18 444
这题有一个难点一个坑点,and 你要细点 难点 你要用DP维护出区间内所有和7无关的数的平方和。 这个需要你做一个数学转换。 因为你用DP预处理一个区间的所有与7无关的数，就是通过固定最高位，枚举次高位来递推。所以这个时候我们 可以试着把一个数拆成最高位和剩余位的和，看看它有没有什么特点，从而得到DP的状态转移方程。 所以一个数可以拆成 t∗10k+xi t*10^{k} + x_{i} t∗10k+xi​ 则有固定最高位，低位,每个数的平方和如下 (x1+t∗10k)2+(x2+t∗10k)2+⋯⋯(
一本通1590恨 7 不成妻
weixin_30315435的博客
02-11 2536
1590：恨 7 不成妻 时间限制: 1000 ms 内存限制: 524288 KB 【题目描述】 单身！ 依然单身！ 吉哥依然单身！ DS 级码农吉哥依然单身！ 所以，他平生最恨情人节，不管是214还是77，他都讨厌！ 吉哥观察了214和77这两个数，发现： 2+1+4=7 7+7=7×2 77=7×11 最终，他发现原来这一切归根到底都...
【数位DP】[LOJ10168] 恨7不成妻
weixin_30302609的博客
09-25 185
还是数位DP... 状态：$f[x][val][sum]$表示当前第x位，当前数字为val，当前各位数字和为sum 观察到$val$,$sum$过大，很套路地模7即可... 每个状态存储三个要用到的值： 　　①cnt:表示满足条件的数字的个数 　　②s:表示满足条件的数字的和 　　③s_2:表示满足条件的数字的平方的和 关于转移： 　　①cnt的转移就和一般的数位DP一样...
[HDU]4507 恨7不成妻 数位Dp好题
Le Petit Prince
10-18 607
吉哥系列故事——恨7不成妻Time Limit: 1000/500 MS (Java/Others) Memory Limit: 65535/32768 K (Java/Others) Total Submission(s): 4851 Accepted Submission(s): 1590Problem Description 　　单身! 　　依然单身！ 　　吉哥依然单身！
HDU 4507 吉哥系列故事——恨7不成妻（数位DP）
窝不是爱酱,喵~~~~
03-22 5882
转载请注明出处，谢谢http://blog.csdn.net/ACM_cxlove?viewmode=contents    by---cxlove Tencent马拉松初赛第二场的鬼畜题。。。 其实还是比较裸的数位DP。 给出的三种定义，都可以用位DP解决。而且是比较基础的。 至于平方和，大概是这样的。 我们假设dfs进下一层的数字是x，那么当前长度是len的话，最高位添加的数字是
DP学习笔记
killer的博客
08-13 368
文章目录记忆化搜索背包**一.01背包Easy\color{green}{Easy}Easy**二.完全背包Easy\color{green}{Easy}Easy三.多重背包Mid\color{orange}{Mid}Mid四.分组背包Easy\color{green}{Easy}Easy线性DP[数字三角形Easy\color{green}{Easy}Easy](https://www.acwing.com/problem/content/900/)[最长上升子序列LISMid\color{orange}
数位dp——某次队内分享任务
07-21
数位动态规划(Digit Dynamic Programming, 简称数位DP)是一种在处理与数字相关的动态规划问题时，通过分析每个数字位上的状态来优化算法的技术。它通常用于解决涉及数字序列特性的数学问题，例如计数特定类型的数...
初探数位dp.pdf
09-16
lazycal的集训队报告：初探数位DP 以HDU 2089，HDU 3652, URAL 1057等题目为例，介绍了数位DP的算法
hdu 4507 吉哥系列故事——恨7不成妻 题解（数位DP（超难））
LightningUZ的博客
05-28 417
原题链接： 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...
【数位DP】B-number HDU - 3652
BlessingXRY的博客
09-27 420
Think: 1知识点：数位DP 2题意：输入r，判断[1, r]区间内有多少个“wqb-number”，“wqb-number”的定义为数位内含有13且可以被13整除，例如13,2613是“wqb-number”，而143和131不是“wqb-number”，143中1和3不是连续存在，131不能被13整除。 3反思： （1）：dp数组的控制条件应与dfs中控制变量相对应 （2）：dp[
HDU 4507 恨7不成妻 （数位DP套路题，详细解析）
pursuit的博客
05-13 2801
不会数位DPDPDP的这里指路一篇介绍非常详细的数位DPDPDP的blogblogblog:点这里。 链接 恨7不成妻 题面 单身! 依然单身！ 吉哥依然单身！ DS级码农吉哥依然单身！ 所以，他生平最恨情人节，不管是214还是77，他都讨厌！ 吉哥观察了214和77这两个数，发现： 2+1+4=72+1+4=72+1+4=7　 7+7=7∗27+7=7*27+7=7∗2 77=7∗1177=7*1177=7∗11 最终，他发现原来这一切归根到底都是因为和7有关！所以，他现在甚至讨厌一切和7有关的
loj题目总览
djrlb735733的博客
02-28 749
——DavidJing提供技术支持 现将今年7月份之前必须刷完的题目列举 完成度[25/34] [195/250] 第 1 章 贪心算法 √ [11/11] #10000 「一本通 1.1 例 1」活动安排 #10001 「一本通 1.1 例 2」种树 #10002 「一本通 1.1 例 3」喷水装置 #10003 「一本通 1.1 例 4」加工生产调度 #10004 「一本通 1.1 ...
loj#10168恨 7 不成妻
quantum11
08-28 320
这一题的题目描述很…无法描述，但还是能马上看出是一道数位dp题。限制的条件有点多，细节处理起来确实很烦。只要记录合法的数的个数，数的和以及数的平方和。 #include&lt;bits/stdc++.h&gt; #define ll long long #define il inline using namespace std; const int mo=1e9+7; #define getch...
atcoder abc372 启发式合并， dp
最新发布
2201_75845839的博客
09-22 372
现在考虑如何存数据。可以用并查集合并两个集合，并且在合并的时候一定要将小的集合合并到大的集合中去，因为小集合合并到大集合中，最后的集合一定大于小集合元素数量的两倍，因此合并的时间复杂度不会超过logn。思路：判断当前变化的位置对结果是否会产生影响，在询问前先计算有多少ABC, 每次变换位置x前后扫描一下[x - 2, x + 2]范围内是否有ABC，变换前有让cnt --， 变换后有让cnt ++思路：三进制转换，可以参考二进制，先把当前可以加入的最大的3的幂次加入,这样一定可以凑出答案。
代码随想录算法训练营Day7
lynyzy的博客
09-18 541
双指针
力扣 16.最接近的三数之和
qq_51352130的博客
09-16 441
设 s=nums[i]+nums[j]+nums[k]，为了判断 s 是不是与 target 最近的数，我们还需要用一个变量 minDiff 维护 ∣s−target∣ 的最小值。分类讨论：如果 s=target，那么答案就是 s，直接返回 s。如果 s>target，那么如果 s−target<minDiff，说明找到了一个与 target 更近的数，更新 minDiff 为 s−target，更新答案为 s。然后和三数之和一样，把 k 减一。
数位类统计问题与数位DP解析
"这篇资源主要介绍了数位动态规划（数位DP）在解决数位类统计问题中的应用，包括基本概念、思想与方法，并通过具体的HDU2089题目进行了实例解析。" 数位动态规划（数位DP）是一种在计算机科学，特别是算法竞赛中...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值