由题引出欧拉函数 洛谷P2158

https://www.luogu.com.cn/problem/P2158

题面：
$我们定义C君坐标(0,0)，我们猜测能否看到后面的人可能会与斜率有关。$
$如果点(i,j)可以被看到，那么连接它(0,0)的直线之间没有点，$$也就是这条线之间没有点的斜率等于它，$$也就是\frac{j}{i}不能约分。也就是\gcd (i,j)=1。$
$那么我们只需要找到每个j，有多少i满足\gcd (i,j)=1。$
$而这刚好就是欧拉函数。$
$Euler函数：\phi (n)表示在1\sim n范围内有多少个数x满足\gcd (n,x)=1。$
$显然,如果n为质数,\phi (n)=n-1。$
$如果n不是质数\dots \dots 不想写了$—_—|||$以后再来补$
$这里给出\log 求欧拉函数的方法：$

int euler(int n) {
    int ans = n;
    for (int i = 2; i * i <= n; ++i) {
        if (n % i == 0) {
            ans = ans / i * (i - 1);
            while (n % i == 0)n /= i;
        }
    }
    if (n > 1)ans = ans / n * (n - 1);
    return ans;
}

$那么这道题的代码就是：$

#include<bits/stdc++.h>

using namespace std;
//    clock_t start, end;
//    start = clock();
//    end = clock();
//    cout << (double) (end - start) / CLOCKS_PER_SEC << endl;
//ios::sync_with_stdio(false);
#define  int long long
#define double long double
#define rep(i, x, y) for(int i=(x);i<=(y);++i)
#define dep(i, x, y) for(int i=(x);i>=(y);--i)
#define gcd(a, b) __gcd(a,b)
const long long mod = 1e9 + 7;
const int maxn = 1e5 + 10;

int lowbit(int x) { return x & -x; }

bool ispow(int n) { return (n & (n - 1)) == 0; }//O(1) 判断是否是 2^k（2的k次方）
//12200160415121876738


int euler(int n) {
    int ans = n;
    for (int i = 2; i * i <= n; ++i) {
        if (n % i == 0) {
            ans = ans / i * (i - 1);
            while (n % i == 0)n /= i;
        }
    }
    if (n > 1)ans = ans / n * (n - 1);
    return ans;
}
int f[maxn];
signed main(){
    //freopen("01.in","r",stdin);
    int n;
    cin>>n;
    if(n == 1){
        cout << 0 << endl;
        return 0;
    }
    int ans = 3;
    for(int i =2;i<n;i++){
        ans += euler(i) * 2;
    }
    cout<<ans<<endl;
   return 0;

}

$线性筛求\phi (n):$

#include<bits/stdc++.h>

using namespace std;
//    clock_t start, end;
//    start = clock();
//    end = clock();
//    cout << (double) (end - start) / CLOCKS_PER_SEC << endl;
//ios::sync_with_stdio(false);
#define  int long long
#define double long double
#define rep(i, x, y) for(int i=(x);i<=(y);++i)
#define dep(i, x, y) for(int i=(x);i>=(y);--i)
#define gcd(a, b) __gcd(a,b)
const long long mod = 1e9 + 7;
const int maxn = 1e5 + 10;

int lowbit(int x) { return x & -x; }

bool ispow(int n) { return (n & (n - 1)) == 0; }//O(1) 判断是否是 2^k（2的k次方）
//12200160415121876738

int vis[maxn],prim[maxn],ff[maxn],cnt,num[maxn];
void f(/*线性筛*/) {
    for (int i = 2; i <= 50000; i++) {
        if (!vis[i])
            prim[cnt++] = i, num[i] = 1, ff[i] = i - 1;
        for (int j = 0; j < cnt && i * prim[j] <= 50000; j++) {
            vis[prim[j] * i] = 1;
            num[prim[j] * i] = num[i] + 1;
            if (i % prim[j] == 0) {
                ff[i * prim[j]] = prim[j] * ff[i];
                break;
            }
            ff[prim[j] * i] = ff[i] * ff[prim[j]];
        }
    }
}
signed main(){
    //freopen("01.in","r",stdin);
    int n;
    cin>>n;
    f();
    if(n == 1){
        cout << 0 << endl;
        return 0;
    }
    int ans = 3;
    for(int i =2;i<n;i++){
        ans += ff[i]*2;
    }
    cout<<ans<<endl;
   return 0;

}