牛客网暑期ACM多校训练营(第三场)H Diff-prime Pairs [素数筛+前缀和]

版权声明: https://blog.csdn.net/weixin_39792252/article/details/82083352

                                            H Diff-prime Pairs

题目:找满足\frac{i}{gcd(i,j)} ,\frac{j}{gcd(i,j)}都是prime的(i, j)对。

题解:如果gcd(i_{1}, j_{1}) \neq gcd(i_{2}, j_{2}),那么(i_{1}, j_{1}) \neq (i_{2}, j_{2}),这样就枚举g\in [1,n],如果两个素数p_{1}, p_{2}不相等,那么\frac{g*p_{1}}{gcd(g*p_{1}, g*p_{2})}\frac{g*p_{2}}{gcd(g*p_{1}, g*p_{2})}一定满足都是素数,这样就是O(n)的复杂度啦。

代码:

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn = 1e7+7;

int n;
bool vis[maxn];
ll ans[maxn];


void init()
{
    memset(vis, 0, sizeof(vis));
    memset(ans, 0, sizeof(ans));

    for(int i = 2; i*i <= maxn; i++)
        if(!vis[i])
            for(int j = i*i; j <= maxn; j+= i)
                        vis[j] = true;
    ll sum = 0;
    for(int i = 2; i < maxn; i++) if(!vis[i]) ans[i] = 1;
    for(int i = 2; i < maxn; i++) ans[i] += ans[i-1];
}


int main()
{
    init();
    while(~scanf("%d", &n))
    {
        ll sum = 0;
        for(ll i = 1; i <= n; i++) sum += ans[n/i] * (ans[n/i] - 1);
        printf("%lld\n", sum);
    }
    return 0;
}

猜你喜欢

转载自blog.csdn.net/weixin_39792252/article/details/82083352
今日推荐