2020年11月_SHOHOKUKU

原创 LeeCode 1515 模拟退火

题意传送门 LeeCode 1515. 服务中心的最佳位置题解二维区域上求到所有客户的欧几里得距离的总和最小位置的服务中心的位置，可以使用模拟退火法求解。在客户区域随机初始化较多的点，那么模拟退火时可以简单地当答案更优时进行转移，而不用依概率转移。class Solution{#define maxn 50#define maxl 100public: struct P { double x, y, d; } ct[maxn], sv[maxn];

2020-11-28 22:41:23 289

原创 P3265 [JLOI2015] 线性基

题意传送门 P3265 [JLOI2015]装备购买题解购买装备的最大数量即线性空间的基，高斯消元求矩阵的秩即可。对于最小的花费，在高斯消元的时候采取贪心策略，取满足当前主元位置非零的且花费最小的行进行消元。简单地证明，线性基与任一向量线性相关，那么可以用一个更小的向量替换其线性基表示中非零系数的一个向量。#include <bits/stdc++.h>using namespace std;#define maxn 505#define eps 1e-5double A[ma

2020-11-27 17:25:18 82

原创 POJ 1830 高斯消元

题意传送门 POJ 1830题解设灯 iii 是否开关为 xix_ixi，灯 jjj 操作时是否会影响到灯 iii 为 ai,ja_{i,j}ai,j，设灯 iii 初始与结束的开关状态分别为 si,tis_i,t_isi,ti，则灯 iii 满足XORj∈[1,N](xj & ai,j)=si⊕tiXOR_{j\in [1,N]}(x_j\ \&\ a_{i,j})=s_i\oplus t_iXORj∈[1,N](xj &&nbsp

2020-11-26 23:33:26 93 1

原创 P4035 [JSOI2008] 高斯消元

题意传送门 P4035 [JSOI2008]球形空间产生器题解设球性坐标为 (x0,x1,…,xn−1)(x_0,x_1,\dots,x_{n-1})(x0,x1,…,xn−1)，在球面上的点 iii 坐标为 (ai,0,ai,1,…,ai,n−1)(a_{i,0},a_{i,1},\dots,a_{i,n-1})(ai,0,ai,1,…,ai,n−1)，满足∑j=0n−1(ai,j−xj)2=R2\sum\limits_{j=0}^{n-1}(a_{i,j}-x_j)^2=R^2j=0

2020-11-26 21:21:13 51

原创 POJ 2891 线性同余方程组

题意传送门 POJ 2891题解线性同余方方程，可以使用拓展欧几里德算法计算逆元求解。对于线性同余方程组，可以通过方程数次拓展欧几里德算法迭代求解。假设方程形式为 ai×x = bi(mod mi)a_i\times x\ = \ b_i(mod\ m_i)ai×x = bi(mod mi)，在题目中为 x = ri(mod ai)x\ = \ r_i(mod\ a_i)x = ri(

2020-11-26 20:26:27 236

原创 POJ 3696 同余 + 欧拉函数

题意传送门 POJ 3696题解8…8=8×1…1=8×(10x−1)/98\dots8=8\times 1\dots 1=8\times (10^x-1)/98…8=8×1…1=8×(10x−1)/9，目标是求解满足下式的最小 xxx8(10x−1)9 mod L=0\frac{8(10^x-1)}{9}\ mod\ L=098(10x−1) mod L=0 对于分母 999，从约数的角度考虑，将其合并进模数；对于分子 888，设 d=gcd(8,L)d=

2020-11-26 16:49:53 116

原创 POJ 3090 欧拉函数

题意传送门 POJ 3090题解可视的点即满足 gcd(x,y)=1gcd(x,y)=1gcd(x,y)=1 的点，因为此时在斜率为 y/xy/xy/x 的射线上，满足以其为斜边，平行于 x,yx,yx,y 轴的边为直角边的直角三角形中，直角边为整数的在射线上非原点的顶点离原点最近的点就是 (x,y)(x,y)(x,y)。观察到这样的点关于 y=xy=xy=x 对称，那么仅处理 y<xy<xy<x 的部分即可，特殊处理 y=xy=xy=x 与 x,yx,yx,y 为 000 的点，

2020-11-26 12:17:02 160

原创 P1072 质因数

题意传送门 P1072 Hankson 的趣味题题解求满足下列方程组的 xxx 个数。 {gcd(x,a)=blcm(x,c)=d\begin{cases}gcd(x,a)=b\\ lcm(x,c)=d \\ \end{cases}{gcd(x,a)=blcm(x,c)=d 考虑质因数 ppp，设 a,b,c,d,xa,b,c,d,xa,b,c,d,x 中因数 ppp 的个数分别为 ma,mb,mc,md,mxm_a,m_b,m_c,m_d,m_xma,mb,mc,md,mx。那么有{

2020-11-26 11:48:53 383

原创 POJ 2054 贪心

题意传送门 POJ 2054 Color a Tree题解本质是在树上遍历顺序的限制下，找到一个染色顺序 chs[i]chs[i]chs[i]，使数列 {1,2,…,n}\{1,2,\dots,n\}{1,2,…,n} 与 {C1,C2,…,Cn}\{C_1,C_2,\dots,C_n\}{C1,C2,…,Cn} 在对应位置相乘，即 ∑i=1ni×Cchs[i]\sum\limits_{i=1}^{n}i\times C_{chs[i]}i=1∑ni×Cchs[i] 取到最小值。假设没有限制

2020-11-26 00:06:38 132

原创 Codeforces 452E 后缀数组 + 并查集

题意传送门 Codeforces 452E. Three strings题解将子串看做后缀的前缀，那么使用两个互不相同且非小写字母的字符将 333 个字符串拼接起来，求解后缀数组以及高度数组。由于后缀数组的有序性，故最长相同前缀不小于某个长度 hhh 的后缀集合在后缀数组中处于连续的位置，且对应位置的高度数组的值不小于 hhh。那么可以使用并查集维护这样连续的位置，从高度数组的最大值（333 个字符串的最长长度）开始处理，使用相邻 sa[i],sa[i+1]sa[i],sa[i+1]sa[i],sa

2020-11-25 21:26:50 155

原创 P2261 [CQOI2007] 整除分块

题意传送门 P2261 [CQOI2007]余数求和题解取模运算可以写为 k mod i=k−⌊k/i⌋×ik\ mod\ i = k-\lfloor k/i\rfloor\times ik mod i=k−⌊k/i⌋×i，那么有G(n,k)=∑i=1n(k−⌊k/i⌋×i)=nk−∑i=1n(⌊k/i⌋×i)G(n,k)=\sum\limits_{i=1}^{n}(k-\lfloor k/i\rfloor\times i)=nk-\sum\limits_{i

2020-11-24 19:51:49 72

原创 P1463 [POI2002][HAOI2007] 因数 + DFS

题意传送门 P1463 [POI2002][HAOI2007]反素数题解答案为 [1,N][1,N][1,N] 中约数个数最多的数中最小的一个。考虑数据范围内 x∈[1,N]x\in [1,N]x∈[1,N] 的约数个数，考虑最小的素数 222，则 xxx 的质因子指数不可能超过 303030；考虑最小的连续素数，则 xxx 的质因数个数不可能超过 999。于是有了搜索范围的初步限制。更进一步的，xxx 为反素数的必要条件是：将 xxx 分解质因数后可写作 2e1×3e2…(e1≥e2⋯≥0)2^{

2020-11-24 12:49:09 100

原创 POJ 3415 后缀数组 + 单调栈 / 并查集

题意传送门 POJ 3415题解子串是原串中连续的一段，也可以定义为前缀的后缀或后缀的前缀。统计分别属于 A,BA,BA,B 的不小于 KKK 的子串个数，那么将 A,BA,BA,B 用一个不属于这两个串的字符拼接起来（避免拼接位置对结果产生影响），构造后缀数组以及高度数组（lcp[i]lcp[i]lcp[i] 为 a[i],sa[i+1]a[i],sa[i+1]a[i],sa[i+1] 的最长公共前缀）。设任意两个后缀在后缀数组中的索引分别为 i,j(i<j)i,j(i<j)i,j

2020-11-24 00:07:26 336

原创 POJ 1509 后缀数组

题意传送门 POJ 1509题解求长度为 nnn 的字符串构成的环任一位置拆分成的 nnn 条链中，字典序最小的，若有多个答案，则返回拆开位置索引最小的。将原串复制一份连接到自身后面，求后缀数组，实现上使用基于倍增的方法。答案为满足长度大于等于 nnn 的后缀的长度为 nnn 的前缀。为了满足拆开位置索引最小的要求，有两种思路：求解后缀数组时，将超出后缀长度的部分赋最大的权值，此时保证按字典序排序时，较短串为较长串前缀时，较短串字典序比较长串更大；或者直接在复制拼接的字符串后面添加一个字典序最大的字

2020-11-21 20:47:09 133

原创 POJ 2018 二分

题意传送门 POJ 2018题解最大化平均值问题考虑二分答案，问题转化为判定是否存在一个长度大于等于 FFF 的子区间，使其平均数不小于二分值。设区间为 [l,r][l,r][l,r]，当前二分值为 xxx，则判定不等式为∑i∈[l,r]C[i]r−l+1≥x\frac{\sum_{i\in [l,r]}C[i]}{r-l+1}\geq xr−l+1∑i∈[l,r]C[i]≥x 左右两边乘上分母，移项后得到∑i∈[l,r](C[i]−x)≥0\sum_{i\in[l,r]}(C[i]-x)\g

2020-11-20 18:59:39 251

原创 Educational Codeforces Round 98 (Rated for Div. 2)

传送门Educational Codeforces Round 98 (Rated for Div. 2)A目标是从 (0,0)(0,0)(0,0) 走到 (x,y)(x,y)(x,y)，横纵坐标的行进需要尽可能地交叉，以减少原地不动的命令次数。假如 x=yx=yx=y，那么无需原地不动的命令；否则设其中较大值为 xxx，那么需要在 x−yx-yx−y 个横轴或纵轴方向的行进中插入 x−y−1x-y-1x−y−1 个原地不动的命令，总命令数为 2x−12x-12x−1。#include <bi

2020-11-20 18:08:08 198

原创 POJ 3889 分形

题意传送门 POJ 3889题解NNN 级城市由 444 个 N−1N-1N−1 级城市构成，总共有 4N4^N4N 个房屋，左上角/左下角为 N−1N-1N−1 级城市顺时针/逆时针旋转 90°90\degree90°，且序号呈水平翻转；右上角/右下角为 N−1N-1N−1 级城市，且序号不变。为了递归求解同一个级别 NNN 的城市，在判断目标房屋在左上/左下/右上/右下哪一个 N−1N-1N−1 级城市中后，需要做相应的坐标变换。#include <algorithm>#inclu

2020-11-20 17:38:54 171

原创 POJ 1740 博弈

题意传送门 POJ 1740 A New Stone Game题解博弈问题，考虑 grundygrundygrundy 数，但每一轮除了选择一堆石子移除至少一颗石子外，其余石子还可以任意转移到其他石子数非零的堆中，状态转移数过多。需要寻找其他方法判断胜负态。在这类游戏当中，做出对称的状态再完全模仿对手的策略常常是有效的。取走最后一堆石子的一方获胜，那么此时必胜态为仅剩一堆石子的状态。考虑 NimNimNim 的异或运算构造必胜态，可以发现甚至无法通过样例。考虑石子堆的奇偶性，当石子堆中石子出现

2020-11-20 15:50:05 178

原创 POJ 3688 博弈 + DP

题意传送门 POJ 3688题解博弈问题，考虑胜负态转移；一般的必胜态只要判断是否能转移到必败态即可，但这个游戏必须要保证必胜态可转移到的非重赛的状态都是必败态。在这个游戏中，存在 444 种状态：必胜态、必败态、胜负都可能的状态、需要重赛的状态。需要求解的是必胜态的个数。可以发现先手或后手获胜的情况只与牌上数字的组合有关，与出牌顺序无关。考虑到先手石子堆没有石子时为必败态，那么有：先手必胜时，石子堆的个数只能由奇数个牌上的数字组成。dp[i][j][k]dp[i][j][k]dp[i][j][k

2020-11-19 11:57:17 188

原创 POJ 1845 逆元 / 分治

题意传送门 POJ 1845题解逆元分解质因数 A=p1e1p2e2…pnenA=p_1^{e_1}p_2^{e_2}\dots p_n^{e_n}A=p1e1p2e2…pnen，则约数和为(1+p1+⋯+p1e1)(1+p2+⋯+p2e2)…(1+pn+⋯+pnen)(1+p_1+\dots +p_1^{e_1})(1+p_2+\dots +p_2^{e_2})\dots (1+p_n+\dots +p_n^{e_n})(1+p1+⋯+p1e1)(1+p2+⋯+p2e2)

2020-11-17 00:01:20 157

原创 HDU 6390 容斥原理

题意传送门 HDU 6390题解若 gcd(a,b)=1gcd(a,b)=1gcd(a,b)=1，则有 ϕ(a)ϕ(b)=ϕ(ab),Gu(a,b)=1\phi(a)\phi(b)=\phi(ab),G_u(a,b)=1ϕ(a)ϕ(b)=ϕ(ab),Gu(a,b)=1 假设 a=pe1,b=pe2a=p^{e1},b=p^{e2}a=pe1,b=pe2，则有 Gu(a,b)=pe1+e2∗(p−1)/p[pe1∗(p−1)/p]∗[pe2∗(p−1)/p]=pp−1G_u(a,b)=\frac{p^

2020-11-16 12:26:48 157

原创 Codeforces Round #682 (Div. 2)

传送门Codeforces Round #682 (Div. 2)A显然元素全部为 111 的数组满足要求。#include <algorithm>#include <cmath>#include <cstdio>#include <cstring>using namespace std;int main(){ int t, n; scanf("%d", &t); while (t--) {

2020-11-14 23:34:23 138

原创 P3131 [USACO16JAN] 前缀和

题意传送门 P3131 [USACO16JAN]Subsequences Summing to Sevens题解求区间和模 777 为 000 的最大区间长度。计算前缀和，记录模 777 的余数在前缀和中第一次与最后一次出现的次数，此时满足前缀之差代表的区间和模 777 等于 000 且对于此余数而言区间长度最大。#include <bits/stdc++.h>using namespace std;#define mod 7#define maxn 50005int N, fs

2020-11-14 15:51:55 172

原创 P3128 [USACO15DEC] LCA + 差分

题意传送门 P3128 [USACO15DEC]Max Flow题解求树上被运输路线覆盖次数最多的点。考虑树上差分，倍增求 LCALCALCA，设路径的两个端点分别为 u,vu,vu,v，差分数组为 ddd，父节点数组为 parparpar，则对于每一条路径du←du+1d_u\leftarrow d_u+1du←du+1 dv←dv+1d_v\leftarrow d_v+1dv←dv+1 dLCA←dLCA−1d_{LCA}\leftarrow d_{LCA}-1dLCA←dLCA−1

2020-11-14 15:37:40 131

原创 LeeCode 514 DP

题意传送门 LeeCode 514. 自由之路题解两串的匹配问题，考虑 DPDPDP。匹配过程中 keykeykey 被匹配的长度是单调递增的，保证 DPDPDP 的第一维是一个 DAGDAGDAG；考虑到不能匹配当前 keykeykey 当前字符的环上编码可以顺时针、逆时针转移，那么为了保证 DPDPDP 的状态转移是一个 DAGDAGDAG，只考虑可以匹配 keykeykey 当前位置的环上编码，可以预处理出 pos[′a′−′z′]pos['a'-'z']pos[′a′−′z′]。dp[i][j

2020-11-13 10:09:40 69

原创 POJ 2068 博弈 + DP

题意传送门 POJ 2068题解轮到某队队员时剩余 000 块石子的情况是必胜态。这是个石子递减的游戏，从小到大枚举石子数量，进行胜负态转移。dp[i][j]dp[i][j]dp[i][j] 代表剩余 iii 块石子时轮到第 jjj 名队员的胜负态dp[i][j] ∣= !dp[i−k][(j+1) mod 2n],1≤k≤M[i]且i−k≥0dp[i][j] \ |=\ !dp[i-k][(j+1)\ mod \ 2n],1\leq k\leq M[i] 且

2020-11-12 23:21:47 108

原创 P4869 线性基

题意传送门 P4869 albus就是要第一个出场题解求某个数在集合子集异或值中按序出现的下标。设全集为 SSS，线性基 BBB，则线性基非空子集异或的值域有 2∣B∣−12^{|B|}-12∣B∣−1 个不同且非零的数，S−BS-BS−B 的元素都唯一存在与线性基异或为零的组合，那么其子集也满足这个性质。则 SSS 的子集非零异或值域的任一个数字都有 2∣S∣−∣B∣2^{|S|-|B|}2∣S∣−∣B∣ 个，则总个数为 (2∣B∣−1)2∣S∣−∣B∣(2^{|B|}-1)2^{|S|-|B|}(

2020-11-12 21:40:56 118

原创 P4151 [WC2011] 线性基 + DFS

题意传送门 P4151 [WC2011]最大XOR和路径题解异或运算和与运算、或运算不同，没有单调性，难以使用一般求最短路的方法求解此问题。考虑环，显然图中环的异或值都是可以取到的。考虑两种情况，环不在路径上，那么可以反复经过某一条通向环的路径，且此路径没有贡献；环在路径上，那么将换上的原路径部分替换为环上非原路径的部分。那么 DFSDFSDFS 求图中所有简单环的异或值，构造线性基。任取一条 1→N1\rightarrow N1→N 的路径，从高位向低位不断尝试异或这一位的线性基，求最大值即可。

2020-11-12 10:14:23 114

原创 HDU 3949 线性基

题意传送门 HDU 3949题解处理集合异或的值域，考虑线性基。线性基的子集的异或即集合的异或非零值域，若集合中出现可以由线性基表达的值，或说线性基的规模小于集合规模，那么集合异或的值域包含零，需要特殊处理。考虑值域第 kkk 大，从高位向低位用类似二分的方法处理。假设已处理好大于 iii 位的线性基，若第 iii 位存在线性基，[0,i−1][0,i-1][0,i−1] 位存在 nnn 个线性基，那么可以构成 2n2^n2n 个包含空集且小于包含第 iii 位线性基的值，根据 kkk 值判断是否使

2020-11-11 23:10:29 146

原创 P3293 [SCOI2016] 贪心 + 主席树

题意传送门 P3293 [SCOI2016]美味题解求某个数字 xxx 与集合中任一数字异或的最大值，考虑从高位向低位贪心的选择与 xxx 在这一位的异或值对应的数。假设当前已处理好比第 iii 位更高的数位，此时区间左界为 lblblb，假如 xxx 的第 iii 位为 111，那么选择的数字应该在区间 [lb+(1<<i),lb+(1<<(i+1)))\big[lb+(1<<i), lb+(1<<(i+1))\big)[lb+(1<<i)

2020-11-11 10:31:05 84

原创 SGU 275 线性基

题意传送门 SGU 275 to xor or not xor题解选择线性基的一组子集，使异或值最大。构造一组线性基，使之第 iii 位为 111 的值至多出现在一个元素中，则线性基的异或值即答案。#include <algorithm>#include <cmath>#include <cstdio>#include <cstring>using namespace std;#define maxl 60typedef long long

2020-11-10 20:32:15 331

原创计蒜客 A1957 DP + 费马小定理

题意传送门计蒜客 A1957 [ACM-ICPC 2018 南京赛区网络预赛] The Great Nim Game题解求 NNN 堆石子任选部分堆进行 NimNimNim 游戏，使各堆石子数量异或非零的选取方案个数。对于异或，有如下性质x⊕x=0,x⊕y≠0(x≠y)x\oplus x=0,x\oplus y \neq 0(x\neq y)x⊕x=0,x⊕y=0(x=y) 考虑按照堆数为奇偶分别处理石子数量相同的堆，dp[i][j]dp[i][j]dp[i][j] 代表考虑前 iii 种

2020-11-10 18:34:35 106

原创 LeeCode 327 坐标离散化 + BIT

题意传送门 LeeCode 327. 区间和的个数题解考虑以 iii 为满足条件的区间右界，设 iii 的前缀和为 preSum[i]preSum[i]preSum[i]，对于每一个右界，目标是统计满足preSum[i]−preSum[j−1]∈[lower,upper]preSum[i]-preSum[j-1]\in [lower, upper]preSum[i]−preSum[j−1]∈[lower,upper] 的左界 jjj，那么坐标离散化前缀和，从前向后枚举右界，用树状数组维护 j,j&l

2020-11-09 22:25:29 308

原创 Codeforces 138D Nim + SG

题意传送门 Codeforces 138D题解发生分割的游戏，考虑 GrundyGrundyGrundy 值。与直接分割平面不同之处在于游戏场地是方格状的，考虑将平面上的格子沿主对角线交替染为黑（坐标和为奇数）、白色（坐标和为偶数），黑色与白色格子切割时是互不影响的。那么将游戏按黑白格子考虑两个独立的分割游戏，用 NimNimNim 游戏的方法判断胜负。为了方便处理，将矩形旋转 454545 度，即同一条主对角线的点横坐标相同，同一条副对角线上的点纵坐标相同，然后将最小值设为零，即 (x,y)(x,

2020-11-08 22:00:43 106

原创 POJ 3537 Nim + SG

题意传送门 POJ 3537题解在博弈中，当前状态的 GrundyGrundyGrundy 值可以看做 NimNimNim 中的一个石子堆，就可以按照 NimNimNim 用异或的方法求解。当前状态的 GrundyGrundyGrundy 值等于除任意一步所能转移到的状态的 GrundyGrundyGrundy 值以外的最小非负整数。此游戏中率先进入 333 个连续的×状态的玩家获胜。设当前游戏规模为 nnn，先手玩家选择位置 iii 为第一步，那么，若后手玩家选择 [i−2,i+2][i-2

2020-11-07 10:33:14 80

原创 POJ 2975 Nim

题意传送门 POJ 2975题解NimNimNim 游戏的必败态是 k1⊕k2⋯⊕kn=0k_1\oplus k_2\dots \oplus k_n=0k1⊕k2⋯⊕kn=0。假设当前为必胜态，有 k1⊕k2⋯⊕kn=xk_1\oplus k_2\dots \oplus k_n=xk1⊕k2⋯⊕kn=x，那么转移的目标是使异或值为零。每一步只能从某一堆石子中取大于等于 111 的石子，假设取第 iii 堆，那么目标是使 ki′=ki⊕xk_{i}'=k_i\oplus xki′=ki⊕

2020-11-06 23:52:25 131

原创 POJ 1082 博弈 + DP

题意传送门 POJ 1082题解设日历上 (y,m,d)(y,m,d)(y,m,d) 的下一个日期为 nxt1(y,m,d)nxt1(y,m, d)nxt1(y,m,d)，下一个月的同一个日子的日期为 nxt2(y,m,d)nxt2(y,m,d)nxt2(y,m,d)。dp[(y,m,d)]dp[(y,m,d)]dp[(y,m,d)] 代表日期 (y,m,d)(y,m,d)(y,m,d) 亚当先手是否可以取胜dp[(y,m,d)] ∣= !dp[nxt1(y,m,d)]&nbsp

2020-11-06 23:12:14 102

原创 POJ 2409 Pólya 计数定理

题意传送门 POJ 2409题解要求将手镯旋转、翻转看做相同的状态，考虑 PoˊlyaPólyaPoˊlya 计数定理。设颜色有 mmm 种。对于旋转，有旋转 [0,n)[0,n)[0,n) 个 2π/n2\pi /n2π/n 的角度的置换，即 nnn 个置换群。考虑旋转 kkk 个位置后和原来相同的染色方案个数，将颜色与第 iii 个石头相同的石头的集合叫做 iii 的轨迹，个数为 n/gcd(k,n)n/gcd(k,n)n/gcd(k,n)。因此，旋转 kkk 个位置后和原来相同的染色方案数就是

2020-11-06 21:38:05 141

原创 POJ 1286 Pólya 计数定理

题意传送门 POJ 1286题解组合问题中，要求把旋转和翻转之后相同的状态看成是本质相同的状态，考虑 PoˊlyaPólyaPoˊlya 计数定理。设颜色有 mmm 种。对于旋转，有旋转 [0,n)[0,n)[0,n) 个 2π/n2\pi /n2π/n 的角度的置换，即 nnn 个置换群。考虑旋转 kkk 个位置后和原来相同的染色方案个数，将颜色与第 iii 个石头相同的石头的集合叫做 iii 的轨迹，个数为 n/gcd(k,n)n/gcd(k,n)n/gcd(k,n)。因此，旋转 kkk 个位置后

2020-11-06 21:11:11 169

原创 P2257 莫比乌斯反演 + 整除分块

题意传送门 P2257 YY的GCD题解求 ∑pmin(n,m)∑i=1n∑j=1m[gcd(i,j)=p]\sum\limits_{p}^{min(n,m)}\sum\limits_{i=1}^{n}\sum\limits_{j=1}^{m}[gcd(i,j)=p]p∑min(n,m)i=1∑nj=1∑m[gcd(i,j)=p]，不易直接求解，考虑 d∣gcd(i,j)d|gcd(i,j)d∣gcd(i,j) 的关系，设F(d)=∑d∣kf(k)F(d)=\sum\limits_{d|k}f

2020-11-06 15:33:03 104

空空如也

空空如也