2017年08月_High_EnergyElectron

原创解题报告：HDU_6123 Destroy the cube (容斥+三元环计数)

题目链接题意及官方题解：记录一下自己的做法：首先如果可以直接跑全部的黑色位置，那就很好写，但是肯定会超时，所以一定要用对称性优化。如果n为偶数，可以完美拆成八个完全一样的小正方体的子问题，很好写如果n为奇数，问题就变得复杂而且要考虑各种细节。。先不考虑最中间的面，解决八个小正方体的子问题把中间面的所有黑色位置建图跑三元环然后发现会有重

2017-08-30 20:25:12 807

原创解题报告：HDU_6129 Just do it (找规律两种做法)

题目链接题意及官方题解：思路：看到另一种做法，要巧妙一点，记录一下解法一（官方）：打出当前位对后面位的贡献表，发现是个斜杨辉三角只有组合数为奇数才用贡献，由Lucas可知组合数C(n,m)为奇数等价于(n&m)==m 这样就可以枚举m(1~n-1)快速更新答案虽然复杂度看上去是O(n^2)，但是满足要求的组合数并不多（不会证明。。），能

2017-08-21 21:24:32 320

原创解题报告：HDU_6122 Color the chessboard (计数)

题目链接题意及官方题解：思路：分析题意可以发现计数只需要维护2*2的矩阵满足题目要求即可将奇数格的颜色翻转发现矩阵只会有三种形式： 1、每一行颜色相同 2、每一列颜色相同 3、全部的颜色相同然后容斥一下即可代码： #include const int mod = 998244353; const int N = 1e3+5; using n

2017-08-21 17:25:41 364

原创解题报告：HDU_6139 Galaxy at War (阶梯博弈)

题目链接题意：一张n*m的表格上有一些格子有一些水晶球，两个人轮流进行游戏每次选择一个有水晶球的格子，选择其中至少一个水晶球将它左移或者下移，不能出界还有一些格子上有M(Meditations)或者P(pollutant sources)，对应的作用为当你选择的格子上有M时，若你选择移动 t 个水晶球，那么会将2*t 个水晶球平分到可以移动到的格子内当t个水晶被移动到有

2017-08-20 21:06:49 426

原创解题报告：HDU_6134:Battlestation Operational (莫比乌斯反演)

题目链接题意：求思路：本来出题人想考的不是反演，但是用反演做意外的简单。。原式：做反演：令易知：（ D( x ) 为 x 的因子个数）那么可在内预处理出 g() ，再在线性时间内得到g()的前缀和每次查询的复杂度，总复杂度代码： #include const int N = 1e6+10; const l

2017-08-18 17:06:13 259

原创解题报告：HDU_6128：Inverse of sum (二次剩余)

