POJ3090.Visible Lattice Points

最新推荐文章于 2021-04-12 17:56:37 发布

Fiveneves

最新推荐文章于 2021-04-12 17:56:37 发布

阅读量106

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/Fiveneves/article/details/104250968

版权

数论专栏收录该内容

13 篇文章 0 订阅

订阅专栏

3090.Visible Lattice Points

题目链接-Visible Lattice Points
在这里插入图片描述
Sample Input

Sample Output

1 2 5
2 4 13
3 5 21
4 231 32549

题目大意
给你一个n*n的网格，任意一点和原点连线，可以形成一条直线，前面的点可以挡住后面的点，问你能看到的点到底有多少个

解题思路

同一条直线上的点，与原点连线的斜率是一样的，所有其他点与原点的连线斜率都是y/x，如果y/x能约分，说明前面肯定有一个点挡住这个点
根据对称性，可将整个正方形分成上三角和下三角两部分，两边可以看到的点的数目肯定一样多
只需要求1到n中互质数对的个数即可求出答案，枚举1到i-1的数字寻找互质数对，累加斜率数量即可
记忆化对代码进行优化

附上代码O(n²)

#include<iostream>
#define int long long
using namespace std;
const int INF=0x3f3f3f3f;
const int N=1e5+5;
const int M=1e9+7;
typedef long long ll;
typedef pair<int,int> PII;
int gcd(int a,int b){//辗转相除求gcd
	int t;
    if(a<b)
	swap(a,b); 
    int ans=1;
    while(b!=0){
    	t=a%b;
    	a=b;
    	b=t;
	}
	return a;
}
int s[1010];
void f(){
	s[1]=2;//斜率为1及斜率不存在时
	for(int i=2;i<=1000;i++){
		s[i]=s[i-1];//加上前面的斜率
		for(int j=1;j<i;j++)
			if(gcd(i,j)==1)
				s[i]+=2;//累加，由于对称性所以是+2
	}
}
signed main(){
	ios::sync_with_stdio(0);
	cin.tie(0);cout.tie(0);
	
	int c;
	cin>>c;
	f();
	for(int i=1;i<=c;i++){
		int n,ans=0;
		cin>>n;
		cout<<i<<" "<<n<<" "<<s[n]+1<<endl;
		//加上中间的斜率(即正方形对角线)
	} 
	return 0;
}

我们写的f()函数其实是找与 i 互质的数字的数目的2倍，这就把问题转化为求小于 n 的与 n 互质的数的个数，即欧拉函数，欧拉函数可以通过质因数分解来快速计算出来，为了进一步提高效率可以再加上线性筛

附上代码O(n)

#include<cstdio>
#define int long long
using namespace std;
const int INF=0x3f3f3f3f;
const int N=1e5+5;
const int M=1e9+7;
typedef long long ll;
typedef pair<int,int> PII;
int phi[1010],prime[1010],tot=0;
int  mark[1010];
void f(){
    memset(mark,0,sizeof(mark));
    phi[1]=1;
    for(int i=2;i<=maxn;i++)
    {
        if(!mark[i]){
            prime[++tot]=i;
            phi[i]=i-1;
        }
        for(int j=1;j<=tot;j++){
            if(i*prime[j]>maxn) break;
            mark[i*prime[j]]=1;
            if(i%prime[j]==0){
                phi[i*prime[j]]=prime[j]*phi[i];
                break;
            }
            else{
                phi[i*prime[j]]=phi[i]*(prime[j]-1);
            }
        }
    }
}
signed main(){
    int t;
    f();
    scanf("%d",&t);
    for(int i=1;i<=t;i++){
        int n,ans=0;
        scanf("%d",&n);
        for(int j=1;j<=n;j++){
            ans+=phi[j];
        }
        printf("%d %d %d \n",time,n,ans*2+1);
    }
    return 0;
}