洛谷p3768 简单的数学题

题目大意

i=1nj=1nijgcd(i,j)
对质数p取膜的结果

n<=10^10 5*10^8<=p<=1.1*10^9

传送门https://www.luogu.org/problem/show?pid=3768

分析

看到这种两层求和带gcd的式子 可以考虑直接用莫比乌斯反演求解。
但注意到i j的上界都是n,所以用 ϕ 推导公式也可以解决这题。(感觉可能会简便一些)

i=1nj=1nijgcd(i,j)

=d=1nd3i=1ndj=1ndij[ij]

=d=1nd3G(nd)

右边的式子是关于 nd 的 可以看成一个函数
如果能够快速计算出G函数 那么这个式子可以分块快速求出
G(n)=i=1nj=1nij[ij]

=i=1nij=1nj[ij]

=2i=1nij=1ij[ij]1

=2i=1ni(iϕ(i)+[i==1]2)1

=i=1n(i2ϕ(i)+[i==1])1

=i=1ni2ϕ(i)

(嗯 看起来很简单吧..)

G(n)显然是一个积性函数的前缀和
如果预处理出G函数 那么可以设计一个O(n)的算法 可以通过60%的数据。

G(n)=ni=1f(i)
考虑将 ϕ 通过狄利克雷卷积消去
已知 ϕId0=Id1

假设 fh=T
T(n)=d|nf(d)h(nd)=d|nd2ϕ(d)h(nd)=d|n(d2h(nd))ϕ(d)
h=Id2 d2h(nd)n2phi(d)n
卷积即可得到 T=Id3

那么就可以用杜教筛快速计算出G(n)了

G(n)=n3i=2ni2G(ni)

总时间复杂度 O(n23)

#include<stdio.h>
#include<math.h>
#include<string.h>
#include<stdlib.h>
#include<map>
#include<set>
#include<algorithm>
#include<queue>
#include<time.h>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
#define fe(x) for(int ii=be[x];ii;ii=e[ii].ne)
using namespace std;
const int mn=20000000;
int p,pri[mn/10],pt,mm;
long long n,p2,p6,ans,g[mn+20],mg[100000];
bool vis[100000],isp[mn+20];
inline long long qpow(int a,int n){
    long long x=a,ans=1;
    while (n){
        if (n&1) ans=ans*x%p;
        x=x*x%p;
        n>>=1;
    }
    return ans;
}
inline long long sum2(long long n){
    n%=p;
    return (n+n+1)%p*(n+1)%p*n%p*p6%p;
}
inline long long sum3(long long n){
    n%=p;
    n=n*(n+1)%p*p2%p;
    return n*n%p;
}
inline long long G(long long x){
    if (x<=mm) return g[x];
    int ms=n/x;
    if (vis[ms]) return mg[ms];
    vis[ms]=1;
    long long i=2,j,te;
    long long &ans=mg[ms];
    ans=sum3(x);
    while (i<=x){
        te=x/i;
        j=x/te;
        ans-=(sum2(j)-sum2(i-1))*G(te)%p;
        i=j+1;
    }
    return ans%=p;
}
int main(){
    scanf("%d%lld",&p,&n);
    mm=mn;if (n<mm) mm=n;
    g[1]=1;
    fo(i,2,mm){
        if (!isp[i]){
            g[i]=i-1;
            pri[++pt]=i;
        }
        fo(j,1,pt){
            long long te=i*(long long)pri[j];
            if (te>mn) break;
            isp[te]=1;
            if (i%pri[j]==0){
                g[te]=g[i]*pri[j];
                break;
            }
            g[te]=g[i]*(pri[j]-1);
        }
    }
    fo(i,2,mm) g[i]=(g[i-1]+g[i]*i%p*i)%p;
    p2=qpow(2,p-2);p6=qpow(6,p-2);
    long long i=1,j,te;
    while (i<=n){
        te=n/i;
        j=n/te;
        ans+=(sum3(j)-sum3(i-1))*G(te)%p;
        i=j+1;
    }
    ans=(ans%p+p)%p;
    printf("%lld\n",ans);
    return 0;
}

洛谷P2791是一道关于幼儿园篮球的题目。题目描述如下: 在幼儿园里,小朋友们正在进行篮球比赛。每个小朋友都有一个篮球,他们按照顺序依次投篮。每个小朋友投篮时,可以选择将篮球投给左边的小朋友或者右边的小朋友。每个小朋友投篮得分的规则如下: 1. 如果一个小朋友左右两边的小朋友都没有投篮过,那么他的得分为1; 2. 如果一个小朋友左边的小朋友投篮过,但右边的小朋友没有投篮过,那么他的得分为左边小朋友的得分加1; 3. 如果一个小朋友右边的小朋友投篮过,但左边的小朋友没有投篮过,那么他的得分为右边小朋友的得分加1; 4. 如果一个小朋友左右两边的小朋友都投篮过,那么他的得分为左边小朋友和右边小朋友得分的最大值加1。 现在给定每个小朋友投篮的顺序,请你计算每个小朋友的得分。 例如,给定投篮顺序为[1, 0, 1, 0, 1],则第一个小朋友的得分为1,第二个小朋友的得分为2,第三个小朋友的得分为1,第四个小朋友的得分为2,第五个小朋友的得分为1。 你可以通过编写程序来解决这个问题。具体的解题思路可以参考以下步骤: 1. 创建一个数组scores,用来存储每个小朋友的得分; 2. 遍历投篮顺序数组,对于每个小朋友,根据上述规则计算他的得分,并将得分存入scores数组中; 3. 最后输出scores数组即可。 希望以上解答对你有帮助!如果你还有其他问题,请继续提问。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值