POJ3090,Visible Lattice Points(GCD)

最新推荐文章于 2019-08-01 20:49:30 发布

香蕉君

最新推荐文章于 2019-08-01 20:49:30 发布

阅读量217

点赞数

分类专栏： # 组合数学

本文链接：https://blog.csdn.net/qq_37957829/article/details/82151128

版权

组合数学专栏收录该内容

5 篇文章 0 订阅

订阅专栏

POJ3090,Visible Lattice Points

题意：

在一个直角坐标系的第一象限上有类似下面的图。

有的点可以从原点直接看到，有的店恰好被其他的点挡住。给一个n，代表N*N的范围。求，范围内有多少个这样的点可以直接被看到。

解：

被遮住的点有一个共同点，假设点为 $(x, y)$ ，那么 $G C D (x, y) = d, d > 1$ ，因为它一定会被 $(x / d, y / d)$ 遮住，所以我们找到2到n之间全部的互质的数有多少，再加上三个点(1,0)(0,1)(1,1)就是答案。答案就是从1到n的欧拉函数的前缀和。公式为：

$2\sum_{i=2}^{n}\phi(i)+3$
巧妙。

####代码：

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
using namespace std;
const int MAXN = 1e5+10;
bool flag[MAXN];///标记数组
long long phi[MAXN];///欧拉函数值，i的欧拉函数值=phi[i]
long long ans[MAXN];
int p[MAXN];///素因子的值
int cnt = 0;
void Get_phi()///筛法求欧拉函数
{
    cnt = 0;
    memset(flag, true, sizeof(flag));
    phi[1] = 1;
    for(int i=2; i<MAXN; i++)///线性筛法
    {
        if(flag[i])///素数
        {
            p[cnt++] = i;
            phi[i] = i-1;///素数的欧拉函数值是素数 - 1
        }
        for(int j=0; j<cnt; j++)
        {
            if(i*p[j] > MAXN)
                break;
            flag[i*p[j]] = false;///素数的倍数，所以i*p[j]不是素数
            if(i%p[j] == 0)///性质：i mod p == 0, 那么 phi(i * p) == p * phi(i)
            {
                phi[i*p[j]] = p[j] * phi[i];
                break;
            }
            else
                phi[i*p[j]] = (p[j]-1) * phi[i];///i mod p != 0, 那么 phi(i * p) == phi(i) * (p-1)
        }
    }
}

int main()
{
    int n;
    Get_phi();
    ans[2]=phi[2]*2;
    for(int i=3;i<MAXN;i++)
        ans[i]=ans[i-1]+(phi[i]<<1);
    int cas=1,T;
    scanf("%d",&T);
    while(T--){
        scanf("%d",&n);
        printf("%d %d %lld\n",cas++,n,ans[n]+3);
    }
    return 0;
}