Description

Tweetuzki 特别喜欢等差数列。尤其是公差为 \(1\) 且全为正整数的等差数列。

显然,对于每一个数 \(s\),都能找到一个对应的公差为 \(1\) 且全为正整数的等差数列各项之和为 \(s\)。这时,Tweetuzki 想知道,满足这样条件的等差数列,最小的首项是多少。

由于 Tweetuzki 的数学非常差,尤其是因式分解,所以请你告诉他结果。

Input

输入仅包含一行一个整数 \(s\)

Output

输出一行两个用空格隔开的整数代表首项和末项

Hint

对于 \(10\%\) 的数据,\(1~\leq~s~\leq~10^6\)。

对于 \(100\%\) 的数据,\(1~\leq~s~\leq~10^{12}\)

Solution

考虑前 \(10\%\) 的点,暴力枚举首项,枚举完首项就可以 \(O(1)\) 判断是否合法了。期望得分 10 pts

考虑剩下的部分:

一个等差数列的长度只有为奇数和偶数两种可能,下面对这两种可能分类讨论:

对于长度为奇数的情况,设这个等差数列共有 \(2x~+~1\) 项,其中中项(最中间)为 \(a_k\)

因为 \(a_{k-i}~+~a_{k+i}~=~2~\times~a_k\),所以 \(\sum~a_i~=~(2x~+~1)~a_k~=~s\)

同理,对于长度为偶数的情况,设共有 \(2x\) 项,其中中间两项为 \(a_k~,~a_{k+1}\),因为公差为\(1\),所以即为 \(a_k~,~a_k~+~1\)

于是 \(\sum~a_i~=~x(a_k+a_k~+~1)~=~x~(2a_k~+~1)~=~s\)

由此我们得到了两种情况的关系式

\[s~=~\begin{cases}(2x~+~1)~a_k \\x~(2a_k~+~1) \end{cases}
\]

其中第一种情况长度为奇数,第二种情况长度为偶数。

注意到这两个式子都是 \(s\) 的因数分解,于是考虑直接枚举 \(x\) 的因数。

在第一种情况下,因为 \(a_k\) 是它的因数,我们考虑枚举 \(a_k\),只要枚举到第一个合法的 \(a_k\) 即可停止,

第二种情况下,因为 \(x\) 是它的因数,我们考虑枚举 \(x\),计算所有合法的 \(a_k\)

另外记得特判一个数是由自己做等差数列的情况

Code

#include <cmath>
#include <cstdio>
#include <algorithm>
#ifdef ONLINE_JUDGE
#define freopen(a, b, c)
#define putchar(o) \
puts("I am a cheater!")
#endif
#define rg register
#define ci const int
#define cl const long long typedef long long int ll; namespace IPT {
const int L = 1000000;
char buf[L], *front=buf, *end=buf;
char GetChar() {
if (front == end) {
end = buf + fread(front = buf, 1, L, stdin);
if (front == end) return -1;
}
return *(front++);
}
} template <typename T>
inline void qr(T &x) {
rg char ch = IPT::GetChar(), lst = ' ';
while ((ch > '9') || (ch < '0')) lst = ch, ch=IPT::GetChar();
while ((ch >= '0') && (ch <= '9')) x = (x << 1) + (x << 3) + (ch ^ 48), ch = IPT::GetChar();
if (lst == '-') x = -x;
} template <typename T>
inline void ReadDb(T &x) {
rg char ch = IPT::GetChar(), lst = ' ';
while ((ch > '9') || (ch < '0')) lst = ch, ch = IPT::GetChar();
while ((ch >= '0') && (ch <= '9')) x = x * 10 + (ch ^ 48), ch = IPT::GetChar();
if (ch == '.') {
ch = IPT::GetChar();
double base = 1;
while ((ch >= '0') && (ch <= '9')) x += (ch ^ 48) * ((base *= 0.1)), ch = IPT::GetChar();
}
if (lst == '-') x = -x;
} namespace OPT {
char buf[120];
} template <typename T>
inline void qw(T x, const char aft, const bool pt) {
if (x < 0) {x = -x, putchar('-');}
rg int top=0;
do {OPT::buf[++top] = x % 10 + '0';} while (x /= 10);
while (top) putchar(OPT::buf[top--]);
if (pt) putchar(aft);
} ll s, ans = 1000000000000ll, ss; int main() {
freopen("1.in", "r", stdin);
qr(s);
for (rg ll i = 2; (i * i) <= s; ++i) if(!(s % i)) { //枚举中项
ll k = s / i;
if (k & 1) {
ll x = k >> 1;
if ((i - x) > 0) {
ans = i - x;
ss = k;
break;
}
}
}
for (rg ll i = 1; (i * i) <= s; ++i) if (!(s % i)) { //枚举x
ll k = s / i;
if (k & 1) {
ll a = k >> 1;
if ((a - i) > 0) {
if (ans > a - i) ans = a - i + 1, ss = i << 1;
}
}
}
if (ans == 1000000000000ll) printf("%lld %lld\n", s, s);
else {
qw(ans, ' ', true);
qw(ans + ss - 1, '\n', true);
}
return 0;
}

