首页 > 其他 > 详细

HDU 1695

时间:2014-11-14 00:08:57      阅读:339      评论:0      收藏:0      [点我收藏+]

http://acm.hdu.edu.cn/showproblem.php?pid=1695

x是[1,b],y是[1,d],求GCD(x,y)=k的对数(x,y无序)

对x,y都除以k,则求GCD(x,y)=1

此时枚举x,问题转化为[1,d]区间内与x互素的数字个数,这个问题是hdu 4135

有一个特殊的地方是x,y无序,对于这点只要保证x始终小于y就可以了

特判k=0

bubuko.com,布布扣
#include <iostream>
#include <cstdio>
#include <cstring>
#include <vector>

using namespace std;

typedef __int64 ll;

vector <int> v;
int b,d,k;

int main(){
    int T;
    scanf("%d",&T);
    for(int cas=1;cas<=T;cas++){
        scanf("%*d%d%*d%d%d",&b,&d,&k);
        if(!k){
            printf("Case %d: 0\n",cas);
            continue;
        }
        b/=k;d/=k;
        if(b>d)swap(b,d);
        int i;
        ll ans=0; 
        for(i=1;i<=b;i++){
            int temp=i;
            v.clear();
            for(int j=2;j*j<=i;j++){
                if(i%j==0){
                    v.push_back(j);
                    while(i%j==0)i/=j;
                }
            }
            if(i>1)v.push_back(i);
            int m=v.size();
            i=temp;
            ll res=0; 
            for(int j=1;j<(1<<m);j++){
                int cnt=0;
                ll val=1;
                for(int k=0;k<m;k++){
                    if(j&(1<<k)){
                        cnt++;
                        val*=v[k];
                    }
                }
                if(cnt&1)res+=d/val-(i-1)/val;
                else res-=d/val-(i-1)/val;
            }
            ans+=(d-i+1)-res;
        }
        printf("Case %d: %I64d\n",cas,ans);
    }
    return 0;
}
View Code

 

HDU 1695

原文:http://www.cnblogs.com/xiaohongmao/p/4096144.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!