Bzoj4916: 神犇和蒟蒻

题面

传送门

Sol

第一问puts(“1”)
第二问, φ(i2)=iφ(i)
ϕ(n)=ni=1iφ(i) 根据杜教筛推的式子

g(1)ϕ(n)=i=1nd|ig(d)(id)φ(id)i=2ng(d)ϕ(nd)

g(i)=i 减号前面原式 =ni=1id|iφ(id)
i|nφ(i)=n 所以就是 ni=1i2=n(n+1)(2n+1)6
跑杜教筛即可

# include <bits/stdc++.h>
# define RG register
# define IL inline
# define Fill(a, b) memset(a, b, sizeof(a))
using namespace std;
typedef long long ll;
const int _(1e6 + 1), Zsy(1e9 + 7), yyb(166666668);

IL ll Read(){
    RG ll x = 0, z = 1; RG char c = getchar();
    for(; c < '0' || c > '9'; c = getchar()) z = c == '-' ? -1 : 1;
    for(; c >= '0' && c <= '9'; c = getchar()) x = (x << 1) + (x << 3) + (c ^ 48);
    return x * z;
}

int prime[_], num;
ll phi[_];
map <int, int> Phi;
bool isprime[_];

IL void Prepare(){
    isprime[1] = 1; phi[1] = 1;
    for(RG int i = 2; i < _; ++i){
        if(!isprime[i]){  prime[++num] = i; phi[i] = (i - 1);  }
        for(RG int j = 1; j <= num && i * prime[j] < _; ++j){
            isprime[i * prime[j]] = 1;
            if(i % prime[j]) phi[i * prime[j]] = 1LL * phi[i] * (prime[j] - 1) % Zsy;
            else{  phi[i * prime[j]] = 1LL * phi[i] * prime[j] % Zsy; break;  }
        }
    }
    for(RG int i = 2; i < _; ++i) phi[i] = 1LL * phi[i] * i % Zsy, (phi[i] += phi[i - 1]) %= Zsy;
}

IL ll Sumphi(RG ll n){
    if(n < _) return phi[n];
    if(Phi[n]) return Phi[n];
    RG ll ans = n * (n + 1) % Zsy * (2 * n + 1) % Zsy * yyb % Zsy;
    for(RG ll i = 2, j; i <= n; i = j + 1){
        j = n / (n / i);
        ans -= 1LL * (j + i) * (j - i + 1) / 2 % Zsy * Sumphi(n / i) % Zsy;
        ans = (ans + Zsy) % Zsy;
    }
    return Phi[n] = ans;
}

int main(RG int argc, RG char* argv[]){
    Prepare(); RG ll n = Read();
    printf("1\n%lld\n", Sumphi(n));
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值