题目链接题意及官方题解：思路：已知公式：转换一下：对于每个y，满足要求的x为：只需要求的y的系数在(%mod)意义下的等价式即可等同于求sqrt(-3)的等价式，也就是求p-3在(%mod)意义下的二次剩余( Cipolla's algorithm ) 注意一些细节：不考虑0、mod=2、mod=3、x==y时的计数以及不存在p-3

2017-08-18 16:39:43 345

原创解题报告：HDU_6127:Hard challenge (极角排序)

题目链接题意及官方题解：补充：我是以到x负半轴的弧度进行排序，然后扫过（0~PI）的弧度，中间每扫过一个点都要可能更新答案代码： #include const double pi=acos(-1.0); using namespace std; class point { public: int x,y,v; double radi

2017-08-18 10:33:52 241

原创解题报告：HDU_6136：Death Podracing (优先队列+循环链表)

题目链接题意： n个人以不同的速度在环上顺时针或逆时针移动，每次相遇，移除下标小的，问最后只剩下一个人的时间的分数形式官方题解及思路：也不是第一次写循环链表的题了，还是写了好久。。注意维护循环链表时要同时更新左右指针代码： #include const int N = 1e5+10; using namespace std; struct no

2017-08-18 10:15:28 512

原创解题报告：CROC 2016 - Elimination Round F. Cowslip Collections (莫比乌斯反演)

题目链接题意：给定一个大小为n(n 每次询问给往集合中插入一个数，求当前集合的贡献集合的贡献为所有大小为k的集合的gcd之和思路：刚开始想歪了令 f[ x ]：集合大小为k且gcd()==x的集合数量 num[ x ]：所有的数中x倍数的个数那么答案为刚开始想歪了：想用树状数组维护后面的部分，均摊复杂度应该能过。。但是维护树

2017-08-14 20:04:57 340

原创解题报告：Codeforces Round #325(Div. 1) E. Present for Vitalik the Philatelist (莫比乌斯反演)

题目链接题意：给n(n 能从中选出一个数x，在从剩下的数中选出一个非空集合S，满足： S的gcd不为1，gcd(S,x)为1 询问满足的方案数闲扯：本来O（nlog(n)）的想法，交上去一直T。。。以为卡常数，优化了还是T。。。发现中间有个部分写成O(n*sqrt(Ai))的级别。。。改了后3000+ms过了。。。思路：

2017-08-13 21:13:18 316

原创解题报告：Codeforces Round #142 (Div. 1) C. Triangles (思维好题)

题目链接题意：给定一个图，求这个图以及它的补图中长度为3的环的数目思路：一道不错的思维题为了方便表述，定义原图给出的边为实边，补图的边为虚边那么所有三角环的组成一共有4种格式为了制定出满足题意的规则，设： a：一条虚边+一条实边的点数 b：两条虚边的点数 c：两条实边的点数根据4种格式的三角环可以得到联立方程组：对应的意思为：

2017-08-13 20:16:58 315

原创解题报告： Codeforces Round #305 (Div. 2) E.Mike and Foam (莫比乌斯反演)

题目链接题意：有一个长度为n(n 序列中的每个位置对应一个数ai(1 每次操作给一个数x，代表挑出的序列中的下标为x的数若对应下标的数已经被挑出，那么就把它放回去每次操作完，询问所有挑出的数中互质的对数刚开始时所有数都没有被挑出思路：刚开始的没有数被挑出，那么初始答案为0。令num[i] ：挑出的数中 i 的个数如果当前操作的数为x，

2017-08-12 22:06:58 238

原创解题报告：Codeforces Round #146 (Div. 1) E. Number Challenge 莫比乌斯反演

题目链接题意：给定n,m,l，求 d(x)：x的约数个数思路：加强版的约数个数和（解题报告：BZOJ_3994 约数个数和莫比乌斯反演学习题）？代入公式得到得到代入常用的公式：得到：这个式子的复杂度为但是可以发现对于一个数k，只需要用到和它互质的部分那么可以在的时间内预测出每个数的互质的所

2017-08-12 20:13:01 513

原创解题报告：HDU_6102 GCDispower (离线处理+莫比乌斯反演)

题目链接题意：给一个全排列，每次询问一个区间的函数官方题解：思路：题解已经说的很详细了Orz....记录一下容斥反演的部分令num(d)为已经扫过的数中d的倍数的个数，预处理出每个数的有效因子（对应的反演系数不为0），就能在均摊O(log(n))的时间得出那么每个当前正在扫的数对它前面的区间的贡献为 A[r]为当前的右端

2017-08-11 11:06:50 540

原创解题报告：BZOJ_3529：数表莫比乌斯反演+离线处理

题目链接题意：一张n*m的表格，( i , j )的格子上的数为gcd( i , j ) 的约数和，询问所有格子的不超过a的数之和，T组数据。 1 思路：这题很多细节不是很容易想到。。为了方便表述，定义： F(x)：x的约数和，可以预处理得到 G(x)：满足题目的gcd==x的组数我们先分析g[i]：得到：这个看上去已经很

2017-08-07 21:25:45 310

原创解题报告：HDU_4947:GCD Array 莫比乌斯反演+树状数组

题目链接题意：长度为n的数组，初始值为0。有两种操作： ①给定a,b,x，将所有下标 i 满足 gcd(a,i)==b的位置加上x ②给定x，查询 1~x 的前缀和思路：对于操作①，显然无法更新所有的满足要求的下标稍作分析可以发现，对于每次更新操作的贡献，只需要更新每个反演式的第一个位置那么查询的时候，每个位置的贡献就是( x/i ) 因为

2017-08-07 15:20:40 332

原创解题报告：HDU_5909:Tree Cutting 树上FWT

题目链接题意：给定一棵树，树上每个结点有它的权值，现在问在树上选一些连通的点，使得这些点的权值的异或结果等于0~m-1的时对应的方案数。思路：对每个子树统计选定根的方案数，那么答案为：因为每个分支可以选也可以不选，那么在进行FWT之前对每个儿子进行dp[son][0]++的操作代码： #include const int mod = 1e9

2017-08-03 20:53:16 224

原创解题报告：CodeForces - 662C:Binary Table FWT(快速沃尔什变换)

题目链接题意：给定一个至多20*100000的01矩阵，可以选择任意行或列，选择行或列的01值改变，问经过操作能到最少的含1的数量。思路：因为行最大为20，考虑将每一列压缩成一个20位的整数，选择的行也最多只有20，同样压缩到一个20位的整数。改变行上的01值，对应的操作为异或，那么就能写出对应的卷积式：根据异或运算的性质，可以交换j,k的位置，得到： C[k

2017-08-03 08:50:21 802

原创解题报告：HDU_6061 RXD and functions NTT

题目链接题意及官方题解：思路：先感谢Q巨指点Qrz... 先求得系数式：拆开组合数：把系数提取出来稍作变换：整理一下：得到：然后就可以直接进行NTT了。。代码； #include const int mod = 998244353; const int N = 4e5+10; con

2017-08-02 10:27:38 1562 2

原创解题报告：BZOJ_2154: Crash的数字表格莫比乌斯反演

题目链接题意：给定n,m( 思路：预处理和的前缀和，这样就能在时间内得出结果总复杂度为代码： #include using namespace std; const int mod = 20101009; const int N = 1e7+10; vectorpr; bool Np[N]; int mu[N]; in

2017-08-01 08:49:22 365

Hee的博客~