bzoj 3309: DZY Loves Math

最新推荐文章于 2019-05-20 09:40:07 发布

qq_20669971

最新推荐文章于 2019-05-20 09:40:07 发布

阅读量674

点赞数

本文链接：https://blog.csdn.net/qq_20669971/article/details/53726793

版权

题目：http://www.lydsy.com/JudgeOnline/problem.php?id=3309
思路：感觉这个dzy系列好有意思的样子，要努力做一下呀。
考虑首先化简式子得到

\sum i = 1 n \sum j = 1 m f (g c d (i, j))

$\sum^{n}_{i =1}\sum^{m}_{j=1}f(gcd(i,j))$

= \sum T = 1 m i n (n, m) ⌊ n T ⌋ ⌊ m T ⌋ \sum d | T f (d) * μ (T d)

$=\sum^{min(n,m)}_{T=1}\lfloor\frac{n}{T}\rfloor\lfloor\frac{m}{T}\rfloor\sum_{d|T}f(d)*\mu{(\frac{T}{d})}$
实际上我们可以做了但后面那个不是积性函数只能

nlogn $n\log{n}$ 求会T
怎么办呢？我们进一步考虑后面那个函数的意义/贡献
设

g(T)=∑d|Tf(Td)∗μ(d) $g(T) = \sum_{d|T}f(\frac{T}{d})*\mu{(d)}$
那么首先我们注意到

μ $\mu$ 函数不为0只有当

d $d$ 的形式是

d=Πpi|dpi $d =\Pi_{pi|d}pi$
进一步考虑，有一个性质就是

∑i为偶数C(n,i)−∑i为奇数C(n,i)=[n=0] $\sum_{i为偶数}C(n,i) - \sum_{i为奇数}C(n,i) = [n = 0]$
根据这个性质我们可以发现只有当

T=(Πipi)k $T = (\Pi_{i}{pi})^k$ 时

g(T)≠0 $g(T) \neq 0$
而此时

g(T)=−μ(Πipi) $g(T) = -\mu(\Pi_{i}{pi})$
根据这个我们能线性的求出

g $g$ 函数，这样就不会

T $T$ 啦
复杂度：

O(Tn√+n) $O(T\sqrt{n} + n)$

代码：

#include<iostream>
#include<cstring>
#include<string>
#include<cstdio>
#define N 10000007
using namespace std;
typedef long long LL;
int cnt,mu[N],g[N],prime[N];
bool not_prime[N];
void in(int &num){
    char c;
    while (!isdigit(c = getchar()));
    num = (c ^ 48);
    while (isdigit(c = getchar())) num = 10 * num + (c ^ 48);
}

void init(){
    cnt = 0;
    memset(not_prime,0,sizeof(not_prime));
    mu[1] = 1; 
    for (int i = 2;i < N - 1; ++i){
        if (!not_prime[i]) { prime[++cnt] = i; mu[i] = -1; }
        for (int j = 1;j <= cnt; ++j)
          if (prime[j] * i >= N - 1) break;
          else {
            not_prime[prime[j] * i] = 1;
            if (i % prime[j]) mu[i * prime[j]] = -mu[i];
            else { mu[i * prime[j]] = 0; break; }
          }
    }
    g[0] = g[1] = 0;
    for (LL i = 2;i < N - 1; ++i)
     if (mu[i]) {
       for (LL j = i;j < N - 1;j *= i)
         g[j] = -mu[i]; }
    for (int i = 2;i < N - 1; ++i) g[i] += g[i - 1];
}

void DO_IT(){
    int T,n,m;
    LL ans;
    scanf("%d",&T);
    while (T--){
        ans = 0;
        scanf("%d%d",&n,&m);
        if (n > m) swap(n,m);
        for (int i = 1,j;i <= n;i = j + 1){
            j = min(n / (n / i),m / (m / i));
            ans += 1LL * (n / i) * (m / i) * (g[j] - g[i - 1]);
        }
        printf("%lld\n",ans);
    }
}

int main(){
    init();
    DO_IT();
    return 0;
}