#include<iostream> #include<cstdio> #include<cstdlib> #include<cstring> using namespace std; ; int prime[maxn]; bool not_prime[maxn]; int main() { ; scanf("%d",&n); memset(not_prime,,sizeof(not_prime)); not_prime[]=tr…

deque<int>q; void spfa(int s) { ;i<=n;i++) d[i]=1e9; d[s]=; q.push_back(s); used[s]=; while(!q.empty()) { int x=q.front(); q.pop_front(); used[x]=; for(int i=first[x];i;i=next[i]) { int u=hh[i].t; if(d[u]>d[x]+hh[i].c) { d[u]=d[x]+hh[i].c; if(…

哎呀大水题..我写了一个多小时..好没救啊.. 数论板子X合一? 注意: 本文中变量名称区分大小写. 题意: 给一个\(n\)阶递推序列\(f_k=\prod^{n}_{i=1} f_{k-i}b_i\mod P\)其中\(P=998244353\), 输入\(b_1,b_2,...,b_n\)以及已知\(f_1,f_2,...,f_{n-1}=1\), 再给定一个数\(m\)和第\(m\)项的值\(f_m\), 求出一个合法的\(f_n\)值使得按照这个值递推出来的序列满足第\(m\)项的值为…

前言 因为loceaner太菜了,他什么东西都不会 所以他打算学一个东西就记录一下 不过因为他很菜,所以他不会写原理-- 而且,他希望在2019CSP之前不会断更 就酱紫,就是写给他自己的--因为他太菜了 别犯些傻逼错误好不好啊!你丢多少分了! 基础算法 小技巧 \[\sum_{i = 0}^{x}C(x, i)* C(y, i) = C(x + y, x)\] 使用负数下标 如何使用负数下标呢? 让数组前面有东西 int y[100]; int *z = y + 50; 这样的话调用\(z[-…

所以要什么时候才不会出现像这样的滚粗记呢? Day-?~Day-4 我也不清楚具体干了什么 不过学了很多东西,至少相对于去年还是个小菜鸡,今年已经变成大菜鸡了 Day-3~Day-1 几乎就是复习前面的内容,顺便搞了下圆方树一道题,但是复习感觉什么都不会了啊qwq Day0 坐上了去广州的高铁,在高铁上写(g)板(t)子(f) 一下高铁 \(\large\text{烫烫烫烫烫烫烫烫烫烫}\) 神tm27℃,直接成仙好吗 然后坐了1个半小时的地铁(站了一个半小时),终于来到了广二附近的酒店.然后想…

题目链接:CF1106F Lunar New Year and a Recursive Sequence 大意:已知\(f_1,f_2,\cdots,f_{k-1}\)和\(b_1,b_2,\cdots,b_k\),且有递推关系 \[ f_i=(\prod_{j=1}^kf_{i-j}^{b_j})\text%p \] 对于所有\(i>k\)均成立,给出\(f_n=m\),求\(f_k\),\(p=998244353\) 分析: 数论板子大集合 首先很显然有\(f_n=f_k^{q}\),于是考…

Preface 老叶说了高中停课但是初中不停的消息后我就为争取民主献出一份力量 其实就是和老师申请了下让我们HW的三个人听课结果真停了 那么还是珍惜这次机会好好提升下自己吧不然就\(AFO\)了 List Luogu P4198 楼房重建 把高度化为斜率,然后就是个动态最长上升子序列的问题了,线段树上二分即可解决,而且可以做到\(O(n\log n)\) NOIP模拟赛10.24 实力翻车,T1主席树裸题切了,T2想了贪心+前缀和+二分正解,最后1min写完发现忘记判边界了炸到60,T3以为很难…

P1125 笨小猴 题目描述 笨小猴的词汇量很小,所以每次做英语选择题的时候都很头疼.但是他找到了一种方法,经试验证明,用这种方法去选择选项的时候选对的几率非常大! 这种方法的具体描述如下:假设maxn是单词中出现次数最多的字母的出现次数,minn是单词中出现次数最少的字母的出现次数,如果maxn-minn是一个质数,那么笨小猴就认为这是个Lucky Word,这样的单词很可能就是正确的答案. 输入输出格式 输入格式: 输入文件word.in只有一行,是一个单词,其中只可能出现小写字母,并且长度…

相关链接: noip2018总结 noip2017是我见过的有史以来最坑爹的一场考试了. 今年北京市考点有一个是我们学校,我还恰好被分到了自己学校(还是自己天天上课的那个教室),于是我同时报了普及提高,一天半的时间都考了. 这次考试总的来说基本上都爆炸了.虽然都拿了一等奖,但这根本不能说明问题,从中可以看出我在敲代码学习上还是问题百出. 下面我分两篇来总结一个kubi的OIer的解题思路及心得,当然包括正解.然而我莫名其妙的打了好长好长……难道我太勤奋了? OIer常识:本文所有log均以2为底…

AC自动机已经足够棒了. 但是,好像有时还是要TLE的. 一般的AC自动还是比较好,如果在某些情况下还是会被卡掉,像是这个水题 考试的感觉 我看到这个题后,我清清楚楚的知道,这是个AC自动机+栈. 经过一番努力,把AC自动机打了出来,然后略加修饰,把栈补在里面. 我一复制,一粘贴,更高兴(?)了,过样例了,噫,我要A了. 但是我又随便写了一坨字符进去,然后就,不对了------ 最后也没调出来,后来我想想,好像是没恢复搜索栈顶元素. 我太菜了. ---以上都是废话--- 更改的过程 我T60之后…

[Noi2002]Savage 时间限制: 5 Sec  内存限制: 64 MB 题目描述 输入 第1行为一个整数N(1<=N<=15),即野人的数目. 第2行到第N+1每行为三个整数Ci, Pi, Li表示每个野人所住的初始洞穴编号,每年走过的洞穴数及寿命值. (1<=Ci,Pi<=100, 0<=Li<=10^6 ) 输出 仅包含一个数M,即最少可能的山洞数.输入数据保证有解,且M不大于10^6. 样例输入 3 1 3 4 2 7 3 3 2 1 样例输出 6 //…

从很早入坑ACM开始,便和各种算法的模板打着交道,虽然kaungbin的模板已经足够强大,但是自己在平常做题中也逐渐有着自己的一些模板,也有一些kuangbin模板中没有的更快的板子,虽然不确定时候以后会用到,但是能够记录下来形成自己的模板是最好的,,这样也对自己所学有个大致的总结,所以在搜寻网络上各种解决思路后,选择一位博主的解决方法: 利用python管理文件夹形式的模板库,生成 latex 的 .tex 文件,然后通过两次编译后得到 .pdf 的模板文件,,这样的好处是显而易见的,首先是代…

Day0 到杭州的时候是下午,休息了一下就吃完饭了. 晚上的时候闲着没事复习了一下几个不太熟的数论板子,\(\mathrm{exgcd}\),\(\mathrm{ExCRT}\),\(\mathrm{BSGS}\),还有高斯消元,线性基什么的,但愿不考. 感觉有点紧张. Day1 好早就到文渊中学了,感觉人特别多,进考场前膜拜了一下\(\mathrm{LHJKY}\)大佬,还从手机上看了一下\(\mathrm{windows}\)系统下的开栈方法,然后就进考场了. 地下室感觉有点闷热,还没有厕所…

As rich as Crassus 题目链接 题目描述 Crassus, the richest man in the world, invested some of his money with the Very Legitimate International Bank. The Bank offered a remarkable interest rate. They promised that given an initial investment of x, the balance…

题意 给定一个三角形和一个点 \(p\),如果该点不在三角形边上直接输出 \(-1\),否则在三角形上找一点 \(q\),使得线段 \(pq\) 平分三角形面积. 思路 看完题想都没想直接二分了. 就是比赛时复制粘贴的时候改错了,贡献了三发罚时. 板子来自 kuangbin 的模板. #include <bits/stdc++.h> using namespace std; typedef double db; const db eps = 1e-8; const db inf = 1e20;…

0 写在前面 本文受 NaVi_Awson 的启发,甚至一些地方直接引用,在此说明. 1 数论 1.0 gcd 1.0.0 gcd $gcd(a,b) = gcd(b,a\;mod\;b)$ 证明:设 $c\mid a$,$c\mid b$,则 $c\mid (b-a)$. 设 $c\nmid a$,则 $c$ 不是 $a,b-a$ 的公因子. 设 $c\mid a$,$c\nmid b$,则 $c$ 不是 $a,b-a$ 的公因子. int gcd(int a,int b){ if(!b) r…

0 写在前面 0.0 前言 由于我太菜了,导致一些东西一学就忘,特开此文来记录下最让我头痛的数学相关问题. 一些引用的文字都注释了原文链接,若侵犯了您的权益,敬请告知:若文章中出现错误,也烦请告知. 该文于 2018.3.31 完成最后一次修改(若有出错的地方,之后也会进行维护).其主要内容限于数论和组合计数类数学相关问题.因为版面原因,其余数学方面的总结会以全新的博文呈现. 感谢你的造访. 0.1 记号说明 由于该文完成的间隔跨度太大,不同时期的内容的写法不严谨,甚至 $LaTeX$ 也有许多…

#include<cstdio> #include<iostream> #include<cstring> #include<cmath> #include<algorithm> using namespace std; const int maxn=3000050; int phi[maxn],prime[maxn],tot; bool notprime[maxn]; void isprime(int x){ if(x<2)return…

Graveyard Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 1707   Accepted: 860   Special Judge Description Programming contests became so popular in the year 2397 that the governor of New Earck — the largest human-inhabited planet of the…

说句实话,我和并查集的缘分还是蛮深的,因为当年学完数论想着找板子题乱做(真是个神奇的找题方式呢),然后就看到了并查集QWQ,看了一会发现是图论就不看了,,,,,,结果还被说是大佬QWQ其实我只是个NaCl Fish而已QAQ 好了现在终于学了并查集,那我们就来总结总结 这一次总共是有三道题要讲 首先我们来看看板子题 P3367 [模板]并查集 并查集,在一些有N个元素的集合应用问题中,我们通常是在开始时让每个元素构成一个单元素的集合,然后按一定顺序将属于同一组的元素所在的集合合并,其间要反复查找…

题目: 1240 莫比乌斯函数 基准时间限制:1 秒 空间限制:131072 KB 分值: 0 难度:基础题 收藏 关注 莫比乌斯函数,由德国数学家和天文学家莫比乌斯提出.梅滕斯(Mertens)首先使用μ(n)(miu(n))作为莫比乌斯函数的记号.(据说,高斯(Gauss)比莫比乌斯早三十年就曾考虑过这个函数). 具体定义如下: 如果一个数包含平方因子,那么miu(n) = 0.例如:miu(4), miu(12), miu(18) = 0. 如果一个数不包含平方因子,并且有k个不同的质因子…

注:转载本文须标明出处. 原文链接https://www.cnblogs.com/zhouzhendong/p/Number-theory.html 数论算法 剩余系相关 学习笔记 (基础回顾,(ex)CRT,(ex)lucas,(ex)BSGS,原根与指标入门,高次剩余,Miller_Robin+Pollard_Rho) 本文概要 1. 基础回顾 2. 中国剩余定理 (CRT) 及其扩展 3. 卢卡斯定理 (lucas) 及其扩展 4. 大步小步算法 (BSGS) 及其扩展 5. 原根与指标入…

嗯,写这个是因为我太弱了\(ORZ\). #\(\mathcal{\color{silver}{1 \ \ Linear \ \ Sieve \ \ Method \ \ of \ \ Prime}}\)线性筛素数 嗯,其实对于这个而言,无非就是一个\(break\)不易理解而已. if(! (i % prime[j])) break ; 那么我们先来分析线性筛的复杂度,嗯,很显然因为是\(O(n)\)才称其为线性筛法.所以也就是说,对于每个合数,我们只让它被筛去一次.那么线性筛是如何保证的呢?…

来自蒟蒻 \(Hero \_of \_Someone\) 的 \(LCT\) 学习笔记 最近学了一波 \(LCT\) , 于是怒刷 \(LCT\) 合集...... $ $ 学的时候借鉴了 Clove_unique的博客 以及 PoPoQQQ的PPT 写得很详细,初学者可以去看看... $ $ 先甩一道板子题... $ $ 辣鸡得不能再辣鸡的题...只有 \(cut\) 和 \(link\) 操作, 然后就是判连通性... $ $ //made by Hero_of_Someone #inclu…

正解:数论 解题报告: 先,放个传送门QwQ 这题难点可能在理解题意,,, 所以我先放个题意QAQ 大概就是说,给定一个整数N,可以被拆成两个质数的成绩p*q,然后给出了一个数e,求d满足e*d=1(mod r),其中r=(p-1)*(q-1),最后还会给定一个c,求dc%N umm就是几个板子题的堆砌昂,,,首先pollard_rho找到pq求出r,然后逆元求出d,最后快速幂走一波 然后就做完辣!over! 然后这里注意一下,就我个人的习惯的话我很喜欢快速幂求逆元,,,因为很简单很无脑,,,但…

传送门 数论基础题. 对于第一种情况用快速幂,第二种用exgcdexgcdexgcd,第三种用bsgsbsgsbsgs 于是自己瞎yyyyyy了一个bsgsbsgsbsgs的板子(不知道是不是数据水了没卡如果有找出错的希望指正谢谢) 下面谈谈我对这个方法的理解. 实际上跟网上说的差不多. 要解ax≡bmod&ThinSpace;&ThinSpace;pa^x\equiv b\mod pax≡bmodp 相当于令p=k∗A+B,0≤B<pp=k*A+B,0\le B<pp=k∗A…

目录 数论 知识点 Exgcd 逆元 gcd 欧拉函数\(\varphi(x)\) CRT&EXCRT BSGS&EXBSGS FFT/NTT/MTT/FWT 组合公式 斯特林数 卡塔兰数 常用数学公式 技巧经验 容斥 组合计数 区间筛 博弈 有趣的式子 gcd有关 数论模板库 黑科技 \(long\ long\)相乘取模 子集枚举 高维前缀和 各种线性筛 高级算法 Exgcd Lucas EXCRT BSGS 高斯消元 线性基 裴蜀定理 FFT 拉格朗日插值 NTT FWT 数论 Tag…

,,,本来自以为,我dp学得还挺好的 然后今天一考发现都不会啊QAQ 连最基础的知识点都不清楚啊QAQ 所以就来写个题解嘛! 先放下板子题 其实我jio得,这题只要大概了解方法就不是很难鸭,,,毕竟是基础算法,还是比较好理解的QAQ 大概说下QwQ 就是开f[i],它是个一个长度为k的数组,表示的是重量为i时的前k优解(物品那一维显然可以滚掉不说 然后转移也很显然,就是f[i+val[j]]=max(f[i+val[j]],f[i]) 注意一下就是这里的取max是指从这k*2个数中选前k大的,彼…

一笔画问题 时间限制:3000 ms  |  内存限制:65535 KB 难度:4   描述 zyc从小就比较喜欢玩一些小游戏,其中就包括画一笔画,他想请你帮他写一个程序,判断一个图是否能够用一笔画下来. 规定,所有的边都只能画一次,不能重复画.   输入 第一行只有一个正整数N(N<=10)表示测试数据的组数.每组测试数据的第一行有两个正整数P,Q(P<=1000,Q<=2000),分别表示这个画中有多少个顶点和多少条连线.(点的编号从1到P)随后的Q行,每行有两个正整数A,B(0&l…

今天讲一讲数论吧(虽然清明讲过了) 进制转换 我们来看10这个数怎么转换成k进制 因为10=2^3+2^1,所以10就是1010 三进制也同理10=3^2+3^0,所以就是101 我们对于一个10进制数,就可以用短除法来求解 比如55的三进制 这里我们把所有的余数向上写一遍,其实代码实现的话就直接写一个栈就可以了,还有一个没啥用的注意事项,就是读的问题 怎么把一个k进制数转成十进制的X? 假设我们有一个k进制数X,总共有n位,那么他就是XnXn-1Xn-2......X1X0 我们第i位其实就对…