题意 给出$n \leq {10} ^ {12}$，求 $$ \sum_{a = 1} ^ n \sum_{b = 1} ^ n \sum_{c = 1} ^ n [\frac{1}{a} + \frac{1}{b} = \frac{1}{c}] [\gcd(a, b, c) = 1] $$ 题解 ...
分类：
其他好文 时间：
2019-10-15 09:36:22
阅读次数：
89
"bzoj" 题目要的连通块个数可以表示为点数$ $所有生成树上的边数.考虑这个生成树边数,我们维护编号最大生成树,按照编号加入边,然后如果加的时候会成环就把环上编号最小的边挤掉,并且当前的第$i$条边的前驱边$pre_i$为刚才被挤掉的第$j$条边,如果没有前驱边就是0 然后对于一个询问,我们只把 ...
分类：
其他好文 时间：
2019-10-13 20:49:19
阅读次数：
74
神炎皇： 题意： 对于一个整数对$(a,b)$，若满足$a+b<=n$且$a+b$是$a*b$的因子，则成为神奇的数对。请问这样的数对共有多少个？（$N<=10^{14}$） 题解： 已知$a+b<=n\\ (a+b)|ab$。 设$d=\gcd(a,b),x=a/d,y=b/d$。 上式为$(x+ ...
分类：
其他好文 时间：
2019-10-13 19:12:33
阅读次数：
92
发现不码题解还是记不清题。 A. 木板 枚举$y_E$，求出$x_F$关于$y_E$的式子，设$y_E$为$x$，发现$Ans=\sum\limits_{x=1}^{n-1} [n|x^2]$ 考场上受《神炎皇》启发，提出$gcd$，设$gcd(x,n)=d$ $n'd|x'^2d^2$ $n'|x ...
分类：
其他好文 时间：
2019-10-13 10:28:06
阅读次数：
66
沙雕的gcd: //1.gcd int gcd(int a,int b) { return b?gcd(b,a%b):a; } 还有个exgcd： //2.扩展gcd )extend great common divisor ll exgcd(ll l,ll r,ll &x,ll &y) { if( ...
分类：
其他好文 时间：
2019-10-12 21:10:22
阅读次数：
78
一.O(logn)代码小证明 我们先来看下面一段代码: 2. 欧几里得算法 3.幂运算 四.$$库里的log函数 在$$库里有log()函数和log2()函数 log()函数的底数默认为自然对数的底数e log2()函数的底数很显然就是2咯qwq include include include in ...
分类：
编程语言 时间：
2019-10-12 10:54:51
阅读次数：
316
T1： 枚举$m$的个数，$O(1)$算出有几个$x$符合条件。 这样不仅效率低下，还会算重。 把$m$一定时的所有结果拍在数轴上发现仅当$ym>=lcm(n,m)$时会算重。 枚举到$lcm$即可。 时间复杂度$O(n)$。 T2： 直接统计复杂度太高，考虑换一个思路。 枚举$gcd$，将所有边权 ...
分类：
其他好文 时间：
2019-10-11 18:14:49
阅读次数：
68
题目描述 神犇YY虐完数论后给傻×kAc出了一题 给定N, M,求1 using namespace std; typedef long long ll; const int maxn = 1e7; int prime[maxn+10]; int vis[maxn+10]; int mu[maxn+ ...
分类：
其他好文 时间：
2019-10-11 12:47:39
阅读次数：
58
1 // 2 /*快速幂*/ 3 inline int qpow(re int x,re int y,re int res=1){ 4 for(;y;y>>=1,x=x*x%mod) if(y&1) res=res*x%mod; return res; 5 } 6 // 7 /*gcd*/ 8 //... ...
分类：
其他好文 时间：
2019-10-10 22:52:18
阅读次数：
101
1、如果$gcd(i,j)==1$，且$i+j==k$，那么这样的数对数就是$\phi(k)$。 也就是$gcd(i,j)==1$导出$gcd(i,k-i)==1$，进而$gcd(i,k)==1$，从而转化为$euler$。 2、https://www.cnblogs.com/henry-1202/ ...
分类：
其他好文 时间：
2019-10-10 22:47:03
阅读次数：
100
