关闭
当前搜索:

[数论] LOJ #508. 「LibreOJ NOI Round #1」失控的未来交通工具

这种非简单路 一般转化成任意一条路加上若干环 这里大概是任意一条路加上若干环长的gcd 任意一条路 可以弄出任意一颗生成树? 但是这里实际上只需要一个带权并查集 详见官方题解#include #include #include using namespace std; typedef long long ll; typedef pair<...
阅读(204) 评论(0)

[数论] LOJ #510. 「LibreOJ NOI Round #1」北校门外的回忆

这个题跟树状数组没有半毛钱关系 首先这是一个最低位翻倍的过程,如果这一位最终会变成 00 ,那么步数是 O(logn)O(\log n)的? 要是这一位不能变成 00 ,也就是在环上跑了,似乎跑到环上的步数也是 O(logn)O(\log n)的? 然后就变成了,修改要是跑到 00 就暴力跑,跑到环上,询问照样 设某个修改最早到环上的起点是w+z×kmw+z\times k^m,那么会对 w...
阅读(208) 评论(0)

[数论] Codeforces 819D R #421 D.Mister B and Astronomers & 516E R #292 E. Drazil and His Happy Friends

两道类似的题 819D考虑一个人应该能够观察的位置 ti,(ti+S)modT,(ti+2S)modT⋯t_i,(t_i+S)\bmod T,(t_i+2S)\bmod T\cdots 这个应该是形成 gcd(S,T)\text{gcd}(S,T) 个环,每个环是长度 TgT\over g 然后把同一个环的一起处理,把点放到环上,那么沿环的方向到下一个点为止应该都是归到这个点答案里面的#inc...
阅读(206) 评论(0)

[反演 数论] 51Nod 1355 斐波那契的最小公倍数

