赛上想写,Track Lost 了属于是. \(\mathscr{Intro}\)   Min_25 筛是用于求积性函数前缀和,同时顺带求出一些"有意思"的信息的筛法.   一些记号约定 \(\mathbb P\) 为素数集,对于以 \(p\) 为记号的数,有 \(p\in\mathbb P\). \(p_i\) 表示第 \(i\) 小的素数.特别地,\(p_0=1\). \(\newcommand{\lpf}[0]{\operatorname{lpf}} \lpf(n)\) 表示…
目录 Preface 数论函数 积性函数 Dirichlet 卷积 Dirichlet 卷积中的特殊函数 Mobius 函数 & Mobius 反演 Mobius 函数 Mobius 反演 基础应用 约数个数 欧拉函数 反演魔法 例一 例二 例三 魔法中的 tricks 线性筛 trick 筛 筛 筛 刷表 trick Conclusion   UPD:修改了 Euler 筛法代码框架.   若无特别说明,\(x,y\) 等形式变量均 \(\in\mathbb N_+\):\(p\) 为素数.…
目录 问题引入 思考 Lagrange 插值法 插值过程 代码实现 实际应用 「洛谷 P4781」「模板」拉格朗日插值 「洛谷 P4463」calc 题意简述 数据规模 Solution Step 1 Step 2 证明 代码 「CF 995F」Cowmpany Cowmpensation 题意简述 数据规模 Solution Step 1 Step 2 证明 代码 「CF 662F」The Sum of the k-th Powers 题意简述 数据规模 Solution 代码 「BZOJ 3…
目录 「CF 750E」New Year and Old Subsequence 「洛谷 P4719」「模板」"动态 DP" & 动态树分治 「洛谷 P6021」洪水 「SP 6779」GSS7 「NOIP 2018」「洛谷 P5024」保卫王国 \(\mathcal{Introduction}\) \(\mathcal{Problem~1}\)   给定序列 \(\{a_n\}\),其中 \(a_i\in\mathbb Z\),求其最大子段和(不能为空).   很显然的 DP…
\(\mathcal{Preface}\)   单位根反演,顾名思义就是用单位根变换一类式子的形式.有关单位根的基本概念可见我的这篇博客. \(\mathcal{Formula}\)   单位根反演的公式很简单: \[[k|n]=\frac{1}k\sum_{i=0}^{k-1}\omega_k^{ni} \] \(\mathcal{Proof}\)   分类讨论: \(k|n\). 那么 \((\forall i)(\omega_k^{ni}=1)\),所以右侧为 \(\frac{1}k\su…
题目链接 https://loj.ac/problem/3069 题解 复数真神奇. 一句话题意:令 \(f(x)\) 表示以原点 \((0, 0)\) 为圆心,半径为 \(x\) 的圆上的整点数量,求 \(\sum_\limits{i = 1}^N f(i)^k \bmod P\) 的值. 令 \(g(x) = \frac{f(x)}{4}\),那么我们需要求 \(\left(4^k\sum_\limits{i = 1}^Ng(i)^k\right) \bmod P\).打表可得 \(g(x)…
传送门 思路 (以下令\(F(n)=f(n)^k\)) 首先肯定要莫比乌斯反演,那么可以推出: \[ ans=\sum_{T=1}^n \lfloor\frac n T\rfloor^2\sum_{d|T}F(d)\mu(T/d) \] 可以整除分块,但后面的东西怎么办呢? 令\(G(T)=F*\mu\),那么就有 \[ ans=\sum_{T=1}^n \lfloor\frac n T\rfloor^2G(T) \] 看到\(\mu\)函数有点烦,考虑用杜教筛的式子消去它. \[ g(1)S(…
传送门. 题解 感觉这题最难的是第一个结论. x/y首先要互质,然后如果在10进制是纯循环小数,不难想到y不是2.5的倍数就好了. 因为十进制下除以2和5是除得尽的. 必然会多出来的什么东西. 如果是k进制,可以类比得gcd(y,k)=1. 证明: 假设纯循环的位数是l 则\(x*k^l\equiv x(mod~y)\) \(k^l\equiv 1(mod~y)\) 要存在l的话,就必须有\(gcd(k,y)=1\),反过来一样. 反演: \(Ans=\sum_{i=1}^n\sum_{j=1}…
\(\mathcal{Description}\)   Link.   对于积性函数 \(f(x)\),有 \(f(p^k)=p^k(p^k-1)~(p\in\mathbb P,k\in\mathbb N_+)\).求 \(\sum_{i=1}^nf(i)\bmod(10^9+7)\).   \(n\le10^{10}\). \(\mathcal{Solution}\)   Min_25 筛是不可能的.   Powerful Number 三步走咯!考虑素数点值: \[f(p)=p^2-p \]…
「学习笔记」Min25筛 前言 周指导今天模拟赛五分钟秒第一题,十分钟说第二题是 \(\text{Min25}​\) 筛板子题,要不是第三题出题人数据范围给错了,周指导十五分钟就 \(\text{AK}​\) 了,为了向 \(\text{AK}​\)王 学习,真诚的膜拜他,接受红太阳的指导,下午就学习了一下 \(\text{Min25}​\) 筛. 简介 如果 \(f(n)\) 是一个积性函数,且 \(f(n)\) 是一个关于 \(n\) 的简单多项式,并可以快速算出 \(f(p^k),\ p\…