奇怪的数学题(51nod1847)——min_25筛+杜教筛+第二类斯特林数

题面

   51nod1847

解析

   定义$f(i)$为$i$的次大因数,特别地$f(1)=0$

    那么我们就是去求$$\sum_{i=1}^{n}\sum_{j=1}^{n}f^{m}(gcd(i, j))$$

  这种东西的套路就是枚举$gcd$然后用欧拉函数处理, \begin{align*}\sum_{i=1}^{n}\sum_{j=1}^{n}f^{m}(gcd(i, j)) &= \sum_{i=1}^{n}\sum_{j=1}^{\left \lfloor \frac{n}{i}\right \rfloor} \sum_{k=1}^{\left \lfloor \frac{n}{i}\right \rfloor}[gcd(j, k)==1]f(i)\\&=\sum_{i=1}^{n}f^{m}(i)((2\sum_{j=1}^{\left \lfloor \frac{n}{i}\right \rfloor}\varphi(j)) - 1)\end{align*}

  求$\varphi(i)$的前缀和用杜教筛即可,问题在于如何求$f(i)$的前缀和

  考虑$min25$筛的过程,它把所有数对前缀和的贡献分为了质数的贡献与合数的贡献,这道题中质数的贡献显然是1,因此,可以构造一个完全积性函数$c(i)=1$处理质数的贡献。接下来处理合数的贡献,考虑在$min25$的第一步处理过程中,每一个合数都会被它最小的质因子$minp$筛去。因此在枚举第$i$个质数$pri[i]$时就会有以$pri[i]$为最小质因子的所有合数的贡献。定义完全积性函数$b(i)=i^{m}$,$g$数组为筛法第一步中处理的函数,$h(i)$表示$1$到$i$中合数对答案的贡献,用公式写来,就是这样:$$g(n, i) = g(n, i)-b(pri[i])*(g(\left \lfloor  \frac{n}{pri[i]}\right \rfloor, i-1) - \sum_{j=1}^{i-1}b(pri[j]))\\h(n)=h(n)+g(\left \lfloor  \frac{n}{pri[i]}\right \rfloor, i-1) - \sum_{j=1}^{i-1}b(pri[j])$$

  $b(i)$的前缀和在线性筛时预处理即可,问题在于如何预处理$g(n, 0)=\sum_{i=1}^{n}i^{m}$,考虑用第二类斯特林数,原式变为:$$\sum_{i=1}^{n}\sum_{j=0}^{m}\begin{Bmatrix}m\\j\end{Bmatrix}\cdot j!\cdot \binom{i}{j}$$

  上式$j$可以直接枚举到$m$, 若$j$大于$m$,$\begin{Bmatrix}m\\j\end{Bmatrix}=0$;若$j$大于$i$,$C_{i}^{j}=0$,两种情况均不会对求和造成影响。

  套路交换求和号,变为:$$\sum_{j=0}^{m}\begin{Bmatrix}m\\j\end{Bmatrix}\cdot j!\sum_{i=0}^{n}\binom{i}{j}$$

  注意到$\sum_{i=0}^{n}\binom{i}{j}=\binom{n+1}{j+1}$,代入得:$$\sum_{j=0}^{m}\begin{Bmatrix}m\\j\end{Bmatrix}\cdot j! \cdot \binom{n+1}{j+1}$$

  展开$C_{i}^{j}=0$, 将$j!$代入,最终变为$$\sum_{j=0}^{m}\begin{Bmatrix}m\\j\end{Bmatrix}\frac{\prod_{i=n-j+1}^{n+1}i}{j+1}$$

  注意$j+1$可能没有逆元,需要先把$j+1$除了才能取模

  最后做一次数论分块统计答案

  在做这道题时才意识到的细节,在数论分块时,存下来的值恰好为每一块区间的右端点...

 代码:

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned int ui;
const int maxn = 1000005;

int n, m, mp[maxn], Mp[maxn], a[maxn];
ui s2[60][60], b[maxn], c[maxn];

ll pri[maxn];
int phi[maxn], cnt;
ui f[maxn], g[maxn], s[maxn];
bool notp[maxn];
map<int, ui> F;

ui qpow(ui x, int y)
{
    ui ret = 1;
    while(y)
    {
        if(y&1)
            ret = ret * x;
        x = x * x;
        y >>= 1;
    }
    return ret;
}

