2016年07月_Danliwoo

原创 HDU 5768 Lucky7 数论中国剩余定理

原题见HDU 5768求[l,r]范围内是7的倍数，同时不满足任意一个给定的同余式的数的个数。如范围为[1,100]，不满足模3余2或模5余3的7的倍数有7,21,42,49,70,84,91 ，故答案为7. 其中除数都为非7的素数(≤105)(\le 10^5)，除数的乘积小于101810^{18}.同余式最多有15个。分析所有[l,r]范围内7的倍数减去满足任意一个同余式的数，即得到了答案。而

2016-07-28 19:03:33 2562 3

原创 CF 429B. Working out DP

原题见CF 429B给一个n*m的方格，A从左上角走到右下角，B从左下角走到右上角，路线交叉处的权值不算，问两条路线权值之和最大值。要求：两条路线只在一点交叉。分析可以枚举交叉点，求该点到四个角落的权值之和。到每个角的权值都可以dp得到最大值。从左上角顺时针得到0,1,2,3四个方向。一开始没注意到题意里的只在一点交叉，妥妥跪了。 CF题解里的图片灰常清晰两幅图的权值之和分别是： dp[

2016-07-26 10:43:04 1330

原创求逆元的几种办法

补了一道CF题，顺便看到有人求逆元的方法叼叼的，怪不得人家过题这么快啊…… 来总结一发……一般求法求a关于N的逆元，即要解同余方程ax≡1(modN)ax\equiv 1(\mod N)的解x. ax≡1(modN)⇔ax+Ny=1ax\equiv 1(\mod N)\Leftrightarrow ax+Ny=1 仅当a与N互质时，存在aa的逆元，利用扩展欧几里得求解。这里N不一定是素数L

2016-07-24 20:14:36 1642

原创 HDU 5728 PowMod 欧拉函数递归

感觉智商被掏空…原题见HDU 5728定义k=∑mi=1φ(i·n) mod 1000000007k=\sum_{i=1}^{m} \varphi (i·n)\ mod\ 1000000007 nn是无质因子平方项的数. 求ans=kkkk...k( mod p)ans=k^{k^{k^{k^{...^k}}}}(\ mod \ p)，其中k有无穷多个数据范围：1≤n,m,p≤1071 \l

2016-07-22 22:24:05 6248

原创 HDU 5135 Little Zu Chongzhi's Triangles 状压

原题见HDU 5135给n（n<13）条边，每条边只能用一次，拼成多个三角形，求三角形面积和的最大值。分析1~2122^{12}只有4096，对每个数二进制分解，第i位为1则用aia_i这条边。预处理出所有的数位和为3的倍数的边。 d[i] i的二进制数位和 f[i] 用了i的二进制所表示的边后，可以组成的三角形面积的最大值。如果d[i] = d[j]+3且i|j <= i（避免0被1与

2016-07-21 22:13:13 630

原创 HDU 5726 GCD 区间gcd查询 MAP RMQ 优化

原题见HDU 5726给N(N≤100000N\le 100000)个数，Q(Q≤100000Q\le 100000)个询问，每次查询输出区间的最大公约数，以及最大公约数为这个数的区间数目。分析查询次数很多，要做预处理，用map<最大公约数,区间数>存下来，实现O(1)的查询。预处理发现对于同一左端点的区间而言，右端点越靠右，区间gcd单调递减。因此可以固定左端点，二分右端点，找到gcd突变的右

2016-07-21 20:26:29 1134

原创 HDU 5724 Chess 从懵逼到学会 SG函数

引入经典的NIM博弈给若干堆石子，Alice和Bob轮流取石子，每次可选择其中的一堆拿走任意多的石子，但不能不拿。最后没有石子可拿的人输。已知：堆数N，每堆石子数目aia_i，Alice先取。求：谁是胜者。博弈基于足够聪明的两者，每次都尽量取到必胜的状态，绝不存在其中一人明明可以胜出却要让别人赢的情况。因此输赢往往由局面本身确定。即存在必胜态、必败态。上面这个游戏的必败态，即玩家面对的是没有

2016-07-20 15:26:33 1193

原创好好学习啊少年

成功了都是经验，失败了都是教训。难道失败者就没有做得对的地方了吗？所以….还是要好好学习啊….今天好困，没刷题，没文章….整个一辣鸡

2016-07-12 19:48:52 731

原创莫比乌斯教我变魔术

莫比乌斯反演详细的不赘述了，都在这里若f(n)f(n)是积性函数，则如下F(x)F(x)也是积性函数F(x)=∑d∣xf(d)F(x)=\sum_{d\mid x}f(d)则莫比乌斯反演公式为：f(x)=∑d∣xμ(d)F(xd)=∑d∣xμ(xd)F(d)f(x)=\sum_{d\mid x}\mu(d)F({x\over d})=\sum_{d\mid x}\mu({x\over d})F(

2016-07-11 16:00:43 2290

原创 HDU 5492 Find a path DP —— 2015 ACM-ICPC Asia Regional Hefei Online

原题见HDU 5492给一个N*M的矩阵，从(1,1)到(N,M)经过的格点分值分别为AiA_i，(路径只能向右或向下走,共N+M-1步)。求各种路径中最大的(N+M−1)∑N+M−1i=1(Ai−Aavg)2(N+M−1)\sum_{i=1}^{N+M−1}(A_i−A_{avg})^2.其中Aavg为Ai的均值A_{avg}为A_i的均值 N，M分别为1~30的整数，AiA_i为不超过30的整

2016-07-11 09:57:38 637

原创 HDU 5476 Explore Track of Point 2015上海网络赛几何题

发霉的题解….（自己已经不会证）题目见HDU 5476给出等腰三角形ABC，AB=AC，M为BC中点。P点为三角形内使min{∠MPB+∠APC,∠MPC+∠APB} 最大的点。求P点轨迹。则容易找到中线AM上的P点都满足使得∠MPB=∠MPC，∠APC=∠APB，则∠MPB+∠APC=∠MPC+∠APB=180° 故轨迹包含中线AM。并且所有满足的P点都应满足：∠MPB+∠APC=∠MP

2016-07-11 09:53:36 892

原创几种常用字符串算法

去年写的，已经发霉，出来晾一晾…KMP模式串P在母串S上的匹配。#define M 100000 int pre[M], n, m; void set(string p) { memset(pre, 0, sizeof(pre)); pre[0] = -1; int m = p.length(); for(int i = 1;i < m;i++) {

2016-07-10 15:42:53 4851

原创线性筛

线性筛用在素数、欧拉函数、莫比乌斯函数的打表上。埃拉托斯特尼筛法一开始最容易理解的筛法是酱紫滴~#define N 1000100 #define LL long long int num[N], prim[N]; int pn = 0; void table(){ memset(num, -1, sizeof(num)); for(int i = 2;i < N;i++) if(n

2016-07-10 15:33:00 1776 1

原创 HDU 5528 Count a - b 欧拉函数素数

原题见HDU 5528题意：f(m)f(m)表示abmodm≠0ab\mod m\neq 0的(a,b)(a,b)的个数，a,ba,b为小于m的非负整数。 G(n)=∑m|nf(m)G(n)=\sum\limits_{m|n}f(m) 给出n(1≤n≤109)n(1 \le n \le 10^9)，求G(n)mod264G(n)\mod 2^{64}.分析mod264\mod 2^{64}最后结

2016-07-10 14:12:49 1216

原创莫比乌斯函数的证明

遗忘是可怕的东西……好记性不如烂笔头讲真……命题现在假设我不知道什么是莫比乌斯函数，只知道F(x)=∑d∣xf(d)F(x)=\sum_{d\mid x}f(d)若已知F(x)F(x)，求f(x)f(x)的表达式。性质从已知的关系，可以得到性质： 1. 若y|x(y<x)y|x(y<x)，则F(y)F(y)包含的所有f(d)f(d)都被F(x)F(x)包含了，F(y)F(y)不能包含f(x)f(x

2016-07-09 15:20:44 4083 2

原创 CF 482 B. Interesting Array 线段树

敲了很久所以mark一下，辣鸡！原题见CF 482B试给出n个数，满足m个条件，每个条件即第l个数到第r个数这连续几个数的&值等于q。n和m范围为1到1051到10^5，q的范围为0到2302^{30}.分析同一个数如果要满足多个条件，则这个数值为q的或。每次满足的条件要对一个区间的值取或（旧值或新值），对于区间操作，利用线段树来操作。所有条件标记完以后，更新到所有结点。再反过来检验是否

2016-07-08 10:30:33 651

原创 POJ 3160 Father Christmas flymouse Tarjon+DP

原题见POJ 3160题意：已知一个有向图，每点均有权值。选择一个起点开始遍历，选择它可以经过的一些点，使得点权值之和最大。一个点可以经过多次，但权值只算一遍。权值可以为负。分析权值为负时，直接标记为0，即不会加这个点的权值，但是照样可以经过。强连通分量缩为一点，求出权值之和为该点权值。构造出若干个有向无环图后，从各个入度为0的点开始遍历。把父亲结点累加到该点的值放在一个额外数组c[]上，它不影响该

2016-07-07 20:10:53 579

原创连通分量涮涮乐

开心刷起题来吧！开涮对象——图论里的小小分支，连通分量概念省略一堆：有向图GG、强连通（两两可达）、强连通分量（极大强连通子图）转置图GTG^T：每个边反向求强连通分量Korasaju算法在图GG随便选择起点dfs，记录访问时间。在图GTG^T根据1得到的访问时间最大的点作为起点dfs，得到若干个树。每棵树在原图里都是一个强连通分量（由对称在转置图中同样也是）。这个算法给人的感觉就是

2016-07-06 15:27:38 1636

原创 HDU 5710 Digit-Sum 数学杂题

原题见HDU 5710题意：定义S(N)S(N)是NN的数位之和，给出一对a,b(0<a,b<101)a,b(0<a,b<101)，求是否存在nn满足a×S(n)=b×S(2n)a\times S(n)=b\times S(2n)分析：S(n)与S(2n)S(n)与S(2n)的关系如何？对于n中的任何一个数位x，若x为0-4，则因为没有进位，所以在S(2n)中贡献为2x;若x为5-9,则由于其超过1

2016-07-05 17:12:43 1859 4

Ugly Gardon