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

UVA 10820 Send a Table 数论 欧拉函数

时间:2017-08-22 17:36:30      阅读:133      评论:0      收藏:0      [点我收藏+]

标签:描述   scanf   直接   net   typedef   logs   end   ide   http   

  题目链接: https://vjudge.net/problem/UVA-10820

  题目描述: 给你一个N, N <= 50000, 让你寻找N之内互素数的个数

  解题思路: 欧拉函数, 由于位置颠倒是两个解, 在小于N的范围内只有(1, 1)x, y相等, 其他的都是不等的, 所以我们只需要算phi(2) + phi(3) + ...... phi(n), 然后 * 2 + 1即可

  代码: 

技术分享
#include <iostream>
#include <cstdio>
#include <string>
#include <vector>
#include <cstring>
#include <iterator>
#include <cmath>
#include <algorithm>
#include <stack>
#include <deque>
#include <map>
#define lson l, m, rt<<1
#define rson m+1, r, rt<<1|1
#define mem0(a) memset(a,0,sizeof(a))
#define meminf(a) memset(a,-0x3f,sizeof(a))
#define fi(n) for(i=0;i<n;i++)
#define fj(m) for(j=0;j<m;j++)
#define sca(x) scanf("%d",&x)
#define scalld(x) scanf("%I64d",&x)
#define print(x) printf("%d\n", x)
#define printlld(x) printf("%I64d\n",x)
#define de printf("=======\n")
#define yes printf("YES\n")
#define no printf("NO\n")
typedef long long ll;
using namespace std;

const int maxn = 50000+100;
int phi[maxn];

void phi_table( int n ) {
    for( int i = 2; i <= n; i++ ) {
        phi[i] = 0;
    }
    phi[1] = 1;
    for( int i = 2; i <= n; i++ ) {
        if( !phi[i] ) {
            for( int j = i; j <= n; j += i ) {
                if( !phi[j] ) phi[j] = j;
                phi[j] = phi[j] / i * (i-1);
            }
        }
    }
}

int main() {
    int n;
    mem0(phi);
    phi_table(50000);
    while( sca(n) == 1 && n ) {
        if( n == 1 ) {
            printf( "1\n" );
            continue;
        }
        else {
            int res = 0;
            for( int i = 2; i <= n; i++ ) {
                res += phi[i];
            }
            printf( "%d\n", 2 * res + 1 );
        }
    }
    return 0;
}
View Code

  思考: 欧拉函数, 复杂度为O(nlognlogn), 本题就是直接打表然后计算和, 要注意这种思想, 倍数可以直接得到就转化成互素这种思想, 好好学数学

UVA 10820 Send a Table 数论 欧拉函数

标签:描述   scanf   直接   net   typedef   logs   end   ide   http   

原文地址:http://www.cnblogs.com/FriskyPuppy/p/7412539.html

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