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

UVa 11064 - Number Theory

时间:2014-10-28 09:25:08      阅读:222      评论:0      收藏:0      [点我收藏+]

标签:blog   io   os   ar   for   sp   数据   on   2014   

题目:求给顶一个数n,的所有的1 ≤ m ≤ n的m,使得gcd(m,n)≠ 1 且 gcd(m,n)≠ m。

分析:数论,素数筛法,欧拉函数。 

            设pi为n的第i个素数因,k1为第i个素数因子的个数,则有:

            1 ≤ m ≤ n,gcd(m,n)= 1 的m的个数为欧拉函数;

            欧拉函数:φ(n)= n *(1 - 1/p1)*(1 - 1/p2)*(1 - 1/p3)*…*(1 - 1/pt);

            1 ≤ m ≤ n,gcd(m,n)= m 的m的个数为n的所有因数的个数;

            因数个数:f(n)= (k1+1)*(k2+1)*...*(kt+1);

            这里利用筛法打表计算出150000内的素数,因为数据范围是20000000000内的,

            所以,不能被前150000内的素数整除的数,也一定是素数,并且每个数n中最多有一个;

            计算输出 n -  φ(n)- f(n)+ 1 即可{ gcd(1,n)计算了两次 }。

说明:UVa10299。

#include <iostream>  
#include <cstdlib>  
#include <cmath>
  
using namespace std;  
  
int fac[32],num[32];  
int prim[150000];  
int used[150000];  
  
int main()  
{  
    for (int i = 0 ; i < 150000 ; ++ i)  
        used[i] = 0;  
    int save = 0;  
    for (int i = 2 ; i < 150000 ; ++ i)  
        if (!used[i]) {  
            prim[save ++] = i;  
            for (int j = 2*i ; j < 150000 ; j += i)  
                used[j] = 1;  
        }  
      
    int n;  
    while (cin >> n && n) {  
        int count = 0,base = 0,m = n;  
        while (n > 1 && base < save) {  
            if (n%prim[base] == 0) {  
                fac[count] = prim[base];
                num[count] = 0;
                while (n%prim[base] == 0) {
                    n /= prim[base];  
                    num[count] ++;
				}
				count ++;
            }  
            base ++;
        }
        if (n > 1) {
			fac[count] = n;
			num[count] = 1;
			count ++;
		}
          
        long long ans = m;
        for (int i = 0 ; i < count ; ++ i)  
            ans = ans/fac[i]*(fac[i]-1);
        
        //gcd(m,n) = m
        int s = 1;
        for (int i = 0 ; i < count ; ++ i)
        	s = s*(num[i]+1);
        
        cout << m-ans-s+1 << endl;  
    }  
    return 0;  
}

UVa 11064 - Number Theory

标签:blog   io   os   ar   for   sp   数据   on   2014   

原文地址:http://blog.csdn.net/mobius_strip/article/details/40535315

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