雯舞

这个题怎么Difficulty rating 80%80\% 啊，送经验的感觉啊D(m,n)=======∑d|m∑k=1nσ0(kd)∑d|m∑k=1n∑a|k∑b|d[(a,b)=1]∑a=1n⌊na⌋∑d|m∑b|d[(a,b)=1]∑a=1n⌊na⌋∑d|m∑b|d∑i|a,i|bμ(i)∑i|mμ(i)×(σ0∗1)(mi)∑id≤n⌊nid⌋∑i|mμ(i)×(σ0∗1)(mi

(k+1)|(n+k)(k+1) | (n+k) 就是 (k+1)|(n−1)(k+1) | (n-1) 这个函数就是最大的 kk 使得 1,2,⋯,k1,2,\cdots,k 都整除n-1吧 随便容斥下咯#include#include#includeusing namespace std;typedef long long ll;inline ll Gcd(ll a,ll

直接考虑每一位的贡献，应该是一个 ai×i×(1+10+⋯+10n−i)a_i\times i\times (1+10+\cdots+10^{n-i}) 的形式，这就是个等比数列求和 然后因为是循环串，还是个等比数列求和，就好了#include#include#includeusing namespace std;typedef long long ll;const int P

二分转化为μ2\mu^2的前缀和。 然后转化为O(n13)O(n^{1\over 3})的运算，但是需要预处理μ\mu的前缀和，大力杜教筛求和。 然后感谢阿爷教我把二分改成了迭代，小范围内一个一个挪，用rho求μ(n)\mu(n)，然后就能卡进去了。 复杂度似乎是萎的吧。#include#include#include#include#include#include#defi

整体二分，然后问题变成，子集加，单点查询，然后像CTSC吉夫特 可以用经典的二进制分高位低位的搞搞。调个参，大概是高5位低12位。 不知道在线怎么做。#include#include#include#include#define pb push_backusing namespace std;typedef long long ll;inline char nc(){

建出后缀树，记录每个子串的出现次数，然后二分下答案在哪个子串中就好了退役选手不会写后缀自动机#include#include#include#include#define pb push_backusing namespace std;typedef long long ll;inline char nc(){ static char buf[100000],*p1=b

原标题 高三志开坑 2017.09.02 23:47 高三，还是记录一下吧。 DATE EVENT HINT 07.22 退役 07.25-07.27 CodeM Final，Beijing 点赞 08.04-08.18 Codechef August Long Challenge 2017 高三可能都不会这么认真的打比...

抱歉，退役了，Blog下的留言我可能不会回答了吧。UPD. 博客已搬 [雯舞 - Manchery's Space](manchery.co) ——2018.05.25

跟 51Nod 1356 代数数的次数 是一样的 不过这里都是质数 也就是就是 2n2^n 关键是输方案 这个不一定有二次剩余感谢sxt 一个一个数加进答案 转化成 已知F(x)F(x)，求F(x+a√)F(x+\sqrt a)和F(x−a√)F(x-\sqrt a)的系数 这个推一下就是一个FFT拷了myy的板子 自己的太慢了QAQ#include<cstdio>#include<cs

大意是把二元运算化成一元运算 感谢sxt的帮助剩下的看代码注释吧 //我觉得蛮清楚了 唯一不爽的是 没想到临时变量这么多 mmp#include<cstdio>#include<cstdlib>#include<algorithm>#include<vector>#define pb push_backusing namespace std;int B;struct abcd{ i

随便搞搞吧#include<cstdio>#include<cstdlib>#include<algorithm>using namespace std;typedef long long ll;#define read(x) scanf("%d",&(x))const int maxn=200005;const int P=1e9+7;const int N=505;int n,a[N]

这是若干个 2xai+12x^{a_i}+1 的东西的卷积 然后这个FWT一下发现每一项只有 −1-1 或 33 那么卷积的FWT每一项就是若干个 −1-1 和 33 的乘积 这个不好求 直接加在一起FWT，那么我们得到了每一项 −1-1 和 33 的和 因为只有这两个取值，可以直接解方程有多少个 −1-1 和 33 然后快速幂再乘起来#include<cstdio>#include<

fi,jf_{i,j}表示长为 ii 的区间 最大值是jj 的答案 转移就枚举最左边的最大值在区间的位置 前缀和优化下就好了好像也可以fi,j,kf_{i,j,k}表示前 ii 个，末尾 KK 个中最大值在 jj ，最大值是 kk，有点复杂？#include<cstdio>#include<cstdlib>#include<algorithm>using namespace std;typ

一个图的 kk 染色数是关于 kk 的 nn 次(?) 多项式 称为色多项式 那么这里模6 我们只要知道模2和模3的值 然后分类讨论下就好了 一张图的0染色数是0，1染色数等于[m=0][m=0]，2染色数与二分图的联通块个数有关#include<cstdio>#include<cstdlib>#include<algorithm>#include<cstring>#define

未经允许搬了过来代码略丑#include<cstdio>#include<cstdlib>#include<algorithm>using namespace std;typedef long long ll;inline char nc(){ static char buf[100000],*p1=buf,*p2=buf; return p1==p2&&(p2=(p1=buf)+f

请教了比利 首先 2n2^n 次项是没有问题的，比如2√+3√\sqrt2+\sqrt3，可以构造f(x)=(x+2√+3√)(x+2√−3√)(x−2√+3√)(x−2√−3√)f(x)=(x+\sqrt2+\sqrt3)(x+\sqrt2-\sqrt3)(x-\sqrt2+\sqrt3)(x-\sqrt2-\sqrt3) 这个不是最优解 最优解应该是构造这样一个矩阵 每一行代表一个质因子

本来想自己再推一下的，但是退役了也就弃坑了 未经允许的搬了搬题人的题解// BEGIN CUT HERE #include<conio.h>#include<sstream>// END CUT HERE #include<cstdio>#include<cstdlib>#include<cstring>#include<cmath>#include<iostream>#in

最小值的最大值小于等于最大值的最小值 那么只要有一行一列最小值等于最大值，那么矩阵就合法 如果我们枚举这一行和这一列，那么最小值最大值一定就是他们的交点 那么再枚举交点最后的值是多少 这样复杂度是O(n4)O(n^4) 换一下 先枚举值 再枚举位置 可以去一个代价最小值 O(n2logn)O(n^2\log n)#include<cstdio>#include<cstdlib>#in

这种非简单路 一般转化成任意一条路加上若干环 这里大概是任意一条路加上若干环长的gcd 任意一条路 可以弄出任意一颗生成树？ 但是这里实际上只需要一个带权并查集 详见官方题解#include<cstdio>#include<cstdlib>#include<algorithm>using namespace std;typedef long long ll;typedef pair<

μ2(n)=∑d2|nμ(d)\mu^2(n)=\sum_{d^2|n} \mu(d) 然后就是xjb推 反正退役了 我也就弃坑了 95分代码 复杂度分析及优化详见官方题解#include<cstdio>#include<cstdlib>#include<cmath>#include<algorithm>using namespace std;typedef long long l