Strange Function(数论，容斥，1600)

最新推荐文章于 2024-10-06 14:10:45 发布

Cambrain_

最新推荐文章于 2024-10-06 14:10:45 发布

阅读量126

点赞数

分类专栏：数学 codeforces 文章标签：算法

本文链接：https://blog.csdn.net/m0_63761643/article/details/127845686

版权

数学同时被 2 个专栏收录

23 篇文章 0 订阅

订阅专栏

codeforces

16 篇文章 0 订阅

订阅专栏

本文解析了Codeforces问题1542C中关于f(x)的计算方法，通过容斥原理和lcm的应用，展示了如何编写代码求解给定n时f(i)=x的个数及其对答案的贡献。关键在于利用lcm和整除性质进行优化。

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

题目链接：

Problem - 1542C - Codeforces

题目大意：

定义f(x),为无法整除x的最小正整数，给定n，求 $\sum f(i)$

思路：

打表未发现明显的规律，考虑到f(i)=x,则

i可以整除lcm(1,2,3,...,x-1),i不可整除lcm(1,2,3,...,x);

则1~n中f(i)=x的个数为：

cnt=n/lcm(1,2,3,...,x-1)-n/lcm(1,2,3,...x) (容斥原理)

对答案的贡献即为x*cnt.不难发现，当lcm(1,2,3,....x-1)大于n时cnt为0，对答案无贡献。因此x大概枚举到100即可.

代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> PII;
const int mod=1e9+7;
ll gcd(ll n,ll m)
{
    return m?gcd(m,n%m):n;
}
ll lcm(ll n,ll m)
{
    return n/gcd(n,m)*m;
}
void solve()
{
    ll n;
    cin>>n;
    ll res=0,now=1,cnt;
    for(int i=2;i<=n+1;i++){
        if(now>n) break;
        cnt=n/now;
        now=lcm(now,i);
        cnt-=n/now;
        res=(res+cnt%mod*i%mod)%mod;
    }
    cout<<res<<'\n';
}
int main()
{
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    cout.tie(nullptr);
    int t;
    cin>>t;
    while(t--) solve();
    return 0;
}