【算法】杜教筛
cz_xuyixuan
这个作者很懒,什么都没留下…
展开
-
【BZOJ3930】【CQOI2015】选数
【题目链接】点击打开链接【思路要点】显然问题可以转化为\(K=1\)的形式。那么,我们实际上要求\(\sum_{i_1,i_2,...,i_N=L}^{R}\epsilon(gcd(i_1,i_2,...,i_N))\)。\(=\sum_{i_1,i_2,...,i_N=L}^{R}\sum_{d/i_1,i_2,...,i_N}\mu(d)\)\(=\sum_{d=1}^{R}\mu(d)(\l...原创 2018-04-01 20:11:56 · 235 阅读 · 0 评论 -
【BZOJ4916】神犇和蒟蒻
【题目链接】点击打开链接【思路要点】第一问答案恒为1。考虑第二问,\(S(N)=\sum_{i=1}^{N}\varphi(i^2)=\sum_{i=1}^{N}i*\varphi(i)\)。注意到\(i^2=\sum_{d/i}d*\varphi( \frac{i}{d} )\)。那么\(\sum_{i=1}^{N}i^2=\sum_{i=1}^{N}\sum_{d/i}d*\varphi( \...原创 2018-04-02 15:04:51 · 310 阅读 · 0 评论 -
【BZOJ4652】【UOJ221】【NOI2016】循环之美
【题目链接】BZOJUOJ【思路要点】通过在十进制下找规律,我们发现分数\(\frac{x}{y}\)在\(k\)进制下为纯循环小数当且仅当\(gcd(y,k)=1\)。稍加分析,我们发现上面这一点并不难证明。那么,原题要求的式子应当是\(\sum_{i=1}^{N}\sum_{j=1}^{M}\epsilon(gcd(i,j))*\epsilon(gcd(j,k))\)。\(=\sum_{i=1...原创 2018-04-03 14:29:54 · 536 阅读 · 0 评论 -
【BZOJ4035】【HAOI2015】数组游戏
【题目链接】点击打开链接【思路要点】补档博客,无题解。【代码】#include<bits/stdc++.h>using namespace std;#define MAXN 100005int n, m;int num, mark[MAXN];int a[MAXN], b[MAXN];int get(int x) { if (x <= m) return a[x];...原创 2018-04-12 21:10:21 · 243 阅读 · 0 评论 -
【BZOJ3944】Sum
【题目链接】点击打开链接【思路要点】杜教筛模板题。 时间复杂度为\(O(TN^{\frac{2}{3}})\)。【代码】#include<bits/stdc++.h>using namespace std;#define MAXN 4000005#define MAXM 10005template <typename T> void ...原创 2018-05-25 20:45:49 · 211 阅读 · 0 评论 -
【51Nod1847】奇怪的数学题
【题目链接】 点击打开链接 【思路要点】 令f(i)=(iMin(i))k(i&amp;gt;1)f(i)=(iMin(i))k(i&amp;gt;1)f(i)=(\frac{i}{Min(i)})^k(i&gt;1),即f(i)f(i)f(i)表示iii次大的因子的kkk次方,特别规定f(1)=0f(1)=0f(1)=0。 那么原式即为∑Ni=1∑Nj=1f(gcd(i,j)...原创 2018-08-14 11:08:46 · 863 阅读 · 2 评论