关闭
当前搜索:

[数学 FFT] Codechef July Challenge 2017 #APRPS Irrational Root

跟 51Nod 1356 代数数的次数 是一样的 不过这里都是质数 也就是就是 2n2^n 关键是输方案 这个不一定有二次剩余感谢sxt 一个一个数加进答案 转化成 已知F(x)F(x),求F(x+a√)F(x+\sqrt a)和F(x−a√)F(x-\sqrt a)的系数 这个推一下就是一个FFT拷了myy的板子 自己的太慢了QAQ#include #include<cs...
阅读(306) 评论(0)

[数学 二分图匹配] SRM 456 div1 FunctionalEquation

本来想自己再推一下的,但是退役了也就弃坑了 未经允许的搬了搬题人的题解// BEGIN CUT HERE #include #include // END CUT HERE #include #include #include #include #include #in...
阅读(234) 评论(0)

[拉格朗日乘数法 二分] BZOJ 2876 [Noi2012]骑行川藏

拉格朗日乘数法 ACdreamers [Math & Algorithm] 拉格朗日乘数法首先那个能量肯定是要花完的,就变成一个限制了,乘上拉格朗日乘子,求偏导,变成了 2λkix2i(xi−vi)=12\lambda k_ix_i^2(x_i-v_i)=1 ∑kisi(xi−vi)2=E\sum k_is_i(x_i-v_i)^2=E发现 x≥vx\ge v 且 x2(x−v)x^2(x-v...
阅读(258) 评论(0)

[最大似然估计 MLE] Codeforces 802DEF Helvetic Coding Contest 2017 D. E. F. Marmots

MLE最大似然估计可以通过观察数据来估计出随机的模型和参数 最大似然估计(Maximum likelihood estimation) Maximum likelihood estimation大概就是从所有可能的参数和模型中,选出概率最大的那一组模型,作为估计的模型和参数相关的还有一个最大后验估计(MAP),我没看过然后就很好写了D#include #include<cstdl...
阅读(329) 评论(0)

[IMO原题] BZOJ 4915 简单的数字题

第52届IMO试题解答 #include #include using namespace std;int main(){ long long l,r; cin>>l>>r; cout<<4<<endl<<max(r/11-l+1,0LL)+max(r/29-l+1,0LL); return 0; }...
阅读(356) 评论(0)

[二进制分组 dsu on tree 二次函数] Codechef KILLER Painting Tree

首先我们可以列一个dp 按深度 fuf_u表示链剖分中选择了从uu到当前深度的祖先这条链所能得到的最小答案 我们发现这条链的贡献是一个关于祖先深度二次函数 具体形式我忘了 而其他贡献则是一路上来兄弟子树中的最大值累加 这样就很明确了 是要求对每一个点维护一个二次函数 还要资瓷子树加 查询子树中x=kx=k最小值首先询问一坨二次函数最大值 有一个经典做法来自 BZOJ2646 我们可以维护nn...
阅读(159) 评论(0)

[可图性判定 Havel-Hakimi定理 构造] Codeforces Gym 100269 NEERC 13 Northern K. Kids in a Friendly Class

黑点和白点的个数可以从小到大枚举 黑白点之间的连边是显然的 那么对于白点或黑点内部 我们就需要连边使得每个点度数相同每次选两个剩余度最大的点相连是错误的 反例 6 2那么应该怎么连 需要Havel-Hakimi定理 由非负数组成的非增序列s:d1,d2,⋯,dn(n≥2,d1≥1)s:d_1,d_2,\cdots,d_n(n\geq2,d1\geq1)是可图的,当仅当序列s1:d2−1,d3...
阅读(633) 评论(0)

[莫比乌斯反演 伯努利数] BZOJ 2627 JZPKIL

ydc的题解 orzz 因为这道题才去学习了伯努利数 可以看这里#include #include #include using namespace std; typedef long long ll;namespace RHO{ int prime[9]={2,3,5,7,11,13,17,19,23}; unsigned long...
阅读(357) 评论(0)

[伯努利数] 51Nod 1258 序列求和 V4

自然数幂和的求法有很多 大概有 由二项式定理得递推式 拉格朗日插值法&牛顿插值法 第一类斯特林数 伯努利数 wiki传送门 ACdreamer的blog伯努利数有两种B+B^+和B−B^- 唯一的争议就在B1B_1上 因为 对于n>1n>1有B2n+1=0B_{2n+1}=0 其实这不是什么大问题 一般个人习惯用B−B^- 以下不注明BB都表示B−B^-然后我就只会伯努利数最裸的应用了...
阅读(944) 评论(0)

[高中数学联赛题] Codeforces Gym 101173 CERC 16 E. Easy Equation

首先(0,1,k)(0,1,k)是个解 然后用主元法 a2−k(b+c)a+b2+c2–kbc−1=0a^2 - k(b+c)a +b^2 + c^2 – kbc -1=0 可得如果(a,b,c)(a,b,c)为解,那么(k(b+c)−a,b,c)(k(b+c)-a,b,c)也为方程的解 直接大力bfs 其实你需要高超的高精度技巧 Python大法好 第一次知道原来set跑的这么快fro...
阅读(669) 评论(0)

[第二类斯特林数 组合计数] 省选模拟赛 2 B. 两弹一星 missile

题目大意 一张无向图的权值定义为 xkx^k,其中 xx 是图中结构为树的连通块个数。给定 nn, kk, 求出所有 nn 个点带标号的的简单无向图的权值和,对 998244353998244353 取模。 令 xix_i 表示树联通块 i 是否存在 , 图的权值为 (∑xi)m(\sum x_i)^m。对于某 kk 个联通块,如果同时出现,那么贡献为 S(m,k)∗k!S(m, k) * k!...
阅读(591) 评论(0)

[Bell数] HDU 4767 Bell & BZOJ 3501 PA2008 Cliquers Strike Back

ACdreamer大牛的详细介绍Bell数 大概是这样的 有递推公式 可以分治FFT求 生成函数很优美 可以用多项式科技求 是第二类斯特林数的和 可以用Bell三角形预处理 神奇的同余性质 可以计算对小质数取模的值 也可以CRT合并 这里的p是不大于100的素数 模素数p的周期 这里就是求bell模一个数的值 可以发现 95041567=31∗37∗41∗43∗4795041567 = 31*...
阅读(214) 评论(0)

[区间GCD预处理 树状数组 离线] HDU 5869 Different GCD Subarray Query

首先确定一个右端点 向左做后缀gcd的值是不超过log的 因为gcd必然递减 每次至少除以2 那么可以对每个右端点预处理出来 然后要求区间不同的gcd个数 这里要用到1878: [SDOI2009]HH的项链的技巧 用树状数组离线处理 按右端点的顺序处理询问 令pos[x]等于x这个值出现的最靠右的位置 显然这里是对计算贡献最有利的 然后用树状数组维护下 对于某个位置 有多少个pos的取值...
阅读(137) 评论(0)

[数学 树 || 环套树] BZOJ 3528 [Zjoi2014]星系调查

瞎扯的题目 直接上大神题解吧 #include #include #include #include using namespace std; typedef long long ll;inline char nc(){ static char buf[100000],*p1=buf,*p2=buf; if (...
阅读(156) 评论(0)

[随机化 || 黄金分割] HihoCoder #1461 Challenge 26 Rikka with Number

这个东西啊 逆过程是一个辗转相减 然后随机化以下最后状态另一个数就好了 比赛的时候被windows的rand坑死 windows下RAND_MAX是32768 被dls教导似乎windows下random_shuffle也有问题? 要不是rand不出来 就拿到奖品了233 #include #include #include #include using namespace std; i...
阅读(394) 评论(3)

[数学] NOI1995 极值问题

m、n为整数,且满足下列两个条件: ①m、n∈{1,2,…,K },(1≤K≤ 109) ②(n2-mn-m2)2=1 编一程序,由键盘输入K,求一组满足上述两个条件的m、n,并且使m2+n2的值最大。 例如,若K=1995,则m=987,n=1597,则m、n满足条件,且可使m2+n2的值最大。 打表可知满足fib数列 K 1...
阅读(170) 评论(0)

[构造 || 斐波那契分解] BZOJ 2912 [Poi1997]Jump

论文:何林--信息学中的守恒法 构造就不摘了 摘录fib分解的做法 #include #include #include #include #include using namespace std; typedef long long ll; const int con=100000000; class Int{ public:long long a[1000];...
阅读(180) 评论(0)

[复数 解析几何 || 随机化] Ural 1046 Geometrical Dreams & POJ 2600 Geometrical dreams

这个可以先设出a1 然后 一圈递推出an+1 an+1==a1 解这个方程就好了 这里用复数来实现 因为复数乘法是天然的旋转  复数旋转就是类似这样的东西 以下来自 POJ Discuss point getnp(point ori, point cen, double angle) //向量cen->ori绕cen旋转逆时针旋转angle弧度 { //复数平面向量旋转 poin...
阅读(154) 评论(0)

[黄金分割比] 黄金分割法求单峰函数最值

论文:杨思雨--美,无处不在——浅谈“黄金分割”和信息学的联系 #include #include #include #include using namespace std; const double eps=1e-4; const double phi=(sqrt(5.0)-1)/2; double A,B,C; double L,R; int cnt=0; i...
阅读(328) 评论(0)

[威佐夫博弈 黄金分割比 数学技巧] POJ 1067 取石子游戏 & HDU 1527 取石子游戏

论文:杨思雨--美,无处不在——浅谈“黄金分割”和信息学的联系 以下摘录证明 #include #include #include #include using namespace std; const double phi=(1+sqrt(5.0))/2.0; int main(){ freopen("t.in","r",std...
阅读(618) 评论(0)
67条 共4页1 2 3 4 下一页 尾页
    个人资料
    • 访问:309413次
    • 积分:12408
    • 等级:
    • 排名:第1311名
    • 原创:969篇
    • 转载:3篇
    • 译文:0篇
    • 评论:54条
    最新评论