参考:http://blog.csdn.net/u014610830/article/details/49493279

这道题做起来感觉非常奇怪啊……头一次见把mu推出来再推没了的……

\[\sum_{i=a}^{b}lcm(i,b)
\]

\[\sum_{i=a}^{b}\frac{i*b}{gcd(i,b)}
\]

\[\sum_{d|b}\sum_{i=a}^{b}[gcd(i,b)==d]\frac{ib}{d}
\]

\[\sum_{d|b}\sum_{i=\left \lfloor \frac{a}{d} \right \rfloor}^{\left \lfloor \frac{b}{d} \right \rfloor}[gcd(i,b)==1]\frac{idb}{d}
\]

\[b\sum_{d|b}\sum_{i=\left \lfloor \frac{a}{d} \right \rfloor}^{\left \lfloor \frac{b}{d} \right \rfloor}[gcd(i,b)==1]i
\]

\[b\sum_{d|b}\sum_{i=\left \lfloor \frac{a}{d} \right \rfloor}^{\left \lfloor \frac{b}{d} \right \rfloor}i\sum_{k|gcd(i,\frac{b}{d})}\mu(k)
\]

\[b\sum_{d|b}\sum_{k|\frac{b}{d}}\mu(k)\sum_{i=\left \lfloor \frac{a}{d} \right \rfloor}^{\left \lfloor \frac{b}{d} \right \rfloor,k|i}i
\]

\[b\sum_{d|b}\sum_{k|\frac{b}{d}}\mu(k)k\sum_{i=\left \lfloor \frac{a}{dk} \right \rfloor}^{\left \lfloor \frac{b}{dk} \right \rfloor}i
\]

\[f(a,b)=\sum_{i=a}^{b}i
\]

\[b\sum_{d|b}\sum_{k|\frac{b}{d}}\mu(k)kf(\left \lfloor \frac{a}{dk} \right \rfloor,\left \lfloor \frac{b}{dk} \right \rfloor)
\]

枚举dk

\[b\sum_{d|b}f(\frac{a}{d},\frac{b}{d})\sum_{k|d}\mu(k)k
\]

其中f可以\( O(1) \)算,设\( g(x)=\sum_{k|x}\mu(k)k \),容易发现这是积性的,显然设p为质数,\( g(p)=1-p,g(p^k)=1-p \),于是这样就把mu推掉了

时间复杂度很不科学但是还是能过

#include<iostream>
#include<cstdio>
#include<vector>
using namespace std;
const int N=100005,m=100000,inv2=500000004,mod=1e9+7;
int T,n,q[N],tot,a,b,p[N],cnt,len[N];
bool v[N];
vector<int>fac;
void dfs(int re,int l)
{
if(l>cnt)
{
fac.push_back(re);
return;
}
int tmp=1;
dfs(re,l+1);
for(int i=1;i<=len[l];i++)
{
tmp*=p[l];
dfs(re*tmp,l+1);
}
}
int main()
{
v[1]=1;
for(int i=2;i<=m;i++)
{
if(!v[i])
q[++tot]=i;
for(int j=1;j<=tot&&i*q[j]<=m;j++)
{
int k=i*q[j];
v[k]=1;
if(i%q[j]==0)
break;
}
}
scanf("%d",&T);
while(T--)
{
scanf("%d%d",&a,&b);
cnt=0;
int now=b;
for(int i=1;i<=tot&&q[i]*q[i]<=now;i++)
{
if(now%q[i]==0)
{
p[++cnt]=q[i];
len[cnt]=0;
}
while(now%q[i]==0)
{
now/=q[i];
len[cnt]++;
}
}
if(now>1)
{
p[++cnt]=now;
len[cnt]=1;
}
fac.clear();
dfs(1,1);//cout<<cnt<<endl;
long long ans=0ll;
for(int i=0;i<fac.size();i++)
{
int v=fac[i];
long long tt=a+v-1,tmp,re=1ll;
tmp=(((tt/v+b/v)%mod)*((b/v-tt/v+1+mod)%mod)%mod)*inv2%mod;
for(int j=1;j<=cnt;j++)
if(v%p[j]==0)
re=re*((1-p[j]+mod)%mod)%mod;
ans=((ans+re*tmp%mod)%mod+mod)%mod;
}
printf("%lld\n",ans*b%mod);
}
return 0;
}

