51nod 1190 最小公倍数之和 V2【莫比乌斯反演】
参考:http://blog.csdn.net/u014610830/article/details/49493279
这道题做起来感觉非常奇怪啊……头一次见把mu推出来再推没了的……
\]
\]
\]
\]
\]
\]
\]
\]
设
\]
\]
枚举dk
\]
其中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【莫比乌斯反演】的更多相关文章
- 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 ...
- 51nod 1222 最小公倍数计数【莫比乌斯反演】
参考:https://www.cnblogs.com/SilverNebula/p/7045199.html 所是反演其实反演作用不大,又是一道做起来感觉诡异的题 转成前缀和相减的形式 \[ \sum ...
- 51nod1238 最小公倍数之和 V3(莫比乌斯反演)
题意 题目链接 Sol 不想打公式了,最后就是求一个 \(\sum_{i=1}^n ig(\frac{N}{i})\) \(g(i) = \sum_{i=1}^n \phi(i) i^2\) 拉个\( ...
- 51nod 1220 约数之和【莫比乌斯反演+杜教筛】
首先由这样一个式子:\( d(ij)=\sum_{p|i}\sum_{q|j}[gcd(p,q)==1]\frac{pj}{q} \)大概感性证明一下吧我不会证 然后开始推: \[ \sum_{i=1 ...
- 51nod 1238 最小公倍数之和 V3
51nod 1238 最小公倍数之和 V3 求 \[ \sum_{i=1}^N\sum_{j=1}^N lcm(i,j) \] \(N\leq 10^{10}\) 先按照套路推一波反演的式子: \[ ...
- 【CJOJ2512】gcd之和(莫比乌斯反演)
[CJOJ2512]gcd之和(莫比乌斯反演) 题面 给定\(n,m(n,m<=10^7)\) 求 \[\sum_{i=1}^n\sum_{j=1}^mgcd(i,j)\] 题解 首先把公因数直 ...
- 51NOD 1238 最小公倍数之和 V3 [杜教筛]
1238 最小公倍数之和 V3 三种做法!!! 见学习笔记,这里只贴代码 #include <iostream> #include <cstdio> #include < ...
- 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 ...
- [51Nod 1238] 最小公倍数之和 (恶心杜教筛)
题目描述 求∑i=1N∑j=1Nlcm(i,j)\sum_{i=1}^N\sum_{j=1}^Nlcm(i,j)i=1∑Nj=1∑Nlcm(i,j) 2<=N<=10102<=N ...
随机推荐
- CSY版最大团,速度快一倍
#include <bits/stdc++.h> using namespace std; #define REP(i, n) for(int i(0); i < (n); ++i) ...
- React学习及实例开发(一)——开始
本文基于React v16.4.1 初学react,有理解不对的地方,欢迎批评指正^_^ 一.构建一个新项目 1.命令行运行如下命令,构建一个新的react项目 npm install -g crea ...
- 某考试 T1 table
我们把每一行看成一个多项式 f[i][1] + f[i][2] * x + f[i][3] * x^2 + ..... + f[i][n] * x^(n-1) 的话,一行转移到下一行就相当于乘上一个{ ...
- spring mvc技术
spring mvc之访问路径 1. @RequestMapping这个注解 在实际项 ...
- Extjs.panel.Panel赋值的问题
初学extjs,很是不爽.也是只有初学者才犯的错误,发出来以免再犯. 先创建一个panel var panel1 = Ext.create('Ext.panel.Panel', { id: 'p1', ...
- MySQL基础笔记(二) 完整性约束
我们知道,一种数据模型必须包含三个基本的部分: 构造机制(数据结构):主要描述数据的类型.内容.性质以及数据间的联系等. 运算机制(数据操作):主要描述在相应的数据结构上的操作类型和操作方式. 约束机 ...
- Android耳机线控具体解释,蓝牙耳机button监听(仿酷狗线控效果)
转载请注明出处:http://blog.csdn.net/fengyuzhengfan/article/details/46461253 当耳机的媒体按键被单击后.Android系统会发出一个广播.该 ...
- 一个很小的C++写的MVC的例子
#include<iostream> #include<vector> //get namespace related stuff using std::cin; using ...
- POJ 1927 Area in Triangle(计算几何)
Area in Triangle 博客原文地址:http://blog.csdn.net/xuechelingxiao/article/details/40707691 题目大意: 给你一个三角形的三 ...
- 2015/12/25 ① 图灵测试 ② 安装jdk出现的问题 ③ 配置环境变量
①图灵测试 1,解释 图灵测试一词来源于计算机科学和密码学的先驱阿兰·麦席森·图灵写于1950年的一篇论文<计算机器与智能>.阿兰·麦席森·图灵1950年设计出这个测试,其内容是,如果电脑 ...