- 博客(8)
- 资源 (1)
- 收藏
- 关注
RSS订阅
原创 uestc 811 GCD 杜教筛 + 自然幂和
传送门:UESTC 811这题思路出的还算快, 但是UESTC OJ评测卡了一天, 代码细节比较多, 调bug的时候都怀疑解法是不是又错了这题好像不好找题解题意1<=N<=10^10, 1<=K<=5求∑Ni=1∑Nj=1gcd(i,j)k\sum_{i=1}^N\sum_{j=1}^Ngcd(i, j)^k题解枚举gcd, 推下式子 ans=∑i=1N∑j=1Ngcd(i,j)k=2∑i=1N
2017-08-22 23:16:00
1067
3
原创 HDU 5608 莫比乌斯反演 + 莫比乌斯函数前缀和
传送门: HDU 5608 题解: 令 G（n） = n ^ 2 - 3 * n + 2 先反演得:f(n)=∑d|nG(d)∗u(nd)f(n) = \sum_{d|n}G(d) * u(\frac{n}{d})令: A(n)=∑i=1nf(i)A(n) = \sum_{i = 1}^nf(i) 则:A(n)=∑i=1nf(i)=∑i=1n∑d|iG(d)∗u(id)A(n) = \sum_
2017-08-06 18:15:37
1217
3
原创 codeforces 837E
题解首先我们知道 gcd(a,b)=gcd(a,a+b)f(a,b)=f(a,b−gcd(a,b))gcd(a, b) = gcd(a, a + b) \\f(a, b) = f(a, b - gcd(a, b))\\令b=k∗a+mb=k*a+m 则gcd(a,b)=gcd(a,m)gcd(a, b) = gcd(a, m) f(a,b)=f(a,k∗a+m−gcd(a,m))=k+f(
2017-08-30 13:48:47
887
原创 2017 CCPC 网络赛 1007
链接:HDU 6156题意给定l, r, L, R, 求∑Ri=L∑rj=lf(i,j)\sum_{i=L}^{R}\sum_{j=l}^{r}f(i,j) 其中f(i,j)={1[i在j进制下不是回文数]j[i在j进制下是回文数]f(i, j) = \left\{\begin{matrix}1[i在j进制下不是回文数] \\j[i在j进制下是回文数]\end{matrix} \righ
2017-08-20 14:08:47
1793
原创 51nod 1227
杜教筛#include <bits/stdc++.h>using namespace std;typedef long long ll;const int N = 1000001;const ll mod = 1e9 + 7;const ll inv = (mod + 1) / 2;const ll _6 = (mod + 1) / 6;const int mo = 2333333;bo
2017-08-07 18:52:24
384
原创 51nod 1237 杜教筛
题解参考1238code:#include <bits/stdc++.h>using namespace std;typedef long long ll;const int N = 1000001;const ll mod = 1e9 + 7;const ll inv = (mod + 1) / 2;const int mo = 2333333;bool isPrime[N];ll p
2017-08-07 18:51:14
428
原创 51nod 1238 杜教筛
传送门：51nod 1238题意求G(N)=∑i−1N∑j=1Nlcm(i,j) G(N) = \sum_{i - 1}^N\sum_{j=1}^Nlcm(i, j)题解首先G(N)=∑i=1N∑j=1Nlcm(i,j)=2∑i=1N∑j=1ilcm(i,j)−∑i=1Nlcm(i,i)=2∑i=1Ni∑d|i∑u=1idu[gcd(u,id)=1]−N(N+1)2=2∑i=1Ni∑d|iidϕ
2017-08-07 16:37:22
792
原创 HDU 6051原根 + 欧拉函数
题解: 主要是原根的幂表示[1, p - 1], 然后转换 hint里的公式: 找了一篇证明:http://blog.csdn.net/Clove_unique/article/details/53152473代码:#include <cstdio>#include <cstring>#include <algorithm>using namespace std;typedef l
2017-08-02 15:55:40
756
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人