题目传送门

分析:

现在我们需要求:

\(~~~~\sum_{i=1}^{n}\sum_{j=1}^{n}lcm(i,j)\)

\(=\sum_{i=1}^{n}\sum_{j=1}^{n}\frac{i ~\cdot ~j}{gcd(i,j)}\)

\(=\sum_{d=1}^{n}d\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum_{j=1}^{\lfloor\frac{n}{d}\rfloor}i\cdot j \cdot [gcd(i,j)=1]\)

这一步推导可以看做是在枚举最小公因数,然后找除去最小公因数后的互质的两个数

又因为一个有趣的结论:

\(~~~~\sum_{i=1}^{n}i\cdot [gcd(i,n)=1]=\frac{n~\cdot ~\varphi(n)}{2}\)

对于一个数x,如果它和n互质,那么n-x也会和n互质,满足对称性

可以用类似等差数列求和的方法求和

我们继续推式子:

\(~~~~\sum_{d=1}^{n}d\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum_{j=1}^{\lfloor\frac{n}{d}\rfloor}i\cdot j \cdot [gcd(i,j)=1]\)

\(=\sum_{d=1}^{n}d\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}i^2\varphi(i)\)

前半截可以O(sqrt(n))算

我们现在要考虑如何快速求:

\(~~~~\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}i^2\varphi(i)\)

首先知道一个公式:

\(~~~~\sum_{i|n}\varphi(i)=n\)

然后:

\(~~~~\sum_{i=1}^{n}\sum_{d|i}\varphi(d)=\frac{n\cdot (n+1)}{2}\)

所以:

\(~~~~\sum_{i=1}^{n}\sum_{d|i}\varphi(d)\cdot i^2=\frac{n^2\cdot (n+1)^2}{4}\)

\(~~~~\sum_{i=1}^{n}\sum_{d|i}\varphi(d)\cdot d^2\cdot (\frac{i}{d})^2=\frac{n^2\cdot (n+1)^2}{4}\)

变换一下:

\(~~~~\sum_{i=1}^{n}\sum_{d=1}^{\lfloor\frac{n}{i}\rfloor}\varphi(d)\cdot d^2\cdot i^2=\frac{n^2\cdot (n+1)^2}{4}\)

\(~~~~\sum_{i=1}^{n}sum(\lfloor\frac{n}{i}\rfloor)\cdot i^2=\frac{n^2\cdot (n+1)^2}{4}\)

将i=1分离出来,得到:

\(~~~~sum(n)=\frac{n^2\cdot (n+1)^2}{4}-\sum_{i=2}^{n}sum(\lfloor\frac{n}{i}\rfloor)\cdot i^2\)

然后也可以O(sqrt(n))算了

整体复杂度就是O(n^(2/3))