我好菜啊 出过一万遍的原题 我怎么第一次看见啊 某乎链接 按照zyz的做法 orzzlcm(fS)==∏T⊆S,T≠∅gcd(fT)(−1)|T|+1∏T⊆S,T≠∅f(−1)|T|+1gcd{T}\begin{eqnarray} \text{lcm}(f_S)&=&\prod_{T\subseteq S,T\neq \emptyset }\text{gcd}(f_T)^{(-1)^{|T|+1...
阅读(567) 评论(1)

[数论 DAG最长路] Codeforces 772C VK Cup 2017 - Round 2 C. Vulnerable Kerbals

ax≡b(modm)ax\equiv b \pmod m 有解就是 (a,m)|b(a,m) | b 也就是 (a,m)|(b,m)(a,m)|(b,m) 那么把数按照和mm的gcd的整除关系 建成DAG 跑一个最长链#include #include #include #include #define pb push_bac...
阅读(138) 评论(0)

[数位DP Lucas定理] 2017 计蒜之道 复赛 E. 商汤智能机器人

阿爷教导我 ∑∞i=0(Ai)∗(A+B−iA),A=x+y2,B=x−y2\sum _{i=0}^\infty \dbinom{A}{i}*\dbinom{A+B-i}{A},A={x+y\over 2},B={x-y\over2} 然后就是数位dp+lucas定理的套路了 注意有减法要处理退位 一开始写的时候考虑的有点问题 调了很久很久 这个可以从低到高 也可以从高到低 因为意识模糊就...
阅读(543) 评论(3)

[期望与方差 Lucas定理] 算法马拉松25 D 小Q的集合

大概是 ∑Cin∗2∗ik∗(ik−(n−i)k)modm\sum C_n^i*2*i^k*(i^k-(n-i)^k) \bmod m 后面那个东西模mm有循环节 同一组的组合数的和用lucas定理也能化成组合数乘上2的幂#include #include #include #include using namespace...
阅读(442) 评论(0)

[数论] 2017 计蒜之道 初赛 第一场 阿里天池的新任务

我们令子串Ss,tS_{s,t}对应的wsw_s为这个子串的bb 因为互质 所以bb是互不相同的 转为求有几个bb满足条件 然后每一位看作一个限制 把所有限制离散化取交即可 注意处理奇偶性 还有一个细节 s≤n−m+1s \leq n-m+1 我们还要倒着把最后m−1m-1个bb中合法的减掉#include #include #include...
阅读(147) 评论(0)

[扩展欧拉定理] BZOJ 4869 [Shoi2017]相逢是问候

欧拉定理 aϕ(n)≡1(modn),(a,n)=1a^{\phi(n)} \equiv 1 \pmod n,(a,n)=1 进而有 ab≡abmodϕ(n)(modn),(a,n)=1a^b \equiv a^{b \bmod \phi(n)} \pmod n,(a,n)=1 那么对于aa和nn不互质 我们有扩展欧拉定理 ab≡abmodϕ(n)+ϕ(n)(modn),(a,n)≠1且b≥ϕ...
阅读(578) 评论(0)

[数论 反演] BZOJ 4833 最小公倍佩尔数

当时比赛时灵机一动 把gigi−1g_i\over g_{i-1}喂给了OEIS 然后就找到了 233 就是这个咯 然后就水过去了 题解?题解我还没看 先挖个坑UPD:跟这个题是一毛一样的咯#include #include #include using namespace std; typedef long long ll;const...
阅读(760) 评论(0)

[组合数取模] BZOJ 4830 [Hnoi2017]抛硬币

习惯性交换aa和bb 令b≥ab\geq a 首先特判a=ba=b 这时答案为22a−Ca2a2{2^{2a}-C_{2a}^a}\over 2 其实就是所有情况减去平局的情况 剩下的不是A赢就是B赢 且是对称的 那么除以2 ∑Cin∗Cin=∑Cin∗Cn−in=Cn2n\sum C_n^i*C_n^i=\sum C_n^i*C_n^{n-i}=C_{2n}^n然后如果b>a 我们考...
阅读(1275) 评论(2)

[离散对数 FFT FastIO] ICPC 2016 Shenyang K & HDU 5958 New Signal Decomposition

特判0 然后将所有下标取对数 那么就是 bk′=∑ah′∗ch′+k′b_{k'}=\sum a_{h'}*c_{h'+k'} 使用FFT解决本来好好的一道题 卡常数+卡空间+卡精度?!#include #include #include #include using namespace std; typedef doubl...
阅读(421) 评论(2)

[反演 莫队算法] HDU 4676 Sum Of Gcd

考虑对一堆数求一个gcd 我们记录每一个约数xx的出现次数num[x]num[x] 那么两两间gcd的和是什么呢 ∑x∗C2num[x]\sum x*C_{num[x]}^2 肯定不是 这样会有重复 猜想f(x)f(x)是关于x的函数 ∑f(x)∗C2num[x]\sum f(x)*C_{num[x]}^2 那么应该满足 ∑d|xf(d)=x\sum_{d|x} f(d)=x 因为...
阅读(156) 评论(0)

[数论模板] BZOJ 3667 Rabin-Miller算法

还是发一下吧 以后找版子容易些 多年前的代码丑 #include #include #include #include #include using namespace std; typedef long long ll; inline char nc() { static char buf[100000],*p1=buf,*p2=buf; if (p1==p2) { p...
阅读(184) 评论(0)

[置换快速幂 中国剩余定理] POJ 1282 庆典的日期

论文:潘震皓--置换群快速幂运算 研究与探讨 置换的快速幂 题解 n2logn #include #include #include #include #define cl(x) memset(x,0,sizeof(x)) using namespace std; typedef long long ll; inline char nc(){ static...
阅读(184) 评论(0)

[斐波那契拆分 乱搞 数学] 51Nod 1350 斐波那契表示

可以证明每次选最大的fib数减去一定是最优解之一 然后预处理G[fib[i]-1]  额 代码横跨noip前后 意识模糊 一片mess #include #include #include #include using namespace std; typedef long long ll; ll fib[205]; int maxn; ll F[205]; inline...
阅读(294) 评论(0)

[数论] 51Nod 1217 Minimum Modular & Codeforces 303C #183 (Div. 1) Minimum Modular

N个不同的数a[1],a[2]...a[n],你可以从中去掉K个数,并且找到一个正整数M,使得剩下的N - K个数,Mod M的结果各不相同,求M的最小值。 来自讨论帖:http://www.51nod.com/question/index.html#!questionId=994 这个题我们可以考虑从小到大枚举m(从max(1,n-k)到max(a[i])+1),然后判断能否...
阅读(401) 评论(0)

[五边形数定理 DP] 51Nod 1259 整数划分 V2 & HDU 4651 Partition

详见:http://blog.csdn.net/acdreamers/article/details/12259815 #include #include #include using namespace std; typedef long long ll; const int N=50005; const int P=1e9+7; ll F[N],f[N];...
阅读(233) 评论(0)

[数论] 51Nod 1225 余数之和

n%i=n-[n/i]*i #include #include #include using namespace std; typedef long long ll; int P=1000000007; int INV=5e8+4; ll n; ll Ans=0; inline ll Sum(ll l,ll r){ return ((l+r)%P)*((r-l+1)%P)%P*...
阅读(163) 评论(0)

[DP 数论] 51Nod 1201 整数划分

dp[i][j] 代表用j个数字的和为i的数的个数 dp[i][j] = dp[i - j][j - 1] + dp[i - j][j](不过好像不太对) 把所有的数都加上1,然后再额外的给出一个1  或 将现有的每一位加1 #include #include #include #include using namespace std; const int N=50005;...
阅读(124) 评论(0)
52条 共3页1 2 3 下一页 尾页
    个人资料
    • 访问:309411次
    • 积分:12408
    • 等级:
    • 排名:第1311名
    • 原创:969篇
    • 转载:3篇
    • 译文:0篇
    • 评论:54条
    最新评论