社区讨论

完美数,玄关

灌水区参与者 2已保存回复 2

讨论操作

快速查看讨论及其快照的属性,并进行相关操作。

当前回复
2 条
当前快照
1 份
快照标识符
@lt9knni5
此快照首次捕获于
2024/03/02 12:17
2 年前
此快照最后确认于
2024/03/02 14:47
2 年前
查看原帖
rt,本人写找完美数的代码,但是long long会爆,球dalao加个高精QWQ
蒟蒻的思路:
2n1为质数,则(n1)n2为完美数若2^n-1为质数,则\frac {(n - 1)n} 2为完美数
CPP
#include <cstdio>
using namespace std;

bool isPrime(long long n)
{
	if (n <= 1)
	{
		return false;
	}
	if (n <= 3)
	{
		return true;
	}
	if (n % 2 == 0 || n % 3 == 0)
	{
		return false;
	}
	for (long long i = 5; i * i <= n; i += 6)
	{
		if (n % i == 0 || n % (i + 2) == 0)
		{
			return false;
		}
	}
	return true;
}

int main()
{
	long long s, x;
	for (int i = 2; ; i++)
	{
		x = (1LL << i) - 1;
		if (isPrime(x))
		{
			s = (x + 1) * x / 2;
			printf("%lld\n", s);
		}
	}
	return 0;
}

回复

2 条回复,欢迎继续交流。

正在加载回复...