LeoJAM的博客

好难啊如果p为质数自然随便做但是不是质数需要扩展Lucas定理说白了就是对唯一分解的质数单独讨论再用中国剩余定理合并#include<bits/stdc++.h> using namespace std; typedef int INT; #define int long long inline void read(int &x){ x=0; int f=1; char ...

额典型的矩阵快速幂 所以毒瘤在每一位都得改变矩阵的权 // luogu-judger-enable-o2 #include<iostream> #include<cstdio> #include<cmath> #include<cstring> #include<algorithm> using namespace std; t...

额这个是卡特兰数 为啥思考选左边一块上面放一块 右边放一堆 这个是卡特兰数递推式 然后我恬不知耻的写的py import math F=[0]*501 F[1]=1 n=int(raw_input()) for i in range(2,n+1): F[i]=F[i-1]*(4*i-2)/(i+1) print(F[n])  ...

、 额这么难得吗 第一道题就是拉格朗日插值法 明显就是多项式插值问题 然后套用拉格朗日插值法就完了 #include<bits/stdc++.h> using namespace std; typedef int INT; #define int long long const int N=101; int x[N]={}; int y[N]={}; int quick...

翻译: 求区间乘积的phi函数值，其中区间中的数的质因数只有60种。带单点修改操作。 明显质因数比较少，勉强可以压一压开一个int64就压住了。 已知   所以易以60的常数求出答案 #include<iostream> #include<cstdio> #include<cmath> #include<cstring> #include&...

矩阵快速幂太暴力了 给出以下证明： 稍微有常识的高考生都会解决这个问题 设  且  解得  设 有  即  根据给出的递推式知 由上式知： 不难看出令     和最开始的式子是相通的 所以就可以快速求出 当然注意特判一下a=1原因我想大家都懂 然后就完了 #include<iostream> #include<cstdio> #incl...

结论题 你可以画一下图 然后发现第K年 可能性为 2^k*n-2^+1-2^k*n 然后快速幂 但是CF有许多神仙数据 第一你不能读进来取mod 第二枚举边界 #include<bits/stdc++.h> using namespace std; const int mod=1e9+7; typedef int INT; #define int long lo...

裴蜀定理到目前为止做的第三个裴蜀定理（HAOI2011向量）这么来看NOIP当时考也不错因为裴蜀定理还是算常见问题转化n个数中选k个求最大的GCDO（n*sqrt sum）完#include<iostream> #include<cstdio> #include<algorithm> #include<cmath> using namespace s...

和礼物差不多但是题目中有一句话：去除零于是在阶乘的时候要搞事情对逆元做阶乘然后扩展Lucas水过#include<bits/stdc++.h> using namespace std; typedef int INT; #define int long long int n,m,k; int mod; int f[11]={}; int fac2[513]={}; int fac5[...

有n个正整数a[i]，设它们乘积为p，你可以给p乘上一个正整数q，使p*q刚好为正整数m的阶乘，求m的最小值。 额简单的数学（我是真的爱数学） 二分答案很明显合法 怎么判断：对每个质数做向下取整除法： #include<iostream> #include<cstdio> #include<cstring> #include<cmath&gt...