LzyRapX-CSDN博客

原创博客搬家QAQ

CSDN越用体验越差。。。。自己也越来越忙。。。就不在这里更新了。。。搬去这里了=>新博客github => github

2017-12-08 22:40:08 2024 2

原创（置顶）POJ题目分类

OJ上的一些水题(可用来练手和增加自信) (poj3299,poj2159,poj2739,poj1083,poj2262,poj1503,poj3006,poj2255,poj3094) http://www.cnblogs.com/kuangbin/archive/2011/07/29/2120667.html (kuangbin的分类)初期: 一.基本算法:

2016-05-03 13:50:59 1817

原创（置顶）HDOJ题目分类

HDOJ题目分类：模拟题, 枚举1002 1004 1013 1015 1017 1020 1022 1029 1031 1033 1034 1035 1036 1037 1039 1042 1047 1048 1049 1050 1057 1062 1063 1064 1070 1073 1075 1082 1083 1084 1088 1106 1107 1113 11

2016-04-01 13:15:47 1399

原创 Personal-Traning题解(一)

1.Codeforces Round #418 E (dp+组合数学)：https://loj.ac/problem/6094 http://codeforces.com/contest/814/problem/E 虽然题目是一样的，但是时限不一样，cf的时限大1s，所以在cf上这题O(n5)O(n^5)可以过，但是loj就不行了， codeforces：把图按照最短路分层，那么这个图需要

2017-10-09 11:02:28 13542

原创 HDU 1512 Monkey King（左偏树）

题目链接： HDU 1512题意：有 nn 只猴子，每只猴子都有一个力量，开始时互相都不认识，它们之间发生 mm 次争斗，每次发生a,ba,b发生争斗时，a,ba,b会从它们认识的猴子中选出一个最强的，并变为这两只猴子发生争斗，打完之后这两个猴子就互相认识，并且力量减半，如果a,ba,b互相认识就输出−1-1，否则输出认识的猴子中最大的力量值。题解：左偏树。AC代码：#include<cs

2017-09-21 18:52:24 966

原创 HDU 5628 Clarke and math （Dirichlet卷积）

题目链接： HDU 5628题意：给你f(i)。(i=1,2,…,n)f(i)。(i=1,2,…,n) 求g(i)=∑i1∣i∑i2∣i1∑i3∣i2⋯∑ik∣ik−1f(ik) mod 1000000007(1≤i≤n)\displaystyle g(i) = \sum_{i_1 \mid i} \sum_{i_2 \mid i_1} \sum_{i_3 \mid i_2} \cdots

2017-09-03 14:08:54 4154

什么叫二次剩余，其实就是对于给定的p(p∈P)和np(p∈P)和n，如果有xx满足x2≡n(modp)x^2≡n(\mod p)，那么nn在模pp意义下就是二次剩余。其实就是模意义下能否开根号。我们先定义Fp\mathbb {F}_p，这是一个数域，其实就是00到p−1p−1这pp个数与模pp意义下加减乘除运算构成的集合。定理1：对于x2≡n(modp)，x^2≡n(\mod p)，总共有p−12\

2017-09-01 13:02:21 6801

原创 ACdream 1121 喵喵的IDE （Trie）

题目链接： ACdream 1121题意：中文题意。。。题解：可以从缓存中，1.拿出字符串进行操作。2.删除末尾字符。3. 在末尾加入任意字符。让你求使得最小步数获得新的字符串并加入缓存中。可以用字典树(Trie)来存储字符串并进行查找。AC代码：#include <bits/stdc++.h>using namespace std;int ch[1000005][27],va

2017-08-30 18:40:25 579

原创 ACdream 1118 Captain Hammer （物理题）

题目链接： ACdream 1118题意：就是一个高中物理模型。给你一个速度，距离，让你求速度与重力加速度的夹角。题解：设角度为 aa。根据物理公式可得 V∗sinag=12∗(Dv∗cosa),\frac{V*sina}{g}=\frac{1}{2}*(\frac{D}{v*cosa}), 可得 sin（2a）=D∗gV2sin（2a）=\frac{D*g}{V^2}，故角度

2017-08-30 11:05:28 496

原创 ACdream 1116 Gao the string! （扩展KMP+矩阵快速幂）