【数学】【P5077】 Tweetuzki 爱等差数列的更多相关文章

  1. [luogu5077][Tweetuzki 爱等差数列]

    题目链接 思路 数学题 首先列出等差数列求和的式子. \[S = \frac{(n + m)(n - m + 1)}{2}(n为末项,m为首项)\] \[S * 2= (n + m)(n - m + ...

  2. 【数学】【P5076】 Tweetuzki 爱整除

    Description 对于一个数 \(k\),找到任意一个 \(x\),满足 \(0~\leq~k~\leq~x~\leq~10^{18}\) 且对于任意一个 \(x\) 进制数,把该数字各数位上的 ...

  3. 【DP】【P5080】 Tweetuzki 爱序列

    Description Tweetuzki 有一个长度为 \(n\) 的序列 \(a_1~,~a_2~,~\dots~,a_n\). 他希望找出一个最大的 \(k\),满足在原序列中存在一些数 \(b ...

  4. 【贪心】【P5078】Tweetuzki 爱军训

    Description Tweetuzki 所在的班级有 \(n\) 名学生,座号从 \(1\) 到 \(n\).有一次,教官命令班上的 \(n\) 名学生按照座号顺序从左到右排成一排站好军姿,其中 ...

  5. 洛谷 P5078 Tweetuzki 爱军训

    题目连接 很明显,1e6的范围,要么nlgn要么O(n) nlgn的话可能会想到借助一些数据结构,我并没有想到这种做法 对于这种题,O(n)的做法要么是线性递推,要么就应该是贪心了 考虑这道题我们怎么 ...

  6. 【hash】【P5079】P5079 Tweetuzki 爱伊图

    Description Input 第一行两个正整数 \(r~,~c\),表示矩阵的行数和列数. 接下来 \(r\) 行,每行输入 \(c\) 个字符,用空格隔开,保证只含有 . 和 # 两种字符.输 ...

  7. 「Luogu P5080 Tweetuzki 爱序列」

    题目大意 给出一些数,需要求出 \(\frac{a_{i+1}}{3}=a_i\) 或 \(a_{i+1}=2 \times a_i\) 时最长的序列 \(a\). 分析 可以发现符合条件的序列 \( ...

  8. 【剑指offer】和为S的连续整数序列

    找到所有和为S的连续整数序列,序列长度>=2 我的思路:数学法,限定首元素范围,计算序列长度. 书上解法:用small和big两个游标记录序列的开始和结束位置,调整游标. 我的解法: /* ​直 ...

  9. Python 正则表达式 flags 参数

    flags参数 re.I IGNORECASE 忽略字母大小写 re.L LOCALE 影响 “w, “W, “b, 和 “B,这取决于当前的本地化设置. re.M MULTILINE 使用本标志后, ...

随机推荐

  1. Python+selenium+pil+tesseract实现自动识别验证码

    一.环境搭建准备: 1.Python下载,安装以及环境配置 2.IDE pycharm 工具下载,安装 3.ie浏览器 4.selenium 5.pil:pil第三方库的下载,win下安装whl文件, ...

  2. 微信小程序和微信H5测试中易出Bug的点和注意事项

    一.微信小程序 易出Bug的点: 小程序的分享转发功能 背景:小程序项目开发基本完毕也都已经测过几轮,功能上基本没有什么问题,但是上线后却被客户发现通过分享转发小程序给别人,别人无法正常打开的情况 原 ...

  3. webpack整体配置结构

    摘自<深入浅出webpack>2.8 const path = require('path'); module.exports = { // entry 表示入口,webpack执行的第一 ...

  4. 20172329 2018-2019《Java软件结构与数据结构》第一周学习总结

    2018-2019-20172329 <Java软件结构与数据结构>第一周学习总结 在这学期就已经大二了,也已经步入了学习专业课的核心时间,在这个阶段,我们应该了解自己的学习情况,针对自己 ...

  5. uc浏览器的用户体验

    用户界面: 我认为,uc浏览器的用户界面还是很招人喜欢的,可以很容易让用户找到自己想看的网页.简单快捷. 记住用户的选择: uc在每次用户访问完网站之后都会记住用户访问的高频网站,以便下次用户可以更好 ...

  6. App接口如何保证安全

    微信开发或者高德地图,百度地图什么的api要使用,使用之前都需要注册一个账号,然后系统会给你一个key,然后调用api的时候把key传给服务器. 平常公司内部开发项目时,直接用mvc为app客户端提供 ...

  7. lintcode-206-区间求和 I

    206-区间求和 I 给定一个整数数组(下标由 0 到 n-1,其中 n 表示数组的规模),以及一个查询列表.每一个查询列表有两个整数 [start, end] . 对于每个查询,计算出数组中从下标 ...

  8. 初识 es6之 const

    const声明一个只读的常量.一旦声明,常量的值就不能改变. 例子: const a=12; a=2;//报错,const 声明的是常量,不能改 const声明的变量不得改变值,这意味着,const一 ...

  9. 【BioCode】删除未算出PSSM与SS的蛋白质序列

    代码说明: 由于一些原因(氨基酸序列过长),没有算出PSSM与SS,按照整理出来的未算出特征的文件,删除原来的蛋白质序列: 需删除的氨基酸文件732.txt(共732条氨基酸): 删除前 氨基酸共25 ...

  10. $(document).click() 在苹果手机上不能正常运行解决方案

    本来是如下一段跳转代码,发现在安卓和微信开发者工具上都能正常运行,但是苹果手机就不行了. $(document).on('click', '.url', function(){ location.hr ...