Codeforces-----233B---Non-square Equation暴力枚举

Non-square Equation

time limit per test
1 second
memory limit per test
256 megabytes
input
standard input
output
standard output

Let's consider equation:

x2 + s(xx - n = 0, 

where x, n are positive integers, s(x) is the function, equal to the sum of digits of number x in the decimal number system.

You are given an integer n, find the smallest positive integer root of equation x, or else determine that there are no such roots.

Input

A single line contains integer n (1 ≤ n ≤ 1018) — the equation parameter.

Please, do not use the %lld specifier to read or write 64-bit integers in С++. It is preferred to use cin, cout streams or the %I64d specifier.

Output

Print -1, if the equation doesn't have integer positive roots. Otherwise print such smallest integer x (x > 0), that the equation given in the statement holds.

Examples
Input
2
Output
1
Input
110
Output
10
Input
4
Output
-1
Note

In the first test case x = 1 is the minimum root. As s(1) = 1 and 12 + 1·1 - 2 = 0.

In the second test case x = 10 is the minimum root. As s(10) = 1 + 0 = 1 and 102 + 1·10 - 110 = 0.

In the third test case the equation has no roots.

s(x)是x的各个位数字之和,因为x*x+x*s(x)-n = 0,即x*x+x*s(x)=n,n <= 1e18,所以x <= 999999999,即s(x) <= 81

且x,n为正整数,所以遍历s(x)即可,再根据一元二次方程求根公式(-b+sqrt(b*b-4*a*c))/ 2 判断s(x)与x是否满足即可

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstring>
#include<string>
#include<algorithm>
#include<map>
#include<queue>
#include<vector>
#define PI acos(-1.0)
#define INF 0x3f3f3f3f
#define CL(a, b) memset(a, b, sizeof(a))
using namespace std;
typedef __int64 LL;
const int MAXN = 5e5+10;
const int MOD = 1e9+7; 
char s[50], str[50];
bool judge(LL a, LL b, LL c){
	return a*a + a*b == c;
}
bool Judge(LL a, LL b){
	LL ans = 0;
	while(a){
		ans += a%10;
		a /= 10;
	}
	return ans == b;
}
int main(){  
    LL t, x, c, ans, kcase = 1;
    scanf("%I64d", &t);
	bool flag = false;
	for(LL i = 1; i <= 81; i++){
		LL v = i*i+4*t;
		if(sqrt(v)*sqrt(v) == v){//x为正整数,v必定是完全平方数
			x = (sqrt(v)-i)/2;
			if(judge(x, i, t) && Judge(x, i)){
				flag = true;
				break;
			}
		}
	}
	printf("%I64d\n", flag ? x : -1);
    return 0;  
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值