2017 ACM/ICPC Asia Regional Qingdao Online 1011 A Cubic number and A Cubic Number
Posted 九月旧约
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了2017 ACM/ICPC Asia Regional Qingdao Online 1011 A Cubic number and A Cubic Number相关的知识,希望对你有一定的参考价值。
A Cubic number and A Cubic Number
Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 65535/32768 K (Java/Others)
Total Submission(s): 0 Accepted Submission(s): 0
For each test case, a line contains a prime number p (2≤p≤1012).
给出一个p,p是素数(2 <= p <=1e12),判断是否存在正整数a,b使得
a*a*a - b*b*b = p
如果存在输出YES,否则输出NO。
我们知道立方差公式:
a*a*a - b*b*b = (a-b)(a*a+a*b+b*b)
我们知道a*a*a - b*b*b = prime(素数)
所以(a-b)(a*a+a*b+b*b) = prime;
由因为prime 只能由1*prime得来。
所以prime = 1*prime
所以(a-b) = 1, (a*a+a*b+b*b) = prime , a,b必是两个相邻的整数。
设b = x,a = x+1
则 (x+1)^3-x^3 = 3*x*x + 3*x + 1 = prime
3*x*x + 3*x = prime-1
x*x + x = (prime-1)/3
x(x+1) = (prime-1)/3
令:(prime-1)/3 = T;
x = int(sqrt(T)) //解出一个x
若x*(x+1) == T
则说明存在一个整数a,是的a^3 - (a-1)^3 等于当前所输入的素数。
这个a就是我们解出的x再加1.
#include <iostream> #include <stdio.h> #include <math.h> using namespace std; int main() { int T; scanf("%d",&T); while(T--) { double n,prime; scanf("%lf",&prime); prime = prime - 1; n = prime/3; int a=(int)(sqrt(n)); int b=a+1; if((double)a*b==n) printf("YES\n"); else printf("NO\n"); } return 0; }
以上是关于2017 ACM/ICPC Asia Regional Qingdao Online 1011 A Cubic number and A Cubic Number的主要内容,如果未能解决你的问题,请参考以下文章
2017 ACM/ICPC Asia Regional Beijing Online
2017 ACM/ICPC Asia Regional Qingdao Online
2017 ACM/ICPC Asia Regional Shenyang Online spfa+最长路
hdu6206 Apple 2017 ACM/ICPC Asia Regional Qingdao Online
HDU 6198(2017 ACM/ICPC Asia Regional Shenyang Online)
HDU - 6215 2017 ACM/ICPC Asia Regional Qingdao Online J - Brute Force Sorting