2018.09.01【POJ3090】Visible Lattice Points (欧拉函数) (线性筛)(前缀和)

传送门


解析:

简单的欧拉函数套上一个线性筛。。。

考虑这个方阵的下三角(反正上面也是对称过去的,乘以2就行了),
i i 列除去x轴上的和对角线上的点,就还剩 i i 个点,而如果有gcd(i,j)!=1,则该点在 (i/gcd(i,j),j/gcd(i,j)) ( i / g c d ( i , j ) , j / g c d ( i , j ) ) 时被收集,所以第 i i 列能够收集到的点的个数就是i以下与 i i <script type="math/tex" id="MathJax-Element-8">i</script>互质的数的个数。

多么明显的欧拉函数。。。
最后加上(1,1)这个在对角线上的东西。。。


代码:

#include<iostream>
#include<cctype>
#include<cstdio>
using namespace std;
#define re register
#define ll long long
#define gc getchar
#define pc putchar
#define cs const
#define st static

inline
ll getint(){
    re ll num=0;
    re char c=gc();
    for(;!isdigit(c);c=gc());
    for(;isdigit(c);c=gc())num=(num<<1)+(num<<3)+(c^48);
    return num;
}

inline
void outint(ll a){
    st char ch[23];
    if(a==0)pc('0');
    while(a)ch[++ch[0]]=(a-a/10*10)^48,a/=10;
    while(ch[0])pc(ch[ch[0]--]);
}

bool mark[1005];
int prime[1005],pcnt;
int phi[1005];
ll sum[1005];

inline
void linear_sieves(int len=1000){
    mark[1]=true;
    phi[1]=1;
    for(int re i=2;i<=len;++i){
        if(!mark[i])prime[++pcnt]=i,phi[i]=i-1;
        for(int re j=1;j<=pcnt&&i*prime[j]<=len;++j){
            mark[i*prime[j]]=true;
            if(i%prime[j]==0){
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
            phi[i*prime[j]]=phi[i]*phi[prime[j]];
        }
    }
}

inline
void init(){
    linear_sieves();
    sum[0]=0;
    for(int re i=1;i<=1000;++i){
        sum[i]=sum[i-1]+phi[i];
    }
}

int T;
int main(){
    init();
    T=getint();
    for(int re kase=1;kase<=T;++kase){
        int n=getint();
        outint(kase),pc(' '),outint(n),pc(' '),outint(sum[n]*2+1),pc('\n');
    }
    return 0;
}
  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值