POJ 3090 Visible Lattice Points

题目链接:http://poj.org/problem?id=3090


题意:有一个n+1 * n+1 的点阵,从(0,0)号点向各个点连线,问一共可以看到多少个点( 比如(2,4)会被(1,2)挡到 )。


思路:考虑(x,y)是否能被看到,仅当x与y互质,否则就会有( x/gcd(x,y) , y/gcd(x,y) )去挡到这个点。所以我们直接把2到n的欧拉函数加一遍即可,对于一个可视点(i,j),i>j,必定存在(j,i)所以答案要乘二.还有特殊的三个点额外计算。


欧拉函数:http://blog.csdn.net/chy20142109/article/details/51263531


#include <cstdio>
#include <cmath>
#include <cstring>
#include <string>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <stack>
#include <map>
#include <set>
#include <vector>
#include <sstream>
#include <queue>
#include <utility>
using namespace std;

#define rep(i,j,k) for (int i=j;i<=k;i++)
#define Rrep(i,j,k) for (int i=j;i>=k;i--)

#define Clean(x,y) memset(x,y,sizeof(x))
#define LL long long
#define ULL unsigned long long
#define inf 0x7fffffff
#define mod %100000007

#define mp make_pair
#define fi first
#define se second

#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define mid int m=(l+r)>>1
const int maxn = 1009;
int mindiv[maxn],phi[maxn];


void getphi(int n)
{
    rep(i,1,n) mindiv[i] = i;
    for(int i = 2; i*i <= n; i++)
        if ( mindiv[i] == i )
        {
            for(int j = i*i; j <= n; j+=i) mindiv[j] = i;
        }
    phi[1] = 1;
    rep(i,2,n)
    {
        phi[i] = phi[i/mindiv[i]];
        if ( ( i / mindiv[i] ) % mindiv[i] == 0 ) phi[i] *= mindiv[i];
        else phi[i] *= ( mindiv[i] - 1 );
    }
}

int cal(int n)
{
    int ans = 3;
    rep(i,2,n) ans += 2*phi[i];
    return ans;
}

int main()
{
    getphi(1000);
    int T,n;
    cin>>T;
    rep(kase,1,T)
    {
        scanf("%d",&n);
        int ans = cal(n);
        printf("%d %d %d\n",kase,n,ans);
    }

    return 0;
}



  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值