2018 ICPC 徐州计蒜客 - Easy Math

最新推荐文章于 2021-08-10 11:53:05 发布

佐理慧

最新推荐文章于 2021-08-10 11:53:05 发布

阅读量412

点赞数 1

分类专栏：刷题小结反演与容斥杜教筛数值计算

本文链接：https://blog.csdn.net/ZLH_HHHH/article/details/82562820

版权

刷题小结同时被 3 个专栏收录

68 篇文章 1 订阅

订阅专栏

反演与容斥

20 篇文章 0 订阅

订阅专栏

杜教筛

12 篇文章 0 订阅

订阅专栏

2018 ICPC 徐州计蒜客 - Easy Math

题目给定 $m<2*10^9,n<10^{12}$ 计算：

a n s = \sum i = 1 m μ (i n)

$ans=\sum_{i=1}^m\mu(in)$
显然：

μ(n)=0 μ ( n ) = 0 $\mu(n)=0$ 时

ans=0 a n s = 0 $ans=0$
当

mu(n)!=0 m u ( n ) ! = 0 $mu(n)!=0$ 时：

a n s = \sum i = 1 m μ (i n) = μ (n) \sum i = 1 m μ (i) [g c d (i, n) = 1]

$ans=\sum_{i=1}^m\mu(in)=\mu(n)\sum_{i=1}^m\mu(i)[gcd(i,n)=1]$
记:

f (k) = \sum i = 1 m μ (i) [g c d (i, n) = = k] F (k) = \sum k | d f (d)

$f(k)=\sum_{i=1}^m\mu(i)[gcd(i,n)==k]\\F(k)=\sum_{k|d}f(d)$
显然,这里

F(k) F ( k ) $F(k)$ 代表下指标遍历所有

gcd g c d $gcd$ 是

k k $k$ 的倍数的情况，那么

k

$k$ 必须是

n n $n$ 的约数才有意义：

F (k) = [k | n] \sum_{i = 1}^{⌊ \frac{m}{k} ⌋} μ (i k)

$F(k)=[k|n]\sum_{i=1}^{\lfloor\frac{m}{k}\rfloor}\mu(ik)$
反演有：

f (k) = \sum k | d μ (d k) F (d)

$f(k)=\sum_{k|d}\mu(\frac{d}{k})F(d)$
则：

a n s = μ (n) f (1) = μ (n) \sum d = 1 m μ (d) [d | n] \sum i = 1 ⌊ m d ⌋ μ (i d) = μ (n) \sum d | n μ (d) \sum i = 1 ⌊ m d ⌋ μ (i d)

$ans=\mu(n)f(1)=\mu(n)\sum_{d=1}^m\mu(d)[d|n]\sum_{i=1}^{\lfloor\frac{m}{d}\rfloor}\mu(id)\\=\mu(n)\sum_{d|n}\mu(d)\sum_{i=1}^{\lfloor\frac{m}{d}\rfloor}\mu(id)$
记：

A (m, n) = a n s = \sum i = 1 m μ (i n)

$A(m,n)=ans=\sum_{i=1}^m\mu(in)$

显然：

$A (m, n) = μ (n) \sum d | n μ (d) A (⌊ m d ⌋, d) A (0, k) = 0 A (k, 1) = M (k) A (1, k) = μ (k)$ $A(m,n)=\mu(n)\sum_{d|n}\mu(d)A(\lfloor\frac{m}{d}\rfloor,d)\\A(0,k)=0\\A(k,1)=M(k)\\A(1,k)=\mu(k)$

其中， $M(n)=\sum_{i=1}^n\mu(i)$ 计算方法参见：https://blog.csdn.net/ZLH_HHHH/article/details/77860249
在计算 $M$ 的同事保留计算A所需的计算信息
计算A的时间可以表示为：

T (m, n) = \sum_{d | n} T (⌊ \frac{m}{d} ⌋, d) < O (\sqrt{m} l o g n)

$T(m,n)=\sum_{d|n}T(\lfloor\frac{m}{d}\rfloor,d)<O(\sqrt{m}log\ n)$
总时间复杂度不超过：

O (m - - \sqrt l o g n) + O (m 0.75)

$O(\sqrt{m}log\ n)+O(m^{0.75})$

#include <algorithm>
#include <stdio.h>
#include <string.h>
#include <cmath>
#define MAXN 2000005
using namespace std;
typedef long long LL;
LL M[MAXN] = { 0,-1 };
LL tmp[MAXN];
LL mu[MAXN + 2];
void clat(LL n, int d)
{
    LL ans = 0;
    int m = (int)sqrt(n) + 1;
    for (int L = 1;L<m;L++)
        ans += M[L] * (n / L - n / (L + 1));
    for (int i = (int)(n / m);i>1;i--)
    {
        LL u = n / i;
        if (u<MAXN)
            ans += M[u];
        else
            ans += tmp[i*d];
    }
    tmp[d] = 1 - ans;
}

LL slove(LL n)
{
    if (n<MAXN)return M[n];
    for (int i = (int)(n / (MAXN - 1));i;i--) clat(n / (LL)i, i);
    return tmp[1];
}

LL DFS(LL m, LL c, LL n, LL mn)
{
    if (n == 1)
    {
        if (m < MAXN)
            return M[m];
        else
            return tmp[c];
    }
    if (m == 1)return mn;
    if (m == 0)return 0;
    LL ans = 0;
    for (LL d = 1;d*d < n;d++)
    {
        if (n%d)continue;
        ans += mu[d] * DFS(m / d, c*d, d, mu[d]);
        int u_ = mn / mu[d];
        LL d_ = n / (LL)d;
        ans += u_*DFS(m / d_, c*d_, d_, u_);
    }
    return ans*mn;
}

int main()
{
    for (int i = 1;i<MAXN;i++)
    {
        M[i] = -M[i];
        mu[i] = M[i];
        for (int j = i << 1;j<MAXN;j += i)M[j] += M[i];
        M[i] += M[i - 1];
    }
    LL n, m;
    scanf("%lld %lld", &m, &n);
    LL mn = 1;
    slove(m);
    if (n >= MAXN)
    {
        LL a = n;
        for (LL i = 2;i*i <= n;i++)
        {
            if (a%i == 0)
            {
                a /= i;
                mn = -mn;
                if (a%i == 0)
                {
                    mn = 0;
                    break;
                }
            }
        }
        if (mn != 0 && a > 1)mn = -mn;
    }
    else
        mn = mu[n];
    if (mn == 0)
    {
        printf("0\n");
        return 0;
    }
    printf("%lld\n", DFS(m, 1, n, mn));
    return 0;
}