- 博客(27)
- 收藏
- 关注
RSS订阅原创 多项式算法7：多项式三角函数
首先证明一下欧拉公式eiθ=cosθ+isinθe^{i\theta}=\cos\theta+i\sin\thetaeiθ=cosθ+isinθ由麦克劳林展开式，我们有ex=∑n=0∞xnn!e^{x}=\sum_{n=0}^{\infty}\frac{x^n}{n!}ex=n=0∑∞n!xn指数乘上虚数单位iiieix=∑n=0∞(ix)nn!=∑n=0∞(i)nxnn!e^{...
2019-06-02 11:55:45
2355
原创 多项式算法6：多项式快速幂
一开始本来想封装多项式乘法，按照普通快速幂的思路去做的Ak=(A⌊k2⌋)2A[(2,n)=1]A^k=\left(A^{\left\lfloor\frac{k}{2}\right\rfloor}\right)^2A^{\left[(2,n)=1\right]}Ak=(A⌊2k⌋)2A[(2,n)=1]但logk\log klogk还是太大，只好考虑将kkk转化为多项式的系数了。首先B...
2019-06-02 11:50:52
1531
原创 多项式算法5：多项式指数函数
先从牛顿迭代讲起。已知多项式函数G(z)G(z)G(z)，求多项式函数F(x)F(x)F(x)满足G(F(x))≡0(modxn)G(F(x))\equiv0 \pmod{x^n}G(F(x))≡0(modxn)考虑用迭代求解，假设我们已经求得F0(x)F_0(x)F0(x)满足G(F0(x))≡0(modx⌈n2⌉)G(F_0(x))\equiv0\pmod{x^{\left\lcei...
2019-06-02 11:47:39
2617
原创 多项式算法4：多项式对数函数
已知多项式F(x)F(x)F(x)，求解G(x)G(x)G(x)使得G(x)=lnF(x)G(x)=\ln F(x)G(x)=lnF(x)模板题，我们知道lnF(x)=∫F′(x)F−1(x) dx\ln F(x)=∫F'(x)F^{-1}(x)\ \text{d}xlnF(x)=∫F′(x)F−1(x) dx由指数函数的求导公式(xα)′=αx...
2019-06-02 11:42:44
638
原创 多项式算法3：多项式除法
已知次数界分别为n,mn,mn,m的多项式F(x)F(x)F(x)与G(x)G(x)G(x)，求解次数界为n−mn-mn−m的两个多项式Q(x)Q(x)Q(x)与R(x)R(x)R(x)使得F(x)=G(x)Q(x)+R(x)F(x)=G(x)Q(x)+R(x)F(x)=G(x)Q(x)+R(x)我们令Ar(x)=xDeg(A)A(x−1)A_r(x)=x^{\text{Deg}(A)}A(...
2019-06-02 11:38:25
4714
原创 多项式算法2：多项式求逆
已知多项式F(x)F(x)F(x)，求满足条件的G(x)G(x)G(x)使得F(x)G(x)≡1(modxn)F(x)G(x)\equiv1\pmod{x^n}F(x)G(x)≡1(modxn)考虑迭代求解令F(x)G0(x)≡1 (modx⌈n2⌉)F(x)G_0(x)\equiv1\ \pmod{x^{\left\lceil\frac{n}{2}\right\rceil}}F...
2019-06-02 11:33:15
1527
原创 多项式算法1：多项式乘法
多项式简介对于数域F\mathbb FF，若有∀i∈{1,2,3,⋯ ,n}\forall i\in\{1,2,3,\cdots,n \}∀i∈{1,2,3,⋯,n}，则f(x)=a0+a1x+a2x2+⋯+anx=∑i=1naixif(x)=a_0+a_1x+a_2x^2+\cdots+a_nx=\sum_{i=1}^na_ix^if(x)=a0+a1x+a2...
2019-06-02 11:27:26
5183
1
原创 2019.05.01【NOIP普及组】模拟赛C组 解题报告
原比赛链接T1：考场思路：原来是暴力模拟，只枚举偶数位的棋子，指针往前移动修改。后来想到用线段树优化，但是要进行大量的单点查询，所以比原算法还要慢。正解：我们可以开三个栈，记录左边界，右边界以及连续的值，最后再将连续区间合并即可，时间复杂度是线性的。T2：考场思路：暴力枚举字符串111的所有子字符串，然后利用判断子串的方式判断其是否为字符串222的子字符串，更新即可。时间复杂度...
2019-05-01 20:34:35
516
原创 题解 P1829 [国家集训队]Crash的数字表格 / JZPTAB
题目背景提示：原 P1829 半数集问题 已经迁移至 P1028 数的计算题目描述今天的数学课上，Crash小朋友学习了最小公倍数(Least Common Multiple)。对于两个正整数a和b，LCM(a, b)表示能同时整除a和b的最小正整数。例如，LCM(6, 8) = 24。回到家后，Crash还在想着课上学的东西，为了研究最小公倍数，他画了一张NM的表格。每个格子里写了一个数...
2019-04-30 21:55:50
290
原创 题解 P4466 [国家集训队]和与积
题目描述给出NNN，统计满足下面条件的数对(a,b)(a,b)(a,b)的个数：1≤a<b≤N1\le a<b\le N1≤a<b≤Na+ba+ba+b整除a×ba×ba×b输入输出格式输入格式一行一个数NNN输出格式一行一个数表示答案思路讲解比较好的莫比乌斯反演题了。ans=∑a=1N∑b=a+1N[(a+b)∣ab]\text{ans...
2019-04-30 21:51:27
454
原创 题解 P3911 最小公倍数之和
题目描述对于A1,A2,...,ANA_1,A_2,...,A_NA1,A2,...,AN，求：∑i=1N∑j=1Nlcm(Ai,Aj)\sum_{i=1}^N\sum_{j=1}^Nlcm(A_i,A_j)i=1∑Nj=1∑Nlcm(Ai,Aj)的值lcm(a,b)lcm(a,b)lcm(a,b) 表示a 和b 的最小公倍数输入输出格式输入格式第1行，1个整数N。第...
2019-04-30 21:47:58
322
原创 题解 P2260 [清华集训2012]模积和
题目背景数学题，无背景。题目描述求∑i=1n∑j=1m(nmod  i)(mmod  j),i≠j\sum_{i=1}^n\sum_{j=1}^m(n\mod i)(m\mod j),i=\not ji=1∑nj=1∑m(nmodi)(mmodj),i≠jmod 19...
2019-04-30 21:42:05
228
原创 题解 JZOJ 1355.珍珠分对
题目描述作为生日礼物,BESSIE得到了N (2 <= N <= 100,000;N%2 == 0)颗(假?)珍珠, 每颗珍珠的颜色为C种颜色之一 (1 <= C <= N).BESSIE发现N是一个偶数,她想把所有珍珠分成N/2对,使得每对的两颗珍珠的颜色都不同.BESSIE知道这样一个分配在给定数据里都是可能的,请帮助她找到这样一个分配. 如果有多个解,任意一个,输出...
2019-04-20 22:03:09
298
原创 题解 JZOJ 1354.土地购买
题目描述农夫John准备扩大他的农场,他正在考虑N (1 <= N <= 5,000) 块长方形的土地. 每块土地的长宽满足(1 <= 宽 <= 1,000,000; 1 <= 长 <= 1,000,000).每块土地的价格是它的面积,但FJ可以同时购买多快土地. 这些土地的价格是它们最大的长乘以它们最大的宽, 但是土地的长宽不能交换. 如果FJ买一块3x5...
2019-04-20 21:56:50
1286
1
原创 题解 JZOJ 1353.渡河问题
题目描述Farmer John以及他的N(1 <= N <= 2,500)头奶牛打算过一条河，但他们所有的渡河工具，仅仅是一个木筏。由于奶牛不会划船，在整个渡河过程中，FJ必须始终在木筏上。在这个基础上，木筏上的奶牛数目每增加1，FJ把木筏划到对岸就得花更多的时间。当FJ一个人坐在木筏上，他把木筏划到对岸需要M(1 <= M <= 1000)分钟。当木筏搭载的奶牛数目从...
2019-04-20 21:54:05
546
原创 题解 JZOJ 1352. 游荡的奶牛
题目描述奶牛们在被划分成N行M列(2 <= N <= 100; 2 <= M <= 100)的草地上游走，试图找到整块草地中最美味的牧草。Farmer John在某个时刻看见贝茜在位置(R1, C1)，恰好T (0 < T <= 15)秒后，FJ又在位置(R2, C2)与贝茜撞了正着。FJ并不知道在这T秒内贝茜是否曾经到过(R2, C2)，他能确定的只是，现在...
2019-04-20 21:50:41
438
原创 题解 TG2013 火柴排队
题目描述：涵涵有两盒火柴，每盒装有n根火柴，每根火柴都有一个高度。 现在将每盒中的火柴各自排成一列， 同一列火柴的高度互不相同， 两列火柴之间的距离定义为： ∑(ai−bi)2\sum(a_i-b_i)^2∑(ai−bi)2其中aia_iai表示第一列火柴中第iii个火柴的高度，bib_ibi表示第二列火柴中第iii个火柴的高度。每列火柴中相邻两根火柴的位置都可以交换，请你通过交换使...
2019-04-14 11:51:08
151
原创 题解 JZOJ 1350.流星雨(meteor)
题目描述：贝茜听说了一个骇人听闻的消息：一场流星雨即将袭击整个农场，由于流星体积过大，它们无法在撞击到地面前燃烧殆尽，届时将会对它撞到的一切东西造成毁灭性的打击。很自然地，贝茜开始担心自己的安全问题。以FJ牧场中最聪明的奶牛的名誉起誓，她一定要在被流星砸到前，到达一个安全的地方（也就是说，一块不会被任何流星砸到的土地）。如果将牧场放入一个直角坐标系中，贝茜现在的位置是原点，并且，贝茜不能踏上一块...
2019-04-14 11:41:18
540
原创 题解 JZOJ 1349.连线游戏(lines)
题目描述：Farmer John最近发明了一个游戏，来考验自命不凡的贝茜。游戏开始的时候，FJ会给贝茜一块画着N (2 <= N <= 200)个不重合的点的木板，其中第i个点的横、纵坐标分别为X_i和Y_i (-1,000 <= X_i <=1,000；-1,000 <= Y_i <= 1,000)。贝茜可以选两个点画一条过它们的直线，当且仅当平面上不存在...
2019-04-14 11:33:31
479
原创 题解 JZOJ 1348.越野跑(racing)
题目描述：为了能在下一次跑步比赛中有好的发挥，贝茜在一条山路上开始了她的训练。贝茜希望能在每次训练中跑得尽可能远，不过她也知道农场中的一条规定：奶牛独自进山的时间不得超过M秒(1 <= M <= 10,000,000)。整条山路被贝茜划分成T个长度相同的小段(1 <= T <= 100,000)，并且，贝茜用S_i表示第i个小段的路况。S_i为u，f，d这3个字母之一，...
2019-04-14 11:26:04
669
1
原创 题解 JZOJ 1347.晚餐队列安排(diningb)
题目描述：为了避免餐厅过分拥挤，FJ要求奶牛们分2批就餐。每天晚饭前，奶牛们都会在餐厅前排队入内，按FJ的设想，所有第2批就餐的奶牛排在队尾，队伍的前半部分则由设定为第1批就餐的奶牛占据。由于奶牛们不理解FJ的安排，晚饭前的排队成了一个大麻烦。第i头奶牛有一张标明她用餐批次D_i(1 <= D_i <= 2)的卡片。虽然所有N(1 <= N <= 30,000)头奶牛排...
2019-04-13 22:22:53
712
原创 题解 JZOJ 1164.【普及组模拟赛】除草(ontherun)
题目描述：一条笔直的路边有N(1 <= N <= 1,000)个草丛，草丛位置互不相同，我们用一个整数表示每个草丛的位置。现在你从某个位置L(1<=L<=1,000^2)出发去除掉所有的草丛，为达到目的你可以来回改变方向去移动，假设你以每秒1个单位距离的速度移动，并且能在到达草丛的瞬间把草除掉。要求计算草丛被清除的时刻总和的最小值。输入：第1行:两个空格隔开的整数N...
2019-04-06 20:21:30
649
原创 题解 JZOJ 1162.【普及组模拟赛】最大杂置(set)
题目描述：令SSS为nnn个元素的集合，则S有2n−12^{n-1}2n−1个子集（除去空集）。现在要你从这2n−12^{n-1}2n−1个子集中选出最多的子集，使这些子集能构成一个杂置。杂置是指任意两个集合没有包含或被包含的关系。例如对于有3个元素的集合{a,b,c}\{a,b,c\}{a,b,c}。{a,b}，{a,c}，{b,c}\{a,b\}，\{a,c\}，\{b,c\}{a,b}...
2019-04-06 18:55:05
215
1
原创 树链剖分
How to change the value of all the nodes passing through the shortest path of two nodes in a tree?Obviously this can be achieved by tree difference.How to find the sum of weights on the shortest pat...
2019-04-06 14:04:39
124
原创 伯努利数
参考文献资料：小蒟蒻yyb的博客百度百科-伯努利数伯努利数(Bernoulli)——学习笔记康复计划#3 简单常用的几种计算自然数幂和的方法Bernoulli number has been found for about 200 years.Its name is named after Bernoulli, its discoverer.Recursive definition：...
2019-02-13 19:06:55
725
原创 莫比乌斯反演
莫比乌斯函数是定义在N+\mathbb N^+N+上的函数μμμ考虑将nnn进行分解质因数：n=∏i=1maipin=∏_{i=1}^ma_i^{p_i}n=i=1∏maipi得到该函数的通解：μ(n)={1n=10∃i∈[1,m]&nbsp;pi≥2(−1)m∀i∈[1,m]&nbsp;pi&amp;lt;2μ(n)=\begin{cases}1&amp;amp;n=1\\0&am...
2019-02-13 10:37:09
103
原创 最小生成树
定义：生成树：在无向连通图G=(V,E)G=(V,E)G=(V,E)中，连通所有点并且没有回路的图为GGG的生成树。若该无向图有nnn个点，则它有n−1n-1n−1条边。若一棵生成树插入一条无向边(u,v)(u,v)(u,v)，则会形成一条回路。最小生成树：(MST)我们规定(u,v)(u,v)(u,v)为一条以uuu为起点，vvv为终点的边，w(u,v)w(u,v...
2019-02-12 21:58:59
242
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人