Codeforces Round #729 (Div. 2) C. Strange Function

最新推荐文章于 2023-01-24 01:31:17 发布

TURNINING

最新推荐文章于 2023-01-24 01:31:17 发布

阅读量331

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/TURNINING/article/details/118529598

版权

数论专栏收录该内容

3 篇文章 0 订阅

订阅专栏

本文介绍了一种高效的算法来解决计算函数f(n)的问题，即第一个不为n的因子的数的和。通过分析因子分布和最小公倍数，作者给出了求解f(i)个数的方法，并提供了C++代码实现。核心思想是利用lcm和因子数量的计算来减少大规模数据的运算复杂度。

摘要由CSDN通过智能技术生成

传送门

题意：函数f(n) = 第一个不为n的因子的数，求f(i)的和。

思路：
n这么大肯定考虑算贡献。关键是求f(i) = x, i 的个数
显然，如果f(i) = x，那么i肯定有因子1,2,3…,(i-1),lcm(1到i-1)就为含有这些因子的最小值，n / lcm(1到i-1)就为1到n中含有这些因子的数的个数，f(i) = x要求x不为i的因子，含有1,2,…i这些因子的最小的数为lcm(1到i)，含有n / lcm(1到i)个。
相减即为要求的个数。

#include<bits/stdc++.h>
using namespace std;
 
typedef long long ll;
typedef unsigned long long ull;
 
typedef pair<int, int> PI;
typedef pair<double, double> PD;
 
const int maxn = 1e3 + 10;
const int max_log_v = 22;
const int INF = 0x3f3f3f3f;
const int mod = 1e9 + 7;
const double eps = 1e-7;
const ull B = 100000007;

ll n;

ll gcd(ll a, ll b) {
    if(b == 0) return a;
    return gcd(b, a % b);
}

ll lcm(ll a, ll b) {
    return a * b / gcd(a, b);
}

void solve() {
    scanf("%lld", &n);
    ll res = 0;
    ll l = 1;
    for(ll i = 2; l <= n; i++) {
        res = (res % mod + i % mod * (n/l - n/lcm(l, i)) % mod) % mod;
        l = lcm(l, i);
    }
    printf("%lld\n", res);
}
 
 
int main() {
    int t = 1; scanf("%d", &t);
    while(t--) {
        solve();
    }
    return 0;
}