51nod 1190 最小公倍数之和 V2【莫比乌斯反演】的更多相关文章

  1. 51nod 1190 最小公倍数之和 V2

    给出2个数a, b,求LCM(a,b) + LCM(a+1,b) + .. + LCM(b,b). 例如:a = 1, b = 6,1,2,3,4,5,6 同6的最小公倍数分别为6,6,6,12,30 ...

  2. 51nod 1222 最小公倍数计数【莫比乌斯反演】

    参考:https://www.cnblogs.com/SilverNebula/p/7045199.html 所是反演其实反演作用不大,又是一道做起来感觉诡异的题 转成前缀和相减的形式 \[ \sum ...

  3. 51nod1238 最小公倍数之和 V3(莫比乌斯反演)

    题意 题目链接 Sol 不想打公式了,最后就是求一个 \(\sum_{i=1}^n ig(\frac{N}{i})\) \(g(i) = \sum_{i=1}^n \phi(i) i^2\) 拉个\( ...

  4. 51nod 1220 约数之和【莫比乌斯反演+杜教筛】

    首先由这样一个式子:\( d(ij)=\sum_{p|i}\sum_{q|j}[gcd(p,q)==1]\frac{pj}{q} \)大概感性证明一下吧我不会证 然后开始推: \[ \sum_{i=1 ...

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

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

  6. 【CJOJ2512】gcd之和(莫比乌斯反演)

    [CJOJ2512]gcd之和(莫比乌斯反演) 题面 给定\(n,m(n,m<=10^7)\) 求 \[\sum_{i=1}^n\sum_{j=1}^mgcd(i,j)\] 题解 首先把公因数直 ...

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

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

  8. 51nod 1363 最小公倍数之和 ——欧拉函数

    给出一个n,求1-n这n个数,同n的最小公倍数的和.例如:n = 6,1,2,3,4,5,6 同6的最小公倍数分别为6,6,6,12,30,6,加在一起 = 66. 由于结果很大,输出Mod 1000 ...

  9. [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 ...

随机推荐

  1. CSY版最大团,速度快一倍

    #include <bits/stdc++.h> using namespace std; #define REP(i, n) for(int i(0); i < (n); ++i) ...

  2. React学习及实例开发(一)——开始

    本文基于React v16.4.1 初学react,有理解不对的地方,欢迎批评指正^_^ 一.构建一个新项目 1.命令行运行如下命令,构建一个新的react项目 npm install -g crea ...

  3. 某考试 T1 table

    我们把每一行看成一个多项式 f[i][1] + f[i][2] * x + f[i][3] * x^2 + ..... + f[i][n] * x^(n-1) 的话,一行转移到下一行就相当于乘上一个{ ...

  4. spring mvc技术

                                                     spring mvc之访问路径        1. @RequestMapping这个注解  在实际项 ...

  5. Extjs.panel.Panel赋值的问题

    初学extjs,很是不爽.也是只有初学者才犯的错误,发出来以免再犯. 先创建一个panel var panel1 = Ext.create('Ext.panel.Panel', { id: 'p1', ...

  6. MySQL基础笔记(二) 完整性约束

    我们知道,一种数据模型必须包含三个基本的部分: 构造机制(数据结构):主要描述数据的类型.内容.性质以及数据间的联系等. 运算机制(数据操作):主要描述在相应的数据结构上的操作类型和操作方式. 约束机 ...

  7. Android耳机线控具体解释,蓝牙耳机button监听(仿酷狗线控效果)

    转载请注明出处:http://blog.csdn.net/fengyuzhengfan/article/details/46461253 当耳机的媒体按键被单击后.Android系统会发出一个广播.该 ...

  8. 一个很小的C++写的MVC的例子

    #include<iostream> #include<vector> //get namespace related stuff using std::cin; using ...

  9. POJ 1927 Area in Triangle(计算几何)

    Area in Triangle 博客原文地址:http://blog.csdn.net/xuechelingxiao/article/details/40707691 题目大意: 给你一个三角形的三 ...

  10. 2015/12/25 ① 图灵测试 ② 安装jdk出现的问题 ③ 配置环境变量

    ①图灵测试 1,解释 图灵测试一词来源于计算机科学和密码学的先驱阿兰·麦席森·图灵写于1950年的一篇论文<计算机器与智能>.阿兰·麦席森·图灵1950年设计出这个测试,其内容是,如果电脑 ...