题目链接： ACdream 1116题意： give you a string, please output the result of the following function mod 10000000071000000007 ∑i=n−1i=0f(a[i])\sum_{i=0}^{i=n-1}f(a[i]) nn is the length of the string。 f(i)f

2017-08-29 18:30:48 467

原创 ACdream 1114 Number theory （莫比乌斯反演）

题目链接： ACdream 1114题意：给你一个序列 aa，求出这个序列中互质数的有多少对。其中所有的整数的都小于等于 222222222222 。题解：f(d)f(d) 为 gcdgcd 恰好为 dd 的数的对数。 F(d)F(d) 为 gcdgcd 为 dd 的倍数的对数。 μ(d)μ(d) 表示莫比乌斯函数。所以， F(d)=∑f(n)。其中（n%d==0)F(d) = \s

2017-08-28 20:53:25 577

原创《数论概论》读书笔记第6章线性方程与最大公约数

这章讲的就是欧几里得算法和exgcdexgcd。原式: ax+by=gcd(a,b)ax+by=gcd(a,b)(假设a≥ba≥b) 当 b=0b=0 时有 gcd(a,b)=agcd(a,b)=a,此时 x=1,y=0x=1,y=0 当 bb 不为 00 时,根据欧几里得定理 gcd(a,b)=gcd(b,a%b)gcd(a,b)=gcd(b,a\%b) 可得ax+by=gcd(a,b)=

2017-08-27 14:24:26 11570 1

原创 LOj #2002. 「SDOI2017」序列计数（容斥+dp+矩阵快速幂）

题目链接： LOj 2002题意：要求得到一个长度为 nn 的序列，序列中的数都是不超过 mm 的正整数，而且这 nn个数的和是 pp 的倍数。这nn个数中，至少有一个数是质数。问你有多少个序列满足要求。题解：根据容斥原理，因为题目要求至少有一个数是素数，用所有方案减去不含质数的方案就是答案。设 dp[i][j]dp[i][j]表示序列前 ii 个数模 pp 的余数为 jj 时的方案

2017-08-19 11:39:27 664

原创 LOj #6164. 「美团 CodeM 初赛 Round A」数列互质 (莫队算法)

题目链接： LOj 6164题意：给出一个长度为n n 的数列，a1,a2,a3,...,ana_1,a_2,a_3,...,a_n以及 mm 组询问 (li,ri,ki)(l_i,r_i,k_i)，求区间[li,ri] [l_i,r_i]中有多少数在该区间中的出现次数与 kik_i互质。题解：莫队算法。时间复杂度：O(n∗n√)O(n*\sqrt{n})。AC代码：#include <b

2017-08-15 11:15:50 985

原创 LOj 2000「SDOI2017」数字表格（莫比乌斯反演）

