【2019寒假集训营2】处女座的测验（一）

题意：

输入样例

本题没有输入

输出样例

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;
}