hdu4135Co-prime 容斥原理水题
2016-03-16 11:30
309 查看
//问一个区间[a,b]与n互素的数的个数
//利用容斥原理可知
//在[a,b] 区间内对n的素数因子
//ans = 被一个数整除的数的个数 - 被两个数的最小公倍数整除的数的个数 + 被三个数的。。。
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std ;
const int maxn = 100010 ;
typedef __int64 ll ;
ll p[maxn] ;int len ;
void get_prime(ll n)
{
len = 0 ;
for(ll i = 2;i*i <= n;i++)
{
if(n%i == 0)p[++len] = i ;
while(n%i == 0)n/=i ;
}
if(n>1)p[++len] = n;
}
ll dfs(int pos , ll n)
{
ll ans = 0 ;
for(int i = pos ;i <= len ;i++)
ans += n/p[i] - dfs(i+1 , n/p[i]) ;
return ans ;
}
int main()
{
ll a , b ,n ;
int T ;
int cas = 0 ;
scanf("%d" ,&T) ;
while(T--)
{
scanf("%I64d%I64d%I64d" , &a , &b , &n);
get_prime(n) ;
ll ans = (b - dfs(1 , b)) - (a - 1 - dfs(1 , a-1)) ;
printf("Case #%d: " ,++cas) ;
printf("%I64d\n" , ans) ;
}
return 0 ;
}
//利用容斥原理可知
//在[a,b] 区间内对n的素数因子
//ans = 被一个数整除的数的个数 - 被两个数的最小公倍数整除的数的个数 + 被三个数的。。。
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std ;
const int maxn = 100010 ;
typedef __int64 ll ;
ll p[maxn] ;int len ;
void get_prime(ll n)
{
len = 0 ;
for(ll i = 2;i*i <= n;i++)
{
if(n%i == 0)p[++len] = i ;
while(n%i == 0)n/=i ;
}
if(n>1)p[++len] = n;
}
ll dfs(int pos , ll n)
{
ll ans = 0 ;
for(int i = pos ;i <= len ;i++)
ans += n/p[i] - dfs(i+1 , n/p[i]) ;
return ans ;
}
int main()
{
ll a , b ,n ;
int T ;
int cas = 0 ;
scanf("%d" ,&T) ;
while(T--)
{
scanf("%I64d%I64d%I64d" , &a , &b , &n);
get_prime(n) ;
ll ans = (b - dfs(1 , b)) - (a - 1 - dfs(1 , a-1)) ;
printf("Case #%d: " ,++cas) ;
printf("%I64d\n" , ans) ;
}
return 0 ;
}
相关文章推荐
- C# 微信支付证书使用
- 如何快速接手一个系统?
- maven工程中dubbo与spring整合
- Axis2创建webservice实例
- C#状态机
- Java 伪随机数
- redis.conf 配置文件详解
- IOS sqlite数据库文件查看(工具)
- Auto Layout 使用心得(六)—— 制造炫酷的下拉刷新动画
- NameValueCollection类
- ESP8266模块使用完整教程
- The method setLatestEventInfo(Context, CharSequence, CharSequence, PendingIntent) from the type
- 查看Oracle 是32位还是64位的方法
- 创建OpenStack外部网络并分配浮动IP
- .net如何判断网页是否由搜索引擎蜘蛛访问?
- iOS开发之runtime的运用-获取当前网络状态
- Html5 android 使用总结
- 转:java.lang.IllegalStateException异常产生的原因及解决办法
- 今天交作业有点慢。
- 移植代码到linux内核中大概流程