hdu1695(莫比乌斯)或欧拉函数+容斥
题意:求1-b和1-d之内各选一个数组成数对。问最大公约数为k的数对有多少个,数对是有序的。(b,d,k<=100000)
解法1: 这个能够简化成1-b/k 和1-d/k 的互质有序数对的个数。
如果b=b/k。d=d/k,b<=d.欧拉函数能够算出1-b与1-b之内的互质对数。然后在b+1到d的数i,求每一个i在1-b之间有多少互质的数。
解法是容斥,getans函数參数的意义:1-tool中含有rem位置之后的i的质因子的数的个数。
在
for(int j=rem;j<=factor[i][0];j++)
ans+=tool/factor[i][j]-getnum(i,tool/factor[i][j],j+1);
这个循环中。ans加的等号后每项表示当前最大的质因子是factor[i][j]的数量,目的是去重。
解法2:莫比乌斯,莫比乌斯数组确实非常实用。事实上也非常easy,mou的位置的含义是,首先假设i有个质因子出现2次或以上。则mou值为0,否则1与-1跟i的质因子奇偶性决定。目的也是容斥。
解法1代码:
/******************************************************
* @author:xiefubao
*******************************************************/
#pragma comment(linker, "/STACK:102400000,102400000")
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <queue>
#include <vector>
#include <algorithm>
#include <cmath>
#include <map>
#include <set>
#include <stack>
#include <string.h>
//freopen ("in.txt" , "r" , stdin);
using namespace std; #define eps 1e-8
#define zero(_) (abs(_)<=eps)
const double pi=acos(-1.0);
typedef long long LL;
const int Max=100010;
const LL INF=0x3FFFFFFF;
int a, b, c, d, k;
int factor[Max][20];
bool rem[Max];
LL oular[Max];
void init()
{
for(int i=1; i<Max; i++)
oular[i]=i;
for(LL i=2; i<Max; i++)
{
if(!rem[i])
{
factor[i][++factor[i][0]]=i;
oular[i]=i-1;
for(LL j=i*2; j<Max; j+=i)
{
factor[j][++factor[j][0]]=i;
rem[j]=1;
oular[j]=oular[j]*(i-1)/i;
}
}
}
for(int i=1; i<Max; i++)
oular[i]=oular[i]+oular[i-1];
}
LL getnum(int i,int tool,int rem)
{
int ans=0;
for(int j=rem;j<=factor[i][0];j++)
ans+=tool/factor[i][j]-getnum(i,tool/factor[i][j],j+1);
return ans;
}
int main()
{
int t;
cin>>t;
init();int kk=1;
while(t--)
{
cin>>a>>b>>c>>d>>k;
printf("Case %d: ",kk++);
if(k==0)
{
cout<<"0\n";
continue;
}
b/=k;
d/=k;
if(b>d)swap(b,d);
LL ans=oular[b];
for(int i=b+1;i<=d;i++)
{
ans+=b-getnum(i,b,1);
}
cout<<ans<<endl;
}
return 0;
}
解法2代码:
/******************************************************
* @author:xiefubao
*******************************************************/
#pragma comment(linker, "/STACK:102400000,102400000")
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <queue>
#include <vector>
#include <algorithm>
#include <cmath>
#include <map>
#include <set>
#include <stack>
#include <string.h>
//freopen ("in.txt" , "r" , stdin);
using namespace std; #define eps 1e-8
#define zero(_) (abs(_)<=eps)
const double pi=acos(-1.0);
typedef long long LL;
const int Max=100010;
const LL INF=0x3FFFFFFF;
int a, b, c, d, k;
bool rem[Max];
int mou[Max];
void init()
{
mou[1]=1;
for(LL i=2; i<Max; i++)
{
if(!rem[i])
{
mou[i]=i;
for(LL j=i*2; j<Max; j+=i)
{
rem[j]=1;
mou[j]=i;
}
}
}
for(int i=2; i<Max; i++)
{
if(i/mou[i]%mou[i]==0) mou[i]=0;
else mou[i]=-mou[i/mou[i]];
}
}
int main()
{
int t;
cin>>t;
init();
int kk=1;
while(t--)
{
cin>>a>>b>>c>>d>>k;
printf("Case %d: ",kk++);
if(k==0)
{
cout<<"0\n";
continue;
}
b/=k;
d/=k;
if(b > d)swap(b,d);
long long ans1 = 0;
for(int i = 1; i <= b; i++)
ans1 += (long long)mou[i]*(b/i)*(d/i);
long long ans2 = 0;
for(int i = 1; i <= b; i++)
ans2 += (long long)mou[i]*(b/i)*(b/i);
ans1 -= ans2/2;
cout<<ans1<<endl;
}
return 0;
}
hdu1695(莫比乌斯)或欧拉函数+容斥的更多相关文章
- HDU 1695 GCD 欧拉函数+容斥定理 || 莫比乌斯反演
GCD Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submiss ...
- hdu 6390 欧拉函数+容斥(莫比乌斯函数) GuGuFishtion
http://acm.hdu.edu.cn/showproblem.php?pid=6390 题意:求一个式子 题解:看题解,写代码 第一行就看不出来,后面的sigma公式也不会化简.mobius也不 ...
- hdu 1695 GCD(欧拉函数+容斥)
Problem Description Given 5 integers: a, b, c, d, k, you're to find x in a...b, y in c...d that GCD( ...
- HDU 1695 GCD 欧拉函数+容斥定理
输入a b c d k求有多少对x y 使得x在a-b区间 y在c-d区间 gcd(x, y) = k 此外a和c一定是1 由于gcd(x, y) == k 将b和d都除以k 题目转化为1到b/k 和 ...
- HDU1695-GCD(数论-欧拉函数-容斥)
GCD Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total Submis ...
- $BZOJ$2818 $gcd$ 莫比乌斯反演/欧拉函数
正解:莫比乌斯反演/欧拉函数 解题报告: 传送门$QwQ$ 一步非常显然的变形,原式=$\sum_{d=1,d\in prim}^{n}\sum_{i=1}^{n}\sum_{j=1}^{n}[gcd ...
- BZOJ2005:[NOI2010]能量采集(莫比乌斯反演,欧拉函数)
Description 栋栋有一块长方形的地,他在地上种了一种能量植物,这种植物可以采集太阳光的能量.在这些植物采集能量后,栋栋再使用一个能量汇集机器把这些植物采集到的能量汇集到一起. 栋栋的植物种得 ...
- Bzoj 2818: Gcd 莫比乌斯,分块,欧拉函数,线性筛
2818: Gcd Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 3241 Solved: 1437[Submit][Status][Discuss ...
- luogu2658 GCD(莫比乌斯反演/欧拉函数)
link 给定整数N,求1<=x,y<=N且Gcd(x,y)为素数的数对(x,y)有多少对. 1<=N<=10^7 (1)莫比乌斯反演法 发现就是YY的GCD,左转YY的GCD ...
随机推荐
- 部署ASP.NET MVC项目
目标:了解部署过程,掌握部署中出现问题该如何处理. 部署网站往往是一件麻烦事,因为在安装部署的过程中,经常有许多步骤要运行,对于许多不太熟悉IIS/SQL的新手来说,部署网站编程一件非常困难且危险的事 ...
- larbin是一种开源的网络爬虫/网络蜘
larbin是一种开源的网络爬虫/网络蜘蛛,由法国的年轻人 Sébastien Ailleret独立开发.larbin目的是能够跟踪页面的url进行扩展的抓取,最后为搜索引擎提供广泛的数据来源.Lar ...
- MapReduce调度与执行原理之任务调度(续)
前言 :本文旨在理清在Hadoop中一个MapReduce作业(Job)在提交到框架后的整个生命周期过程,权作总结和日后参考,如有问题,请不吝赐教.本文不涉及Hadoop的架构设计,如有兴趣请参考相关 ...
- IOS的处理touch事件处理(按照手指的移动移动一个圆,开发环境用的ios7,storyboard)
先看下页面的效果图: 首先定义这个ball它有两个属性和两个方法: @property(nonatomic) CGPoint location; @property(nonatomic) CGFloa ...
- sharePoint常用命令
New-SPStateServiceDatabase -Name "StateServiceDatabase" | New-SPStateServiceApplication -N ...
- 幻世(OurDream)2D图形引擎使用教程11——播放媒体文件(1)
声明:本教程版权归Lizcst Software Lab所有,欢迎转载,但是转载必须保留本段声明文字,并注明文章来源:http://blog.csdn.net/kflizcst 谢谢合作! 播放媒体是 ...
- javascript每日一练(十)——运动二:缓冲运动
一.缓冲运动 实现原理:(目标距离-当前距离) / 基数 = 速度(运动距离越大速度越小,运动距离和速度成反比) (500 - oDiv.offsetLeft) / 7 = iSpeed; 需要注意: ...
- JVM调优总结(十)-调优方法
JVM调优工具 Jconsole,jProfile,VisualVM Jconsole : jdk自带,功能简单,但是可以在系统有一定负荷的情况下使用.对垃圾回收算法有很详细的跟踪.详细说明参考这里 ...
- POJ 3986 Math teacher's homework
题目 给出\(n,m_1,m_2,...,m_n\),求\(x_1 xor x_2 xor ... xor x_n=k (0 \leq x_i \leq m_i)\)的解的数量.二进制位数小于\(32 ...
- linux i2c驱动架构-dm368 i2c驱动分析
linux i2c驱动架构-dm368 i2c驱动分析 在阅读本文最好先熟悉一种i2c设备的驱动程序,并且浏览一下i2c-core.c以及芯片提供商的提供的i2c总线驱动(i2c-davinc ...