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

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

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值