HDU 3501 更相减损术 + 欧拉函数

题意

传送门 HDU 3501

题解

求小于 N N N 且与 N N N 不互素的正整数之和,考虑小于 N N N 的正整数之和减去小于 N N N 且与 N N N 互素的正整数之和。根据更相减损术,若 g c d ( N , x ) = 1 gcd(N,x)=1 gcd(N,x)=1,则有 g c d ( N , N − x ) gcd(N,N-x) gcd(N,Nx)=1,那么考虑成对的与 N N N 互素的数,则答案为 ( N − 1 ) × N / 2 − e u l e r ( N ) ∗ N / 2 (N-1)\times N/2 - euler(N)*N/2 (N1)×N/2euler(N)N/2

#include <algorithm>
#include <cmath>
#include <cstdio>
#include <cstring>
using namespace std;
typedef long long ll;
const int mod = 1000000007;
int N;

int euler_phi(int n)
{
    int res = n;
    for (int i = 2; i * i <= n; ++i)
    {
        if (n % i == 0)
        {
            res = res / i * (i - 1);
            while (n % i == 0)
                n /= i;
        }
    }
    if (n != 1)
        res = res / n * (n - 1);
    return res;
}

int main()
{
    while (~scanf("%d", &N) && N)
    {
        ll res = 1LL * N * (N - 1) / 2;
        res -= 1LL * euler_phi(N) * N / 2;
        printf("%lld\n", (res % mod + mod) % mod);
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值