The Blog of Creed_

摘要： https://loj.ac/problem/6062 直接套用霍尔定理。 由于A有多个选择，考虑维护B是否合法。 首先B数组的顺序显然是没有用的，可以直接排序。 然后每个A就都变成了向一个后缀连边。 对于B，原本需要check每一个集合是否满足|u|=0 因此，只需要维护F【i】 的最小值即可。 阅读全文

摘要： 先是几个比较简单的定理 最大匹配=最小点覆盖=n 最小边覆盖=n 最大独立集 最小点覆盖和最大独立集互补 霍尔定理 对于任意集合a，都满足|a| 阅读全文

摘要： 非常神仙的一道题！ 题意：给出某n个数字跑完全背包m容量的dp数组，求满足要求的字典序最小的n个元素，不知道n是多少。 首先考虑付公主的背包这个题。 对dp数组求一个ln，设它为F。 已知 e^(G1+G2+G3)=e^F，其中Gi是第i个物品的生成函数求ln。（重量为i的物品的Gi=∑ 1/i ✖ 阅读全文

摘要： https://loj.ac/problem/6072 虽然结合了很多算法，但是一步一步地推一下还不算太难的一道题。 首先考虑枚举枚举有用的苹果的集合，然后去算生成树个数。 先考虑怎么计算生成树个数。 发现可以使用matrix tree。 所有有用点可以和有用点以及坏点连边，所有不是坏点的无用点只能 阅读全文

摘要： 好像没什么好说的。 生成树方案=abs【（度数矩阵-邻接矩阵）删掉任意一行一列】 设A为一个矩阵。 abs(A)等于对它跑高斯消元然后把对角线上的数字乘起来。 顺便提一下高斯消元的辗转相除法。（用于模数不为质数的情况） 在用第i行对第j行消元的时候，设A为a[i][i]，B为a[j][i]。 然后用 阅读全文

该文被密码保护。 阅读全文

摘要： 做法挺多的一道题 第一种做法就是直接考虑min max容斥。 然后可以转成一个n/logn sqrt(n)的dp。 不多做赘述。 第二种，ans=sigema x p(len =x)。 令f(i)=sigema x p（前x个数的gcd是i）。 则ans=sigema i=2 n f(i) 考虑对f 阅读全文

摘要： 终于找到了一个只会用[gcd(i,j)==1] = sigema d|gcd(i,j) mu(d) 做不了的题。 考虑枚举gcd后。 此时，ans可以表示为一个 sigema x f(x)的形式。 考虑对反演f(x)。 然后发现f(x)也很容易在nlogn的复杂度内算出来，就做完了。 cpp inc 阅读全文

摘要： 又一道没有做出来的题 首先做最基础的mobius反演，可以得到，答案就是一堆fi^k乘起来，k是下面这个东西。 套路的用T代替id，并枚举T。 指数的式子变成 π(i|T) mu（T/i）✖（n/T）✖（m/T） 考虑先枚举指数再枚举i（没想到这一步！） 令dp[T]=π (i|T) f[i]^mu 阅读全文

摘要： https://www.luogu.org/problemnew/show/P3312 比较nb的一个题。 首先考虑没有a的限制怎么搞。 把d提到前面 mobius反演一下 令T=dx,把T挪到前面（没想到这一步！） 设F(x)=d(x)和mu(x)的狄利克雷卷积。 然后考虑怎么带上a的限制，只要能 阅读全文