[BZOJ2005]-线性筛欧拉函数

说在前面

这破题me在BZoj上wrong answer5次，之前真的是没发现什么错，对拍也没有问题。
然后发现是爆int了······
（搬电脑摔！）

---

题目

题目传送门：BZOJ2005

题目大意

虽然是中文题，还是简述一下
就是要让我们求：$ans = \sum_{i=1}^n\sum_{j=1}^n(2*\gcd(i,j)-1)$

---

题目分析

我们需要知道这些

1. 一个点阵中，以坐标$(0,0)$和$(x,y)$ 为两个端点的线段经过的整数点个数是 $\gcd(x,y)-1$ 。

   证明：设$x/ \gcd(x,y)$为k，$y/ \gcd(x,y)$为q，那么整数点即为$(k,q), (2k,2q), (3k,3q)，\cdots ，(x,y)$。又因为点(x,y)是端点，因此减1。点数即为 $\gcd(x,y)-1$ 。

2. $n = \sum_{d|n}\varphi(d)I(\frac{n}{d})$。
   简单理解一下吧，I函数恒等于1（定义）。那么上式就是$n = \sum_{d|n}\varphi(d)$。对于$1$到$N$的每一个数字$i$，都与$N$有唯一一个最大公因数$d$。那么所有$d$的$\varphi(d)$就一定不重不漏的包含了$1$到$N$的每一个数字

于是由题意我们可以立刻写出$ans = \sum_{i=1}^n\sum_{j=1}^n(2*\gcd(i,j)-1)$
但是上面的公式，i和j都要枚举到1e5，时间过不去。
于是开始推导
    $\sum_{i=1}^n\sum_{j=1}^n(2*\gcd(i,j)-1)$
=$2*\sum_{i=1}^n\sum_{j=1}^n\gcd(i,j)-n*m$

对于$\sum_{i=1}^n\sum_{j=1}^n\gcd(i,j)$

$=\sum_{i=1}^n\sum_{j=1}^n\sum_{d|\gcd(i,j)}\varphi(d)$

$=\sum_{i=1}^n\sum_{j=1}^n\sum_{d=1}^n[d|i][d|j]*\varphi(d)$     (其中[]是bool，成立为1不成立为0)

$=\varphi(d)*\sum_{d=1}^n\sum_{i=1}^n[d|i]\sum_{j=1}^n[d|j]$

$=\varphi(d)*\sum_{d=1}^n* \lfloor\frac{n}{d}\rfloor*\lfloor\frac{m}{d}\rfloor$

于是就化为了一个$O(n)$内就可以出答案的式子。

---

实现方案&&代码

先线性筛求出欧拉函数，再把d从1~N for一边即可求出答案

代码：
bzoj ：   M:6876 kb    T:4 ms     L:931 B

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int limit , pri[300005] , phi[1000005] , cntp ;
bool is_notp[1000005] ;

void eular () {
    phi[1] = 1 ;
    for( int i = 2 ; i <= limit ; ++ i ){
        if( !is_notp[i] ){
            pri[ ++cntp ] = i;
            phi[i] = i - 1 ;
        }
        for(int j = 1 ; j <= cntp ; ++ j ){
            if( pri[j] * i > limit )    break ;
            is_notp[ pri[j] * i ] = true ;
            if( i % pri[j] == 0 ){
                phi[ i * pri[j] ] = phi[i] * pri[j] ;
                break;
            }
            else 
                phi[ i * pri[j] ] = phi[i] * ( pri[j] - 1 ) ;
        }
    }
}

int main(){
    int n , m ;
    long long ans = 0 ;
    scanf( "%d%d" , &n , &m ) ;
    limit = min( n , m ) ;
    eular( ) ;
    for( int i = 1 ; i <= limit ; ++ i )
        ans += ((long long)phi[i]) * ( n / i ) * ( m / i ) ;
    printf( "%lld" , (long long)- n * m + ans * 2  ) ;
    return 0 ;
}