题意:
输入样例
本题没有输入
输出样例
2000行,每行一个正整数
输出的每个整数都必须在1~4*10^8之间
如果有多组答案,输出任意一组即可。
思路:
积性函数:对于任意互质的整数a和b有性质f(ab)=f(a)f(b)的数论函数。
完全积性函数:对于任意整数a和b有性质f(ab)=f(a)f(b)的数论函数。
有了上述概念后,那些是积性函数呢?
φ(n) -欧拉函数
μ(n) -莫比乌斯函数,关于非平方数的质因子数目
gcd(n,k) -最大公因子,当k固定的情况
d(n) -n的正因子数目
σ(n) -n的所有正因子之和
因此这个题是求
σ(x*y)大于10,所以转化成
σ(x)*
σ(y)>10就行了。
显然一个质数的正因子个数为2,那么两个质数的积就是4,表示正因子个数为4,所以枚举出前4000个质数,两两组合成2000个数,这样2000个数任意两两相乘就是16>10了。
AC代码:
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <string>
#include <map>
#include <unordered_map>
#include <set>
#include <queue>
#define MAX_INF 0x3f3f3f
#define ll long long
#define ull unsigned long long
const int MAXN = 2e5+10;
const int MOD = 1e9+7;
using namespace std;
int primes[MAXN],cnt =0;
bool is_prime[MAXN];
void init(){
is_prime[1] = true;
for(int i=2;i<=MAXN-5;i++){
if(!is_prime[i])primes[cnt++]=i;
for(int j=0;j<cnt&&i*primes[j]<=MAXN-5;j++){
is_prime[i*primes[j]] = true;
if(i%primes[j] == 0){
break;
}
}
}
}
int main(){
init();
ll ans = 1;
int j = 3999;
for(int i=0;i<2000;i++){
ans = 1LL * primes[j--]*primes[i];
printf("%lld\n",ans);
}
return 0;
}