【算法】莫比乌斯反演
关注
分享
扫一扫
文章平均质量分 87
cz_xuyixuan
这个作者很懒，什么都没留下…
展开
-
【CodeForces】Educational Codeforces Round 37 题解
【比赛链接】点击打开链接【题解链接】点击打开链接【A】Water The Garden【思路要点】按照题意模拟，或者简单计算一下均可。时间复杂度\(O(N)\)。【代码】#include<bits/stdc++.h>using namespace std;const int MAXN = 5005;template <typename T> void read(T &...原创 2018-02-12 21:08:20 · 470 阅读 · 0 评论 -
【BZOJ3994】【SDOI2015】约数个数和
【题目链接】点击打开链接【思路要点】补档博客，无题解。【代码】#include<bits/stdc++.h>using namespace std;#define MAXN 50005template <typename T> void read(T &x) { x = 0; int f = 1; char c = getchar(); for (; !...原创 2018-05-30 20:36:05 · 213 阅读 · 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 · 931 阅读 · 2 评论 -
【省内训练2018-12-23】String
【思路要点】考虑无问号的情况，分为两种：111 、 S=TS=TS=T ，那么 AAA 和 BBB 取任意字符串均可，贡献为 arbitrary=∑i=1N2i∑j=1N2jarbitrary=\sum_{i=1}^{N}2^i\sum_{j=1}^{N}2^jarbitrary=∑i=1N2i∑j=1N2j 。222 、 S≠TS\ne TS̸=T ，那么要求 AAA 和 BBB...原创 2018-12-23 17:47:10 · 320 阅读 · 0 评论 -
【CodeForces】Hello 2019 (Div. 1 + Div. 2) 题解
【比赛链接】点击打开连接【题解链接】点击打开链接**【A】**Gennady and a Card Game【思路要点】按照题意模拟。时间复杂度 O(1)O(1)O(1) 。【代码】#include<bits/stdc++.h>using namespace std;const int MAXN = 2e5 + 5;typedef lo...原创 2019-01-09 13:01:13 · 957 阅读 · 2 评论 -
【省内训练2019-06-04】树
【思路要点】关键点的限制显然可以对每个非关键点的联通块去除统计的影响。考虑如何计算 ∏i=1N∏j=iNDis(i,j)\prod_{i=1}^{N}\prod_{j=i}^{N}Dis(i,j)i=1∏Nj=i∏NDis(i,j)记 S={x∣x=pq,p is a prime,q∈N+}S=\{x|x=p^q,p\ is\ a\ prime,q\in...原创 2019-06-07 12:04:16 · 223 阅读 · 0 评论 -
【CodeForces】CodeForces Global Round 3 题解
【比赛链接】点击打开连接【题解链接】点击打开链接【A】 Another One Bites The Dust【思路要点】答案为 2c+2min{a,b}+[a≠b]2c+2min\{a,b\}+[a\ne b]2c+2min{a,b}+[a̸=b]。时间复杂度 O(1)O(1)O(1) 。【代码】#include<bits/stdc++.h>...原创 2019-06-03 13:48:09 · 766 阅读 · 0 评论 -
【校内训练2019-07-09】积性函数求和
【思路要点】所求的式子是约数和的形式，考虑枚举一个约数 iii ，计算其被算了多少次，则答案为∑i=1N∑j=1⌊Ni⌋[gcd(i,j)=1]i=∑i=1N∑j=1⌊Ni⌋∑g∣i,g∣jμ(g)i=∑g=1N∑i=1⌊Ng2⌋i⌊Ng2i⌋\sum_{i=1}^{N}\sum_{j=1}^{\lfloor\frac{N}{i}\rfloor}[gcd(i,j)=1]i\\=\sum_{...原创 2019-07-09 14:35:05 · 362 阅读 · 0 评论 -
【BZOJ4407】于神之怒加强版
【题目链接】点击打开链接【思路要点】补档博客，无题解。【代码】#include<bits/stdc++.h>using namespace std;#define MAXN 5000005#define P 1000000007template <typename T> void read(T &x) { x = 0; int f = 1; char c...原创 2018-05-30 18:11:44 · 247 阅读 · 0 评论 -
【BZOJ2190】【SDOI2008】仪仗队
【题目链接】点击打开链接【思路要点】补档博客，无题解。【代码】#include<bits/stdc++.h>using namespace std;#define MAXN 50005template <typename T> void read(T &x) { x = 0; int f = 1; char c = getchar(); for (; !...原创 2018-05-28 17:05:12 · 283 阅读 · 0 评论 -
【BZOJ2154】Crash的数字表格
【题目链接】点击打开链接【思路要点】补档博客，无题解。【代码】#include<bits/stdc++.h>using namespace std;#define MAXN 10000005#define P 20101009template <typename T> void read(T &x) { x = 0; int f = 1; char c ...原创 2018-05-27 16:04:56 · 243 阅读 · 0 评论 -
【BZOJ3309】DZY Loves Math
【题目链接】点击打开链接【思路要点】补档博客，无题解。【代码】#include<bits/stdc++.h>using namespace std;#define MAXN 10000005int tot, prime[MAXN], f[MAXN], miu[MAXN], cnt[MAXN];int g[MAXN], sum[MAXN];bool mark[MAXN];in...原创 2018-03-08 13:00:25 · 275 阅读 · 0 评论 -
【BZOJ2440】【中山市选2011】完全平方数
【题目链接】点击打开链接【思路要点】补档博客，无题解。【代码】#include<bits/stdc++.h>using namespace std;#define MAXN 100005int tot, prime[MAXN], f[MAXN], miu[MAXN];void init() { for (int i = 2; i < MAXN; i++) { if ...原创 2018-03-10 11:26:23 · 207 阅读 · 0 评论 -
【BZOJ1101】【POI2007】Zap
【题目链接】点击打开链接【思路要点】补档博客，无题解。【代码】#include<bits/stdc++.h>using namespace std;#define MAXN 100005int tot, prime[MAXN], f[MAXN], miu[MAXN], sum[MAXN];void init() { for (int i = 2; i < MAXN; i...原创 2018-03-12 13:13:42 · 227 阅读 · 0 评论 -
【BZOJ2301】【HAOI2011】Problem b
【题目链接】点击打开链接【思路要点】补档博客，无题解。【代码】#include<bits/stdc++.h>using namespace std;#define MAXN 100005int tot, prime[MAXN], f[MAXN], miu[MAXN], sum[MAXN];void init() { for (int i = 2; i < MAXN; i...原创 2018-03-12 19:31:16 · 175 阅读 · 0 评论 -
【BZOJ2820】YY的GCD
【题目链接】点击打开链接【思路要点】补档博客，无题解。【代码】#include<bits/stdc++.h>using namespace std;#define MAXN 10000005int tot, prime[MAXN], f[MAXN], cnt[MAXN], type[MAXN];long long sum[MAXN];void init() { for (i...原创 2018-03-13 11:04:17 · 240 阅读 · 0 评论 -
【BZOJ3529】【SDOI2014】数表
【题目链接】点击打开链接【思路要点】补档博客，无题解。【代码】#include<bits/stdc++.h>using namespace std;#define MAXQ 20005#define MAXN 100005struct query {int n, m, a, home; };struct info {long long value; int home; };...原创 2018-03-13 12:41:09 · 226 阅读 · 0 评论 -
【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 · 287 阅读 · 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 · 607 阅读 · 0 评论 -
【BZOJ2005】【NOI2010】能量采集
【题目链接】点击打开链接【思路要点】补档博客，无题解。【代码】#include<bits/stdc++.h>using namespace std;#define MAXN 100005template <typename T> void read(T &x) { x = 0; int f = 1; char c = getchar(); for (; ...原创 2018-05-27 14:42:42 · 288 阅读 · 0 评论 -
【CodeForces1205E】Expected Value Again
【题目链接】点击打开链接【思路要点】一个字符串存在 borderborderborder iii 等价于其存在周期 N−iN-iN−i 。记 p(i)p(i)p(i) 表示是否存在周期 iii ，由期望的线性性，答案即为 ∑i,jE(p(i)p(j))\sum_{i,j}E(p(i)p(j))∑i,jE(p(i)p(j)) 。考虑枚举 i,ji,ji,j ，如何计算 E(p...原创 2019-08-24 15:18:06 · 568 阅读 · 0 评论