void Euler()
{
    phi[1] = 1;
    for(int i = 2; i <= 1000000; ++i)
    {
        if(!notp[i])
        {
            pri[++cnt] = i;
            phi[i] = i - 1;
            s[cnt] = s[cnt-1] + qpow(i, m);
        }
        for(int j = 1; j <= cnt; ++j)
        {
            if(pri[j] * i > 1000000)    break;
            notp[pri[j]*i] = 1;
            if(i % pri[j] == 0)
            {
                phi[pri[j]*i] = pri[j] * phi[i];
                break;
            }
            phi[pri[j]*i] = (pri[j] - 1) * phi[i];
        }
    }
    for(int i = 1; i <= 1000000; ++i)
        f[i] = f[i-1] + phi[i];
}

void init()
{
    s2[0][0] = 1;
    for(int i = 1; i <= m; ++i)
        for(int j = 1; j <= i; ++j)
            s2[i][j] = s2[i-1][j-1] + j * s2[i-1][j];
}

int get_id(int x)
{
    return x <= 1000000? mp[x]: Mp[n/x];
}

int gcd(int a, int b)
{
    return b == 0? a: gcd(b, a % b);
}

int stak[60];

ui get_mi(int x)
{
    ui ret = 0, mul;
    int sj = min(x, m), t, gc;
    for(int i = 1; i <= sj; ++i)
    {
        mul = s2[m][i];
        for(int j = 1; j <= i + 1; ++j)
            stak[j] = x + 2 - j;
        t = i + 1;
        for(int j = 1; j <= i + 1; ++j)
        {
            gc = gcd(stak[j], t);
            stak[j] /= gc;
            t /= gc;
            if(t == 1)    break;
        }
        for(int j = 1; j <= i + 1; ++j)
            mul = mul * stak[j];
        ret += mul;
    }
    return ret;
}

ui dfs(int x)
{
    if(x <= 1000000)    return f[x];
    if(F.find(x) != F.end())    return F[x];
    ll tmp = (1LL * x * (x + 1)) >> 1;
    ui ret = (ui)tmp;
    for(int l = 2, r; l <= x; l = r + 1)
    {
        r = x / (x / l);
        ret -= (r - l + 1) * dfs(x / l);
    }
    return F[x] = ret;
}

int main()
{
    scanf("%d%d", &n, &m);
    Euler();
    init();
    int t, tot = 0, id, now;
    for(int l = 1, r; l <= n; l = r + 1)
    {
        t = n / l;
        r = n / t;
        a[++tot] = t;
        if(t <= 1000000)
            mp[t] = tot;
        else
            Mp[n/t] = tot;
        b[tot] = get_mi(t) - 1;
        c[tot] = t - 1;
    }
    for(int i = 1; i <= cnt && pri[i] * pri[i] <= n; ++i)
        for(int j = 1; j <= tot && pri[i] * pri[i] <= a[j]; ++j)
        {
            id = get_id(a[j] / pri[i]);
            b[j] -= (b[id] - s[i-1]) * (s[i] - s[i-1]);
            g[j] += b[id] - s[i-1];
            c[j] -= c[id] - i + 1;
        }
    ui ans = 0;
    for(int l = 1, r; l <= n; l = r + 1)
    {
        r = n / (n / l);
        now = get_id(r);
        id = get_id(l - 1);
        ans += ((g[now] + c[now]) - (g[id] + c[id])) * (2 * dfs(n / l) - 1);
    }
    printf("%u", ans);
    return 0;
}

原文地址:https://www.cnblogs.com/Joker-Yza/p/11969074.html

时间: 2024-09-28 10:12:55

奇怪的数学题(51nod1847)——min_25筛+杜教筛+第二类斯特林数的相关文章

EOJ Monthly 2019.11 E. 数学题(莫比乌斯反演+杜教筛+拉格朗日插值)