辣鸡式子推半天,写的时候有几个点没开long long调半天。。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
#include<iostream>
#include<map>
#include<string> #define maxn 5000005
#define INF 0x3f3f3f3f
#define MOD 1000000007
#define inv2 500000004
#define inv4 250000002
#define inv6 166666668 using namespace std; inline long long getint()
{
long long num=0,flag=1;char c;
while((c=getchar())<'0'||c>'9')if(c=='-')flag=-1;
while(c>='0'&&c<='9')num=num*10+c-48,c=getchar();
return num*flag;
} long long N;
long long pri[maxn],cnt,np[maxn],phi[maxn];
long long sum[maxn];
map<long long,long long>M;
long long ans; inline void init()
{
phi[1]=1;
for(int i=2;i<maxn;i++)
{
if(!np[i])pri[++cnt]=i,phi[i]=i-1;
for(int j=1;j<=cnt&&i*pri[j]<maxn;j++)
{
np[i*pri[j]]=1;
if(i%pri[j]==0){phi[i*pri[j]]=phi[i]*pri[j];break;}
phi[i*pri[j]]=phi[i]*(pri[j]-1);
}
}
for(int i=1;i<maxn;i++)sum[i]=(phi[i]*i%MOD*i%MOD+sum[i-1])%MOD;
} inline long long cal(long long x)
{return x*(x+1)%MOD*(2*x+1)%MOD*inv6%MOD;} inline long long getans(long long x)
{
if(x<maxn)return sum[x];
if(M.count(x))return M[x];
long long tmp=x%MOD;
long long num=tmp*tmp%MOD*(tmp+1)%MOD*(tmp+1)%MOD*inv4%MOD;
for(long long i=2,j;i<=x;i=j+1)
{
j=x/(x/i);
num+=MOD-(cal(j%MOD)-cal((i-1)%MOD)+MOD)%MOD*getans(x/i)%MOD;
num%=MOD;
}
return M[x]=num;
} int main()
{
N=getint();![](https://img2018.cnblogs.com/blog/1891964/202001/1891964-20200113163955776-377036624.jpg)
init();
for(long long i=1,j;i<=N;i=j+1)
{
j=N/(N/i);
ans+=(i+j)%MOD*((j-i+1)%MOD)%MOD*inv2%MOD*getans(N/i)%MOD;
ans%=MOD;
}
printf("%lld\n",ans);
}

51Nod 1238 最小公倍数之和V3的更多相关文章

  1. 51nod 1238 最小公倍数之和 V3

    51nod 1238 最小公倍数之和 V3 求 \[ \sum_{i=1}^N\sum_{j=1}^N lcm(i,j) \] \(N\leq 10^{10}\) 先按照套路推一波反演的式子: \[ ...

  2. 51NOD 1238 最小公倍数之和 V3 [杜教筛]

    1238 最小公倍数之和 V3 三种做法!!! 见学习笔记,这里只贴代码 #include <iostream> #include <cstdio> #include < ...

  3. 51nod 1238 最小公倍数之和 V3 【欧拉函数+杜教筛】

    首先题目中给出的代码打错了,少了个等于号,应该是 G=0; for(i=1;i<=N;i++) for(j=1;j<=N;j++) { G = (G + lcm(i,j)) % 10000 ...

  4. 51Nod 1238 - 最小公倍数之和 V3(毒瘤数学+杜教筛)

    题目 戳这里 推导 ∑i=1n∑j=1nlcm(i,j)~~~\sum_{i=1}^{n}\sum_{j=1}^{n}lcm(i,j)   ∑i=1n​∑j=1n​lcm(i,j) =∑i=1n∑j= ...

  5. 【51nod】1238 最小公倍数之和 V3 杜教筛

    [题意]给定n,求Σi=1~nΣj=1~n lcm(i,j),n<=10^10. [算法]杜教筛 [题解]就因为写了这个非常规写法,我折腾了3天…… $$ans=\sum_{i=1}^{n}\s ...

  6. 51 NOD 1238 最小公倍数之和 V3

    原题链接 最近被51NOD的数论题各种刷……(NOI快到了我在干什么啊! 然后发现这题在网上找不到题解……那么既然A了就来骗一波访问量吧…… (然而并不怎么会用什么公式编辑器,写得丑也凑合着看吧…… ...

  7. 51 Nod 1238 最小公倍数之和 V3 杜教筛

    题目链接:http://www.51nod.com/Challenge/Problem.html#!#problemId=1238 题意:求$\sum_{i=1}^{n}\sum_{j=1}^{n}l ...

  8. [51Nod 1238] 最小公倍数之和 (恶心杜教筛)

    题目描述 求∑i=1N∑j=1Nlcm(i,j)\sum_{i=1}^N\sum_{j=1}^Nlcm(i,j)i=1∑N​j=1∑N​lcm(i,j) 2<=N<=10102<=N ...

  9. 【学术篇】51nod 1238 最小公倍数之和

    这是一道杜教筛的入(du)门(liu)题目... 题目大意 求 \[ \sum_{i=1}^n\sum_{j=1}^nlcm(i,j) \] 一看就是辣鸡反演一类的题目, 那就化式子呗.. \[ \s ...

随机推荐

  1. vue-learning:7-template-v-bind-with-class-and-style

    绑定元素样式的指令v-bind:class 和v-bind:style 在HTML元素结构中,class和style特性(attribute)是非常突出的,可以为元素添加样式属性(property). ...

  2. python常见关键字的使用

    常见关键字 在循环中常见的关键字使用方法 continue:结束本次循环,继续执行下一次循环 break:跳出一个循环或者结束一个循环 例 使用用户名密码登录(有三次机会)count=0while c ...

  3. 第四阶段:1.从零打造一款社区web产品

    ---恢复内容开始--- 熟人关系:微信 陌生人关系:微博 1.把各种竞品罗列起来形成一个分析池.分析其目标用户是哪些.这些产品满足了用户什么需求.可以从时间角度分析趋势.针对每一类竞品画一个商业模式 ...

  4. 机器学习——Bagging与随机森林算法及其变种

    Bagging算法:  凡解:给定M个数据集,有放回的随机抽取M个数据,假设如此抽取3组,3组数据一定是有重复的,所以先去重.去重后得到3组数据,每组数据量分别是s1,s2,s3,然后三组分别训练组合 ...

  5. Kafka学习笔记(四)—— API使用

    Kafka学习笔记(四)-- API使用 1.Producer API 1.1 消息发送流程 Kafka的Producer发送消息采用的是异步发送的方式.在消息发送的过程中,涉及到了两个线程--mai ...

  6. nginx和keeplive实现负载均衡高可用

    一. Keeplive服务介绍 Keeplive期初是专门为LVS设计的,专门用来监控LVS集群系统中各个服务节点的状态,后来又加入VRRP的功能,因此除了配合LVS服务以外,也可以作为其他服务(ng ...

  7. 原生JS数组方法实现(一)————push()、unshift()、pop()和shift()

    push 向数组末尾添加一个或多个元素,并返回数组新的长度 ```javascript function push(){ for(let i=0;i<arguments.length;i++){ ...

  8. 【证明与推广与背诵】Matrix Tree定理和一些推广

    [背诵手记]Matrix Tree定理和一些推广 结论 对于一个无向图\(G=(V,E)\),暂时钦定他是简单图,定义以下矩阵: (入)度数矩阵\(D\),其中\(D_{ii}=deg_i\).其他= ...

  9. 【题解】CF986E Prince's Problem(树上差分+数论性质)

    [题解]CF986E Prince's Problem(树上差分+数论性质) 题目大意: 给定你一棵树,有点权\(val_i\le 10^7\).现在有\(m\)组询问给定参数\(x,y,w\)问你对 ...

  10. $CH$ $0x50$ & $0x51$ 做题记录

    [X]$Mr.Young's\ Picture\ Permutations$ 前面这儿写了挺多道辣,,,懒得写辣$QAQ$ (后面所有同上都是同这个$QwQ$ [X]$LCIS$ 做过了,看这儿 $u ...