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

【POJ 3090】 Visible Lattice Points

时间:2018-07-09 14:17:49      阅读:165      评论:0      收藏:0      [点我收藏+]

标签:com   通过   坐标   code   limits   clu   前缀和   mes   return   

【题目链接】

            http://poj.org/problem?id=3090

【算法】

              通过观察发现,在这个平面直角坐标系中,除了(1,1),(1,0)和(0,1),所有可见点的横纵坐标互质

              那么,问题就转化为了求 2 * (phi(1) + phi(2) + ... + phi(n)) + 3

              预处理phi的前缀和,O(1)回答询问,即可

【代码】

             

#include <algorithm>  
#include <bitset>  
#include <cctype>  
#include <cerrno>  
#include <clocale>  
#include <cmath>  
#include <complex>  
#include <cstdio>  
#include <cstdlib>  
#include <cstring>  
#include <ctime>  
#include <deque>  
#include <exception>  
#include <fstream>  
#include <functional>  
#include <limits>  
#include <list>  
#include <map>  
#include <iomanip>  
#include <ios>  
#include <iosfwd>  
#include <iostream>  
#include <istream>  
#include <ostream>  
#include <queue>  
#include <set>  
#include <sstream>  
#include <stdexcept>  
#include <streambuf>  
#include <string>  
#include <utility>  
#include <vector>  
#include <cwchar>  
#include <cwctype>  
#include <stack>  
#include <limits.h>
using namespace std;
#define MAXN 1010

int i,n,TC,T;
int sum[MAXN];

inline void init()
{
        int i,j,tmp,cnt = 0;
        static int f[MAXN],prime[MAXN],phi[MAXN];
        for (i = 2; i < MAXN; i++)
        {
                if (!f[i]) 
                {
                        f[i] = prime[++cnt] = i;
                        phi[i] = i - 1;
                }
                for (j = 1; j <= cnt; j++)
                {
                        tmp = i * prime[j];
                        if (tmp >= MAXN) break;
                        f[tmp] = prime[j];
                        if (prime[j] == f[i])
                        {
                                phi[tmp] = phi[i] * prime[j];
                                break;        
                        }    else phi[tmp] = phi[i] * (prime[j] - 1);
                }        
        }        
        for (i = 1; i < MAXN; i++) sum[i] = sum[i-1] + phi[i];
}

int main() 
{
        
        init();
        scanf("%d",&T);
        while (T--)
        {
                scanf("%d",&n);
                printf("%d %d %d\n",++TC,n,2*sum[n]+3);
        }
        return 0;
    
}

 

【POJ 3090】 Visible Lattice Points

标签:com   通过   坐标   code   limits   clu   前缀和   mes   return   

原文地址:https://www.cnblogs.com/evenbao/p/9283344.html

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