P1147 连续自然数和
题目描述
对一个给定的自然数M,求出所有的连续的自然数段,这些连续的自然数段中的全部数之和为M。
例子:1998+1999+2000+2001+2002 = 10000,所以从1998到2002的一个自然数段为M=10000的一个解。
输入输出格式
输入格式:包含一个整数的单独一行给出M的值(10 <= M <= 2,000,000)。
输出格式:每行两个自然数,给出一个满足条件的连续自然数段中的第一个数和最后一个数,两数之间用一个空格隔开,所有输出行的第一个按从小到大的升序排列,对于给定的输入数据,保证至少有一个解。
输入输出样例
输入样例#1:combo.in输出样例#1:
10000
combo.out
18 142
297 328
388 412
1998 2002
小水题,等差数列公式推一推就好了,注意浮点数误差。
1 #include <bits/stdc++.h>
2 const int INF = 0x3f3f3f3f;
3 inline int min(int a,int b){return a > b? b : a;}
4 inline int max(int a,int b){return a < b? b : a;}
5 inline int read(long long &x){
6 x = 0;char ch = getchar();char c = ch;
7 while(ch > '9' || ch < '0')c = ch, ch = getchar();
8 while(ch <= '9' && ch >= '0')x = x * 10 + ch - '0',ch = getchar();
9 if(c == '-') x = -x;
10 }
11 long long m;
12 int main(){
13 read(m);
14 for(long long x = 1;x < m;x ++)
15 {
16 long long k = x * x - x + 2 * m;
17 long long y = (long long)(sqrt(1/4.0 - x + x * x + 2 * m) - 1/2.0);
18 if(y < x)continue;//防止取重
19 //防止误差,分别考察y,y-1,y+1
20 if(y * y + y == k)printf("%lld %lld\n", x, y);
21 y ++;
22 if(y * y + y == k)printf("%lld %lld\n", x, y);
23 y -= 2;
24 if(y * y + y == k)printf("%lld %lld\n", x, y);
25 }
26 return 0;
27 }