【数论】 HDOJ 5407 CRB and Candies

最新推荐文章于 2017-08-23 14:27:58 发布

yysys

最新推荐文章于 2017-08-23 14:27:58 发布

阅读量557

点赞数

CC 4.0 BY-SA版权

分类专栏：数论文章标签：数论

本文链接：https://blog.youkuaiyun.com/blankcqk/article/details/47830339

数论专栏收录该内容

21 篇文章

订阅专栏

本文探讨了数论中的一个结论，通过编程实现计算1到n+1的最小公倍数除以n+1的过程。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

数论结论题。。。就是1-n+1的lcm除以n+1。。。。不知道结论怎么做？----没法做。。。

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;

const int maxn = 1000005;
const int N = 1000000;
const int mod = 1e9+7;

int val[maxn];
LL f[maxn];
LL g[maxn];
int p[maxn];

LL powmod(LL a, LL b)
{
	LL res = 1, base = a;
	while(b) {
		if(b % 2) res = res * base % mod;
		base = base * base % mod;
		b /= 2;
	}
	return res;
}

void init()
{
	for(int i = 2; i < maxn; i++) if(p[i] == 0)
		for(int j = i + i; j < maxn; j += i) p[j] = 1;
	
	for(int i = 2; i < maxn; i++) if(p[i] == 0)
		for(LL j = (LL)i; j < maxn; j *= i)
			val[j] = i;
	f[1] = 1;
	for(int i = 2; i < maxn; i++) 
		if(val[i]) f[i] = f[i-1] * val[i] % mod;
		else f[i] = f[i-1];
	for(int i = 1; i <= N; i++) g[i] = f[i+1] * powmod(i+1, mod-2) % mod;
}

int main()
{
	init();
	int _;
	scanf("%d", &_);
	while(_--) {
		int x;
		scanf("%d", &x);
		printf("%lld\n", g[x]);
	}
	
	return 0;
}