题目链接： LOj 2000题意：给你 nn 和 mm。让你求 ans(n,m)=∏ni=1∏mj=1fib[gcd(i,j)]ans(n,m)= \prod_{i=1}^{n} \prod_{j=1}^{m} fib[gcd(i,j)]。题解：因为要求：ans(n,m)=∏ni=1∏mj=1fib[gcd(i,j)]ans(n,m)= \prod_{i=1}^{n} \prod_{j=1

2017-08-09 22:36:31 1368

原创 2017 Multi-University Training Contest-第三场 Kanade's convolution（FWT）

题目链接： HDU 6057题意：给你 A[0..2 m −1] A[0..2^m-1] 和 B[0..2 m −1] B[0..2^m-1]。先让你求 C[0..2 m −1] C[0..2^m-1]。 C[0..2 m −1] C[0..2^m-1]的运算方式为： C[k]=∑ i and j=k A[i xor j]∗B[i or j] C[k]=\sum_{i~and~j=k}

2017-08-03 11:41:40 659

原创【概率问题】开保险箱（结论题）

题目： 500 500个保险箱，500 500把钥匙，每把钥匙恰好能打开一个保险箱。每个保险箱也只有一把钥匙能打开。钥匙锁在保险箱里，每个保险箱里一把钥匙，放法是随机的。先打破两个箱子取出钥匙。问其余的保险箱均可用钥匙打开的概率是多少？题解：在钥匙已经放好之后，把箱子排成一排，使得箱子里面的钥匙刚好可以开下个箱子。很明显，这样排好之后有可能出现几个环，例如有三个箱子的钥匙是互相可以开的，但

2017-07-31 17:06:55 1404

原创 BZOJ 2820 YY的GCD （莫比乌斯反演）

题目链接： BZOJ 2820 权限题Description 神犇YY虐完数论后给傻×kAc出了一题。给定N,M N, M,求1<=x<=N,1<=y<=M 1<=x<=N, 1<=y<=M且gcd(x,y) gcd(x, y)为质数的(x,y) (x, y)有多少对？kAc这种傻×必然不会了，于是向你来请教……多组输入 Input 第一行一个整数T 表述数据组数接下来T行，每行两个正整

2017-07-30 18:13:17 615

原创 BZOJ 2301 [HAOI2011]Problem b （莫比乌斯反演）

题目链接： BZOJ 2301题意：对于给出的n n个询问，每次求有多少个数对(x,y) (x,y)，满足a≤x≤b，c≤y≤d，且gcd(x,y)=k a≤x≤b，c≤y≤d，且gcd(x,y) = k。 100% 100\%的数据满足：1≤n≤50000，1≤a≤b≤50000，1≤c≤d≤50000，1≤k≤50000 1≤n≤50000，1≤a≤b≤50000，1≤c≤d≤50000

2017-07-29 15:55:26 763

原创 ACdream 1104 瑶瑶想找回文串 (Splay+二分+hash)

题目链接： ACdream 1104题意：中文题面….题解： Splay+二分+hash。AC代码：/** this code is made by LzyRapx* Problem: 1104* Verdict: Accepted* Submission Date: 2017-07-28 18:18:18* Time: 708MS* Memory: 8116KB*/#inc

2017-07-28 18:37:51 439

原创 2017 Multi-University Training Contest - 第二场 09 TrickGCD （容斥+莫比乌斯）

题目链接： HDU 6053题意：给出你一个长度为n n的A A数组，让你构造出一个长度也为n n的B B数组，且B B数组要满足对于所有的1<=i<=n,A[i]<=B[i] 1<=i<=n,A[i]<=B[i]，且对于B B数组，任意一个区间的gcd>=2 gcd>=2，求满足的方案数。题解：容斥原理+莫比乌斯函数。答案就是∑ min(A) i=2 −mu(i)∗∏(⌊cnt i

2017-07-28 11:34:17 479

原创 2017 Multi-University Training Contest - 第二场 07 If the starlight never fade （数论）

题目链接： HDU 6051题意： f(i) f(i) 表示满足(x+y) i ≡x i %p (x+y)^{i}≡x^{i}\% p 中 (x,y) (x,y) 的对数，其中，1≤x≤p−1,1≤y≤m. 1≤x≤p−1,1≤y≤m. 给你一个非负数 m m 和素数 p p，让你求∑ p−1 i=1 if(i) \sum_{i=1}^{p-1}if\left(i\right)，答案mo

2017-07-27 23:00:02 680

原创 ACdream 1079 Walking in the Rain （线性dp）

题目链接： ACdream 1079题意：有 n n个地板，你可以从 i i 跳到i+1 i +1 也可以跳到i+2 i + 2 ,我们的任务是从 i i 跳到 n n 下雨了，第i i个地板能够坚持的天数为a[i] a[i] 。问你最多多少天以后我们还可以跳到第n n个地板。题解： dp dp。跳到 n n 点的有 n−1 n-1 点和 n−2 n

2017-07-26 14:31:00 536 6

原创 ACdream 1112 Alice and Bob （SG函数）

题目链接： ACdream 1112题意：给你n n堆石子，每次要按照规则将该堆石子的个数分为两堆或者将该堆石子的个数减少，谁不能继续操作了谁就输。规则如下：若n=a∗b n=a*b 且a>1，b>1 a>1，b>1，可以将石子分为 a，b a，b 两堆；或者可以减少为 a a 或 b b 。比如n=6 n=6，可以分为(2,3)，2,3 (2,3)，2,3这三种情况，但是不能分为(1

2017-07-26 12:18:27 430 2

原创 2017 Multi-University Training Contest - 第一场 04 Division Game （NTT+数学）

题目链接： HDU 6036题解： (官方：题解)显然每个石子堆最多做 ∑ m i=1 e i \sum_{i = 1}^{m}{e_i} （记为 w w ）次操作。此外，如果定义一个堆做 x x 次操作恰好变为 1 1 的方案数为 f(x) f(x) ，显然每个数字做少于 x x 次操作不变为 1 1 的方案数也是 f(x) f(x) 。为了统计结束于石子堆i i的情况数，我们可以枚举这是它

2017-07-26 10:33:41 866

原创 ACdream 1110 True love （多重背包+dp）

题目链接： ACdream 1110题意：给你一些物品的体积和对应的数量，求可以拿走多少体积不一样的物品，且不超过背包的容量。题解：多重背包呗。 dp dp 题。设 dp[i] dp[i] 表示容量为 i i 的种类数量。那么容易得到转移方程： dp[j+a[i]]=dp[j]+1 dp[ j + a[ i ] ] = dp[j] + 1 且 j+a[i]<cap j + a[

2017-07-23 10:41:27 451

原创 ACdream 1107 瑶瑶迷上新游戏（计算几何+确认抛物线）

题目链接： ACdream 1107题意：给你一系列的点，问你这些点是否可以确认一条抛物线。题解：先将这些点按横坐标从小到大排序。我们设抛物线：y=ax 2 +bx+c y=ax^2+bx+c。三点可以求出一条抛物线。再求出 a a，b b，c c。最后将剩余的点代入到已求出的抛物线，验证一下即可。注意：给出的点不能太靠拢，否则求出的抛物线就不准确，就No No。也就是

2017-07-22 13:24:56 457

原创 ACdream 1078 Cutting Figure （dfs）

题目链接： ACdream 1078题意：给你一个n∗m n*m的图， # 表示存在的点，. . 表示无，#被保证全都连在一起，求最少需要把多少个#转换成，使得它分成两个部分。题解：其实只有三种情况。第一种是‘#’数量少于3 3的情况，这时候不可分割。第二种是只需要割一个点就可以的。最后一种情况是只需要割两个点就可以。然后判断‘#’个数，如果小于3 3，直接输出−1 -1

2017-07-21 10:42:59 369

原创 ACdream 1076 XXX的机器人（dp+置换+康拓展开）

题目链接： ACdream 1076题意：中文题面。。。题解： dp dp 题。设dp[i][j] dp[i][j]表示指令为i i 时，全排列状态为j j 。因为每条指令可以选择执行或不执行。那么 dp dp 的转移方程就是： dp[i+1][x]=min(dp[i+1][x],dp[i][j]+|t−s|+1) dp[i+1][x] = min(dp[i+1][x],dp[i

2017-07-20 10:50:49 474

原创 ACdream 1072 Kill The Monster （dp）

题目链接： ACdream 1072题意：中文题面….题解： dp dp 题。设dp[t][u][i][j] dp[ t ][ u ][ i ][ j ]表示在时间 t t 时在地图上 u u 位置上用了 i i 次技能并且这时的技能还需要 j j 的冷却时间。然后推一下转移方程就可以了。注意有坑：会有输入相同的时间和位置。AC代码：/** this code is made

2017-07-17 10:48:20 331

原创 ACdream 1063 平衡树（Trie树）

题目链接： ACdream 1063题意：让你设计一种数据结构，可以支持插入一个数，可以在这个数据结构里查询其中哪个数和给定的数的异或值最大或最小。题解： trie树。建立一个trie树，按01进制存进树里面，在查询异或最值的时候，我们可以在trie树上一直走就可以了。AC代码：/** this code is made by LzyRapx* Problem: 1063* Ver

2017-07-15 23:36:44 424

原创 ACdream 1060 递推数（矩阵快速幂+找循环节）

题目链接： ACdream 1060题意：已知A(0)=0,A(1)=1,A(n)=3∗A(n−1)+A(n−2)(n≥2) A(0) = 0 , A(1) = 1 , A(n) = 3 * A(n-1) + A(n-2) (n ≥ 2). 求 A(A(A(A(N))))Mod(1e9+7) A(A(A(A(N)))) Mod (1e9 + 7).题解：这种题一看就是用矩阵乘法来做啊。

2017-07-14 11:58:59 486

原创 ACdream 1119 瑶瑶的动感光波（加强版）(LCA)（背包dp）

题目链接： ACdream 1109题意：中文题意….题解：如果这题按照 ACdream 1102 题解去做，肯定会TLE。我们先DFSDFS预处理出全部结点的父节点，深度所有点对的LCALCA（最近公共祖先）。枚举以每个结点开始到根结点这段路径上每一段长度上背包容量为0 ~ 50的选取情况的最优价值。对于每次询问，可以拆分成两部分，一部分是xx到zz且包括zz的一段路，即：dep[x]−d

2017-07-12 19:28:16 473

原创 ACdream 1103 瑶瑶正式成为CEO （树链剖分+最小费用最大流）

题目链接： ACdream 1103题意：中文题面….题解：因为 C C 修改操作要在树上修改边权，我们用树链剖分去解决就好了。然后对于 Q Q 查询操作，我们可以用最小费用最大流去求最小费用就可以了。建图问题：假设每条边都跑了 0 0 的流量，我们先算出跑了 0 0 的流量的费用，然后对于一条边，当跑了小于 a i a_i 的流量的时候，每次增加一点流量，就相当于减小了 c i

2017-07-11 19:41:34 438

原创 ACdream 1102 瑶瑶的动感光波（树形dp）（背包dp）

题目链接： ACdream 1102 题意：自己点链接进去看吧…中文题面。。。题解：先预处理求出每个结点的父节点，然后在这棵树上做背包 dp dp 就可以了。这题还有加强版。题目链接： ACdream 1119 加强版题解： ACdream 1119 题解做法有点复杂。这题代码交上去完全会TLE。等等写一下加强版吧。AC代码：/** this code is made by L

2017-07-10 12:07:28 446

原创 ACdream 1101 瑶瑶想要玩滑梯（线段树）（区间替换+区间查询LCIS）

题目链接： ACdream 1101题意：给你一个区间，每次修改可以修改一段区间的值，也可以询问一个区间，问你这个区间中（严格上升）最长连续上升的列长度(LCIS) (LCIS)有多长？题解：线段树。其实就是区间替换和区间查询。调了好久，太久不写会变菜。这是个不争的事实。AC代码：/** this code is made by LzyRapx* Problem: 1101* V

2017-07-09 15:38:18 374

原创 ACdream 1070 神奇的%系列二 (数学+二分)

题目链接： ACdream 1070题解: 直接用lower_bound()和upper_bound()去二分。AC代码：/** this code is made by LzyRapx* Problem: 1070* Verdict: Accepted* Submission Date: 2017-07-09 13:15:09* Time: 2508MS* Memory: 8264

2017-07-09 13:24:36 404

原创 ACdream 1055 Crayon （离散化+树状数组）

题目链接： ACdream 1055题意：给你 nn 个操作，这些操作可以画一条区间大小为 [L,R][L,R] 的线段，可以删去第 ii 条加入的线段，也可以询问区间 [L,R][L,R] 中有多少个整数点是这些线段共有的。题解：先离线处理一下，再用树状数组维护一下区间变化即可。详细看代码。AC代码：#include <bits/stdc++.h>using namespace

2017-07-04 20:39:09 339

原创 BZOJ 3224 Tyvj 1728 普通平衡树 (Splay)

题目链接： BZOJ 3224题意：让你实现一棵树，实现插入，删除，查询xx数的排名，查询排名为xx的数，求xx的前驱(前驱定义为小于xx，且最大的数)，求xx的后继(后继定义为大于xx，且最小的数)的功能。题解：套平衡树Splay。AC代码：#include <cstdio>#include <cstring>#include <iostream>using namespa

2017-06-27 13:39:05 452

原创 ACdream 1006 Mengzhu （数学推导）

题目链接： ACdream 1006题意：题解：纯数学推导。因为： log2(x+y)=Alog2(x+y)=A log2(x−y)=Blog2(x-y)=B 所以， 2A=x+y2^{A} = x+y 2B=x−y2^{B} = x-y相加一下，2x=2A+2B2x=2^{A}+2^{B} 所以，x=2A+2B2=2A−1+2B−1x=\frac{2^{A}+2^{B

2017-06-23 15:51:33 378

空空如也

VR（Virtual Realty）所要面对的挑战有哪些? VR在未来的发展将会如何?