码迷,mamicode.com
首页 > 其他好文 > 详细

HDU 2841

时间:2014-10-02 18:06:23      阅读:133      评论:0      收藏:0      [点我收藏+]

标签:blog   io   os   for   sp   div   c   log   amp   

明显,当(X,Y)=1时,是可以看见的。

这题,记得POJ 上也有类似的一题。。。

不过比较奇怪的是,我以为会超时,因为范围达到了100000,但竟然直接枚举没超时。。。。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#define LL __int64
#define N 100002
using namespace std;

bool isprime[N];
LL prime[N],np;
LL fac[100],fp;


void initial(){
	memset(isprime,true,sizeof(isprime));
	np=0;
	for(LL i=2;i<N;i++){
		if(isprime[i]){
			prime[np++]=i;
			for(LL j=i*i;j<N;j+=i)
			isprime[j]=false;
		}
	}
}

void dfs(LL i,LL num,LL p,LL &ans,LL s,LL &n){
	if(i>=s){
		if(num==0)
		ans=n;
		else if(num&1)
		ans-=(n/p);
		else ans+=(n/p);
		return ;
	}
	dfs(i+1,num,p,ans,s,n);
	dfs(i+1,num+1,p*fac[i],ans,s,n);
}

LL work(LL m,LL n){
	fp=0;
 	for(LL i=0;i<np&&prime[i]*prime[i]<=m;i++){
	 	if(m%prime[i]==0){
	 		while(m%prime[i]==0){
	 			m/=prime[i];
	 		}
	 		fac[fp++]=prime[i];
	 	}
 	}
 	if(m>1) fac[fp++]=m;
 	
 	LL ans;
 	dfs(0,0,1,ans,fp,n);
 	return ans;
 	
}

int main(){
	int T;
	initial();
	scanf("%d",&T);
	while(T--){
		LL m,n;
		scanf("%I64d%I64d",&m,&n);
		LL ans=n;
		for(LL i=2;i<=m;i++)
			ans+=work(i,n);
		printf("%I64d\n",ans);
	}
	return 0;
}

  

HDU 2841

标签:blog   io   os   for   sp   div   c   log   amp   

原文地址:http://www.cnblogs.com/jie-dcai/p/4004096.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!