P1147 连续自然数和(枚举,数学/二分)

2024-03-19 17:32

本文主要是介绍P1147 连续自然数和(枚举,数学/二分),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

https://www.luogu.org/problemnew/show/P1147
题目描述
对一个给定的自然数M,求出所有的连续的自然数段,这些连续的自然数段中的全部数之和为M。

例子:1998+1999+2000+2001+2002 = 10000,所以从1998到2002的一个自然数段为M=10000的一个解。

输入输出格式
输入格式:
包含一个整数的单独一行给出M的值(10≤M≤2,000,000)。

输出格式:
每行两个自然数,给出一个满足条件的连续自然数段中的第一个数和最后一个数,两数之间用一个空格隔开,所有输出行的第一个按从小到大的升序排列,对于给定的输入数据,保证至少有一个解。

输入输出样例
输入样例#1:

10000

输出样例#1:

18 142 
297 328 
388 412 
1998 2002

/*
way1:数据不大,可以直接考虑枚举+二分
way2:更快的肯定就是数学推导了
*/
way1(O(nlogn)):

#include <iostream>
using namespace std;
typedef long long LL;
int main()
{LL m;cin>>m;for(int i = 0; i < m; i++){LL s = i,e  = -1;LL low = s+1,high = m-1,mid;while(low <= high){mid = (low + high) / 2;LL temp = (s + mid) * (mid - s + 1);//这里int不够if(temp < 2 * m){low = mid + 1;}else if(temp > 2 * m){high = mid - 1;}else{e = mid;break;}}if(e != -1){cout<<s<<" "<<e<<endl;}}return 0;
}

way2 (O(n)):
数学推导:
根据等差数列求和:

sum = (s+e) * (e-s+1) / 2
(s+e) *(e-s+1) = 2 * sum
设 s + e = x1,e-s+1 = x2;解上述方程组:e = (x1+x2-1)/2, s = (x1-x2+1)/2;再由第二个表达式:x1 * x2 = 2*sum可知2*sum一定为偶数,那么x1,x2可能有两种情况:1.两个都是偶数2.一奇一偶根据 e = (x1+x2-1)/2, s = (x1-x2+1)/2; 并且要是自然数解,那么x1+x2-1与x1-x2+1必定要为偶数,so可以排除第一种情况,所以x1,x2必定一奇一偶根据x1,x2必定一奇一偶,枚举其中一个,根据x1 * x2 = 2*sum得出另一个,再根据 e = (x1+x2-1)/2, s = (x1-x2+1)/2就可以求出结果了。比如枚举x2,根据x1 * x2 = 2*sum,确定x2的上限:up <= sqrt(2*sum),后续看代码~
#include <iostream>
#include <math.h>
using namespace std;
typedef long long LL;//其实int就够
int main()
{LL m;cin>>m;LL s,e,x1,x2;LL up= sqrt(2*m);for(int i = up; i > 1; i--)//逆序枚举x2,因为影响s,e起主导作用的是x1{x2 = i;x1 = 2 * m / x2;if(((x2&1) && !(x1&1)) || (!(x2&1) && (x1&1))){if(x1 * x2 == 2*m){s = (x1 - x2 + 1) >> 1;e = (x1 + x2 - 1) >> 1;cout<<s<<" "<<e<<endl;}}}return 0;
}

这篇关于P1147 连续自然数和(枚举,数学/二分)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/826750

相关文章

hdu2241(二分+合并数组)

题意:判断是否存在a+b+c = x,a,b,c分别属于集合A,B,C 如果用暴力会超时,所以这里用到了数组合并,将b,c数组合并成d,d数组存的是b,c数组元素的和,然后对d数组进行二分就可以了 代码如下(附注释): #include<iostream>#include<algorithm>#include<cstring>#include<stack>#include<que

hdu2289(简单二分)

虽说是简单二分,但是我还是wa死了  题意:已知圆台的体积,求高度 首先要知道圆台体积怎么求:设上下底的半径分别为r1,r2,高为h,V = PI*(r1*r1+r1*r2+r2*r2)*h/3 然后以h进行二分 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#includ

poj2406(连续重复子串)

题意:判断串s是不是str^n,求str的最大长度。 解题思路:kmp可解,后缀数组的倍增算法超时。next[i]表示在第i位匹配失败后,自动跳转到next[i],所以1到next[n]这个串 等于 n-next[n]+1到n这个串。 代码如下; #include<iostream>#include<algorithm>#include<stdio.h>#include<math.

uva 10014 Simple calculations(数学推导)

直接按照题意来推导最后的结果就行了。 开始的时候只做到了第一个推导,第二次没有继续下去。 代码: #include<stdio.h>int main(){int T, n, i;double a, aa, sum, temp, ans;scanf("%d", &T);while(T--){scanf("%d", &n);scanf("%lf", &first);scanf

uva 10025 The ? 1 ? 2 ? ... ? n = k problem(数学)

题意是    ?  1  ?  2  ?  ...  ?  n = k 式子中给k,? 处可以填 + 也可以填 - ,问最小满足条件的n。 e.g k = 12  - 1 + 2 + 3 + 4 + 5 + 6 - 7 = 12 with n = 7。 先给证明,令 S(n) = 1 + 2 + 3 + 4 + 5 + .... + n 暴搜n,搜出当 S(n) >=

uva 11044 Searching for Nessy(小学数学)

题意是给出一个n*m的格子,求出里面有多少个不重合的九宫格。 (rows / 3) * (columns / 3) K.o 代码: #include <stdio.h>int main(){int ncase;scanf("%d", &ncase);while (ncase--){int rows, columns;scanf("%d%d", &rows, &col

hdu 2489 (dfs枚举 + prim)

题意: 对于一棵顶点和边都有权值的树,使用下面的等式来计算Ratio 给定一个n 个顶点的完全图及它所有顶点和边的权值,找到一个该图含有m 个顶点的子图,并且让这个子图的Ratio 值在所有m 个顶点的树中最小。 解析: 因为数据量不大,先用dfs枚举搭配出m个子节点,算出点和,然后套个prim算出边和,每次比较大小即可。 dfs没有写好,A的老泪纵横。 错在把index在d

【生成模型系列(初级)】嵌入(Embedding)方程——自然语言处理的数学灵魂【通俗理解】

【通俗理解】嵌入(Embedding)方程——自然语言处理的数学灵魂 关键词提炼 #嵌入方程 #自然语言处理 #词向量 #机器学习 #神经网络 #向量空间模型 #Siri #Google翻译 #AlexNet 第一节:嵌入方程的类比与核心概念【尽可能通俗】 嵌入方程可以被看作是自然语言处理中的“翻译机”,它将文本中的单词或短语转换成计算机能够理解的数学形式,即向量。 正如翻译机将一种语言

XTU 1233 n个硬币连续m个正面个数(dp)

题面: Coins Problem Description: Duoxida buys a bottle of MaiDong from a vending machine and the machine give her n coins back. She places them in a line randomly showing head face or tail face o

poj 2976 分数规划二分贪心(部分对总体的贡献度) poj 3111

poj 2976: 题意: 在n场考试中,每场考试共有b题,答对的题目有a题。 允许去掉k场考试,求能达到的最高正确率是多少。 解析: 假设已知准确率为x,则每场考试对于准确率的贡献值为: a - b * x,将贡献值大的排序排在前面舍弃掉后k个。 然后二分x就行了。 代码: #include <iostream>#include <cstdio>#incl