传送门 题意: 统计\(k\)元组个数\((a_1,a_2,\cdots,a_n),1\leq a_i\leq n\)使得\(gcd(a_1,a_2,\cdots,a_k,n)=1\). 定义\(f(n,k)\)为满足要求的\(k\)元组个数,现在要求出\(\sum_{i=1}^n f(i,k),1\leq n\leq 10^9,1\leq k\leq 1000\). 思路: 首先来化简一下式子,题目要求的就是: \[ \begin{aligned} &\sum_{i=1}^n\sum_{j=1

hdu6607 min25筛+杜教筛+伯努利数求k次方前缀和

推导过程类似https://www.cnblogs.com/acjiumeng/p/9742073.html 前面部分min25筛,后面部分杜教筛,预处理min25筛需要伯努利数 //#pragma GCC optimize(2) //#pragma GCC optimize(3) //#pragma GCC optimize(4) //#pragma GCC optimize("unroll-loops") //#pragma comment(linker, "/stack

数论入门——莫比乌斯函数,欧拉函数,狄利克雷卷积,线性筛,莫比乌斯反演,杜教筛

一个菜鸡对数论的一点点理解... 莫比乌斯函数 定义函数\(\mu(n)\)为: 当n有平方因子时,\(\mu(n)=0\). 当n没有平方因子时,\(\mu(n)=(-1)^{\omega(n)}\),\(\omega(n)\)表示n不同质因子的个数. 性质1: \(\sum_{d|n}\mu(d)=[n=1]\) 证明:我们把n分解质因数,则原式\(=(-1+1)^{\omega(n)}=0\). 因为对于不同的质因子,只有选和不选两种方案,这是一个组合数相加的形式,偶数加奇数减,根据二项式

【数论】狄利克雷卷积及其快速计算方法及杜教筛

目录(假的 狄利克雷卷积基础知识 数论函数 狄利克雷卷积定义 狄利克雷卷积性质 常用卷积 卷积计算方法 最暴力的暴力 稍好的暴力 优美的暴力 莫比乌斯反演(待填坑) 杜教筛 经典杜教筛 第二种杜教筛 第三种杜教筛 背景 本人即将去CTS&APIO2019,由于一些特殊原因,发现自己数论突然变得很菜. 就决定在去的前一天,翻出来以前的数论学习资料看一看.翻到了czgj的校内狄利克雷卷积课件,发现其中提到了的任意数列\(f(n)\)和\(g(n)\)的狄利克雷卷积\((f*g)(n)\)(从1到n,

【Luogu3768】简单的数学题(莫比乌斯反演,杜教筛)

[Luogu3768]简单的数学题(莫比乌斯反演,杜教筛) 题面 洛谷 \[求\sum_{i=1}^n\sum_{j=1}^nijgcd(i,j)\] $ n<=10^9$ 题解 很明显的把\(gcd\)提出来 \[\sum_{d=1}^nd\sum_{i=1}^n\sum_{j=1}^nij[gcd(i,j)==d]\] 习惯性的提出来 \[\sum_{d=1}^nd^3\sum_{i=1}^{n/d}\sum_{j=1}^{n/d}ij[gcd(i,j)==1]\] 后面这玩意很明显的来一发

loj#6229. 这是一道简单的数学题 (??反演+杜教筛)

题意:给定\(n\le 10^9\),求:\(F(n)=\sum_{i=1}^n\sum_{j=1}^i\frac{\mathrm{lcm}(i,j)}{\mathrm{gcd}(i,j)}\),对1e9+7取模 推式子: \(F(n)=\sum_{i=1}^n\sum_{j=1}^i\frac{\mathrm{lcm}(i,j)}{\mathrm{gcd}(i,j)}\) \(=\sum_{i=1}^n\sum_{j=1}^i\frac{ij}{\gcd^2(i,j)}\) \(=\sum_{

杜教筛 学习总结

看了看唐老师的blog,照猫画虎的做了几道题目,感觉对杜教筛有些感觉了 但是稍微有一点难度的题目还是做不出来,放假的时候争取都A掉(挖坑ing) 这篇文章以后等我A掉那些题目之后再UPD上去就好啦 由于懒得去写怎么用编辑器写公式,所以公式就准备直接copy唐老师的啦 首先积性函数和完全积性函数什么的就不再多说了 列举常见的积性函数: 1.约数个数函数和约数个数和函数 2.欧拉函数phi 3.莫比乌斯函数mu 4.元函数e 其中e(n)=[n==1] 5.恒等函数I 其中I(n)=1 6.单位函数

【bzoj 4176】 Lucas的数论 莫比乌斯反演(杜教筛)

Description 去年的Lucas非常喜欢数论题,但是一年以后的Lucas却不那么喜欢了. 在整理以前的试题时,发现了这样一道题目“求Sigma(f(i)),其中1<=i<=N”,其中 表示i的约数个数.他现在长大了,题目也变难了. 求如下表达式的值: 一行一个整数ans,表示答案模1000000007的值. Sample Input 2 Sample Output 8 HINT 对于100%的数据n <= 10^9. 题解: 解锁新技能:杜教筛. 再复习一下: 若$F(n)=\s

【51nod-1239&amp;1244】欧拉函数之和&amp;莫比乌斯函数之和 杜教筛

题目链接: 1239:http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1239 1244:http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1244 杜教筛裸题,不过现在我也只会筛这俩前缀和... $$s(n)=\sum _{i=1}^{n}f(i)$$ 那么就有: $$\sum_{i=1}^{n}f(i)\lfloor \frac{n}{i} \