【codevs4355】王的对决(简单数论) 莫比乌斯反演
题目描述 Description
/*在神秘而遥远的T国度生活着一群神奇的物种——DD!
DD们有着独特的生活习惯,他们总是喜欢居住在tree上,而且只有上层住了DD下层才能再住。并且每层住的DD最多不超过2n2个。。。
这一天,DD跳下来了!(也就是说刚才的条件一点用都没有~_~)*/
他们遇到了一只自称为S的dd。但是本身性质为S的DD们不服,于是他们决定进行摔跤,摔跤获胜者成为王!
(由于DD和dd们的精力无限,所以每一只DD&&dd都会同另一部落的每一个人进行一次摔跤)两阵营分别有n,m只DD和dd,并且每阵营按照1~n和1~m进行体力值分配(体力值同序号相同)
进行比赛时,双方中体力值较大一方会减去另一方的体力值使之作为自己的新体力值,但由于DD和dd们实力雄厚,因此只有当双方的体力值相同时才会停止摔跤。如果此时双方的体力值是一个纯洁的数字,DD和dd会很高兴地退场,并且将这个数作为这场比赛的得分~~~,但如果只是一般的数字,他们会不开心地退场导致这场比赛没有得分。因为DD和dd作为T团队不可或缺的一员,所以获得的分数会加入T团队的总分(一场只有一个得分,不要以为你们是两个人,就可以得两遍分啊!!!)
输出T团队最后获得的总分,数据有多组询问。
注:纯洁的数字即最后的这个值没有次幂高于1的因子(SFN,无平方因子数)
一句话题意:i从1到n,j从1到m,如果gcd(i,j)为sfn,那么就加起来,求最终答案……
输入描述 Input Description
第一行一个正整数N表示一共N组数据,接下来N行每行两个正整数n,m
输出描述 Output Description
输出有N行,每行一个正整数,表示T团队的总分
样例输入 Sample Input
1
5 5
样例输出 Sample Output
33
数据范围及提示 Data Size & Hint
对于样例:
1 代表一组数据;
下面1行两个数5 5
1和1~5摔跤;2和1~5摔跤;3和1~5摔跤;4和1~5摔跤;5和1~5摔跤;
1和任何人摔跤返回的得分总是1,因为1是最小的。
2和1摔跤得1,2和2摔跤得2,2和3摔跤得1,2和4摔跤得2,2和5摔跤得1;
3和1摔跤得1;3和2摔跤得1;3和3摔跤得3;3和4摔跤得1;3和5摔跤得1;
4和1摔跤得1;4和2摔跤得2;4和3摔跤得1;4和4摔跤得4;5和4摔跤得1;
5和1摔跤得1;5和2摔跤得1;5和3摔跤得1;5和4摔跤得1;5和5摔跤得5;
最后把满足square free number的数加起来(4不是square free number)
30%的数据:N<=10,n<=m<=1000
60%的数据: N<=300,n<=m<=600000
100%的数据: N<=10000,n<=m<=1500000
请尽可能水分~
2015.5.1我们在tyvj出的比赛的题,Orzseavot&&rainheart
话说题面还是黑的我QAQ
设
f(n)= \begin{cases} n & \text {如果n是sfn} \\ 0 & \text {else} \\ \end{cases}
那么喜闻乐见,原式等于:
\sum_{i
=\sum_{d
其中F是f反演后的式子。
然后筛,把三种情况带进去看看就行了,推一下 F(pk) F(p^k)就行了…
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;typedef long long LL;const int SZ = 2000010;
const int MAXN = 1500000;bool vis[SZ];
int pri[SZ],f[SZ],F[SZ],a[SZ];
LL sum[SZ];void shai()
{a[1] = f[1] = F[1] = 1;for(int i = 2,tot = 0;i <= MAXN;i ++){if(!vis[i]) {pri[++ tot] = i;f[i] = i;F[i] = i - 1;a[i] = i;}for(int j = 1,m;j <= tot && (m = i * pri[j]) <= MAXN;j ++){vis[m] = 1;if(i % pri[j] == 0){f[m] = 0;F[m] = -F[i / a[i]] * f[a[i]];a[m] = a[i] * pri[j];break;}else{f[m] = f[i] == 0 ? 0 : m;F[m] = F[i] * F[pri[j]];a[m] = pri[j];}}}for(int i = 1;i <= MAXN;i ++)sum[i] = sum[i - 1] + F[i];
}LL ask(int n,int m)
{LL ans = 0;if(n > m) swap(n,m);for(int i = 1,r;i <= n;i = r + 1){r = min(n / (n / i),m / (m / i));ans += (sum[r] - sum[i - 1]) * (n / i) * (m / i);}return ans;
}int main()
{int T;scanf("%d",&T);shai();while(T --){int n,m;scanf("%d%d",&n,&m);printf("%lld\n",ask(n,m));}return 0;
}
【codevs4355】王的对决(简单数论) 莫比乌斯反演相关推荐
- 【数论+莫比乌斯反演】Exclusive Multiplication | gym103688E
[数论+莫比乌斯反演]Exclusive Multiplication | gym103688 前言 题意 思路 代码 前言 [ f ( x ) ] [f(x)] [f(x)] 表示艾弗森括号,若 f ...
- [数论]莫比乌斯反演1
索引 莫比乌斯反演1 定理 莫比乌斯反演2 证明 莫比乌斯反演3 技巧 前言 本篇内容全部为定理,无证明 定义 莫比乌斯函数的符号为\(\mu\),通俗的来讲 \[ \mu(n) = \left\{ ...
- 数论 —— 莫比乌斯反演
[反演] 假设我们手头有个数列 F,通过某种变换 H,可以得到函数 G.,即: 但现在只有函数 G,需要求 F,那么我们就需要寻找一种变换 ,使得 G 在经过这种变换后能够获得 F,这个过程即为反演, ...
- 【bzoj4176】Lucas的数论 莫比乌斯反演+杜教筛
题目描述 去年的Lucas非常喜欢数论题,但是一年以后的Lucas却不那么喜欢了. 在整理以前的试题时,发现了这样一道题目"求Sigma(f(i)),其中1<=i<=N" ...
- 数学--数论--莫比乌斯反演
一.莫比乌斯反演涉及知识 1.莫比乌斯函数 2.莫比乌斯的线性筛法 3.狄利克雷卷积 4.莫比乌斯反演详解 5.整除法分块 6.杜教筛 二.μ 莫比乌斯函数定义 μ(n)={1n=1(−1)kn= P ...
- 模板 - 数学 - 数论 - 莫比乌斯反演 - 2
新东西: 求$\sum\limits_{i=1}^{n}|\mu(i)|$ 根据莫比乌斯函数的性质,实际上就是求$\sum\limits_{i=1}^{\lfloor\sqrt{n}\rfloor}\ ...
- 【bzoj3601】一个人的数论 莫比乌斯反演+莫比乌斯函数性质+高斯消元
Description Sol 这题好难啊QAQ 反正不看题解我对自然数幂求和那里是一点思路都没有qwq 先推出一个可做一点的式子: \(f(n)=\sum_{k=1}^{n}[(n,k)=1]k^d ...
- 【笔记】莫比乌斯反演(1)
上一篇:莫比乌斯反演(0) 文章目录 莫比乌斯反演 关于反演 莫比乌斯函数 定义 性质 莫比乌斯反演公式 公式1 公式2 整除分块 引入 关于整除分块 基础推导 简单扩展 莫比乌斯反演的应用 例1:证 ...
- matlab狄利克雷函数,数论入门1——莫比乌斯函数,欧拉函数,狄利克雷卷积,线性筛,莫比乌斯反演,杜教筛...
数论入门1 一个菜鸡对数论的一点点理解... 莫比乌斯函数 定义函数$\mu(n)$为: 当n有平方因子时,$\mu(n)=0$. 当n没有平方因子时,$\mu(n)=(-1)^{\omega(n)} ...
- 《算法竞赛进阶指南》数论篇(3)-组合计数,Lucas定理,Catalan数列,容斥原理,莫比乌斯反演,概率与数学期望,博弈论之SG函数
文章目录 组合计数 例题:Counting swaps Lucas定理 Cnm≡Cnmodpmmodp∗Cn/pm/p(modp)C_n^m\equiv C_{n\ mod\ p}^{m\ mod\ ...
最新文章
- 潘石屹首次Python考试成绩 99 分,失分点:画完图后忘了隐藏画笔的箭头
- 自动化工具之二:win32gui
- LVS技术浅析-proc参数
- html偷拍代码,一段植入木马的html代码
- Java EE 6 VS Spring 3:Java EE杀死了Spring? 没门!
- [你必须知道的异步编程]——异步编程模型(APM)
- 微信公众号自定义菜单设置管理
- 命令行_Laravel-admin artisan 命令行脚本使用
- android 页面计时器,Android中页面添加计时器进行的刷新
- mysql盲注ascii_[翻译]关于通过对8bit的ascii做右位移提高mysql盲注效率
- 1602: [Usaco2008 Oct]牧场行走
- 同学,那些个证书没屌用!
- Unity3D AssetBundle相关
- c++采集声卡输出_windows上面捕获声卡数据
- 教育行业运维审计用什么堡垒机好?有什么作用?
- 国家区块链漏洞库 《区块链漏洞定级细则》发布
- 腾讯bugly升级注意事项
- 亲身经历,大龄程序员找工作,为什么这么难...
- 精准填报志愿、一分不浪费……靠谱吗?
- mysql查询当前用户中所有的表空间_oracle 查看用户所在的表空间