hdu5942 Just a Math Problem

题意

f(i) 为i的不同质因子的个数,求 2f(i)[1<=i<=n]

思路

如果 i=pc11pc22...pcf(i)f(i) ,则满足 PQ=igcd(P,Q)=1 (P,Q) 有序对的个数为 2f(i) ,因为每个 pckk 要么属于 P ,要么属于Q。那么原题转化为求满足 1<=PQ<=ngcd(P,Q)=1 (P,Q) 有序对的个数,即

[1<=P,Q<=n][1<=PQ<=n][gcd(P,Q)=1]
[1<=P,Q<=n]=[1<=P<Q<=n]+[1<=Q<P<=n]+[1<=P=Q<=n]
所以
[1<=P,Q<=n][1<=PQ<=n][gcd(P,Q)=1]=1+2[1<=P<Q<=n][1<=PQ<=n][gcd(P,Q)=1]
记它为S。
枚举 1<=P<n ,则 P<Q<=nP ,则
S=[1<=P<n][P<Q<=np][gcd(P,Q)=1]
=[1<=P<n]([1<=Q<=np][gcd(P,Q)=1][1<=Q<=P][gcd(P,Q)=1])
,所以问题就是要求1~n与P互质的数的个数。
假设 P 的质因子集合为T={p1,p2,...pt},则“与P互质”等价于“不是 p1 的倍数,且不是 p2 的倍数,…,且不是 pt 的倍数”,取反即“是 p1 的倍数,或是 p2 的倍数,…,或是 pt 的倍数”,容斥原理即可求之。

代码
#include <bits/stdc++.h>
using namespace std;
#ifdef LOCAL
#include "local.h"
#endif // LOCAL
typedef long long ll;
const ll mod = 1e9 + 7;
const ll N = 1e6 + 7;
vector<ll> prime, fac[N];
ll minfac[N], phi[N];
bool used[N];
void init() {
  phi[1] = 1;
  for (ll i = 2; i < N; i++) {
    if (!used[i]) {
      minfac[i] = i;
      phi[i] = i - 1;
      prime.push_back(i);
    }
    for (ll j = 0; j < prime.size(); j++) {
      if (prime[j] * i >= N) break;
      used[prime[j] * i] = 1;
      minfac[prime[j] * i] = prime[j];
      if (i % prime[j] == 0) {
        phi[prime[j] * i] = prime[j] * phi[i];
        break;
      }
      else {
        phi[prime[j] * i] = (prime[j] - 1) * phi[i];
      }
    }
  }
  for (ll i = 2; i < N; i++) {
    ll x = i;
    while (x > 1) {
      ll buf = minfac[x];
      fac[i].push_back(buf);
      while (x % buf == 0) x /= buf;
    }
  }
}
inline ll lowbit(ll x) {
  return x & -x;
}
ll Log2[N];
ll calc(vector<ll> &fac, ll n) {
  ll ans = n, sz = fac.size();
  ll mul[1 << sz] = {1};
  for (ll i = 1; i < (1 << sz); i++) {
    ll lb = lowbit(i);
    mul[i] = mul[i - lb] * -fac[Log2[lb]];
    ans = (ans + n / mul[i] + mod * mod) % mod;
  }
  return ans;
}
ll solve(ll n) {
  for (ll i = 0; (1 << i) < N; i++) {
    Log2[1 << i] = i;
  }
  ll ans = 0;
  for (ll p = 1; p * p < n; p++) {
    ans = (ans + calc(fac[p], n / p) - phi[p] + mod) % mod;
  }
  return (ans * 2 + 1) % mod;
}
int _, cas;
ll n;
int main() {
  init();
  cin >> _;
  while (_--) {
    printf("Case #%d: ", ++cas);
    cin >> n;
    cout << solve(n) << endl;
  }
  return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值