Combination bzoj-2982 题目大意：求$C_n^m/%10007$。 注释：$1\le n,m\le 2\cdot 10^9$。 想法：裸卢卡斯定理。 先处理出$mod$数之内的阶乘和阶乘的逆元。 然后用$Lucas$直接算即可。 Code： 小结：刷水有益身心健康。 ...
分类：
其他好文 时间：
2018-12-10 11:47:27
阅读次数：
167
const int maxn=3000000; ll inv[maxn+5]; void marktable_inv(int p){ inv[0]=inv[1]=1; for(int i=2;i<maxn;i++) inv[i]=(p-p/i)*inv[p%i]%p; } ...
分类：
其他好文 时间：
2018-12-08 19:14:54
阅读次数：
125
1 利用扩展欧几里得求解 2 利用欧拉函数 ...
分类：
其他好文 时间：
2018-12-08 18:27:46
阅读次数：
159
题意： 求n行三角形中等边三角形个数，图二的三角形也算。 n<=1e9 思路： 打表找下规律，打表方法：把所有点扔坐标系里n^3爆搜即可 打出来为 1，5，15，35，70，126，210.. 没感觉，作差 4, 10, 20, 35, 56, 84 还是没感觉，作差 6, 10, 15, 21, ...
分类：
其他好文 时间：
2018-11-29 22:28:50
阅读次数：
1141
这是我做题史上摔得最惨的一道黄题，15条记录转眼化为泪水。o(╥﹏╥)o 这道题目从10.12开始尝试，随机跳题跳到了这题，一看就是卡特兰数，因为样例太像了。。 然后小心证明 $dp[i]=\sum{dp[j]*dp[i-j-2]}$ 这个就是catalan的式子啊，， 然后打了一发阶乘+逆元求组合 ...
分类：
其他好文 时间：
2018-11-29 14:56:32
阅读次数：
207
逆元定义先摆上来 对于正整数a和m，如果a*x≡1(mod m)，那么把这个同余方程中x的最小正整数解叫做a模m的逆元。 求解方法： 1.扩展欧几里得 利用欧几里得求x 先将方程转化为 ax-my=1 此时求解x和y 最后利用返回的gcd(a,m)==1 如果成立，则x为逆元存在，否则不存在 注意最 ...
分类：
其他好文 时间：
2018-11-19 18:11:50
阅读次数：
703
#include using namespace std; typedef long long ll; void exgcd(ll a,ll b,ll& d,ll& x,ll& y) { if(!b) { d = a; x = 1; y = 0; } else{ exgcd(b, a%b, d, y... ...
分类：
其他好文 时间：
2018-11-07 21:39:05
阅读次数：
176
数论基本糙作： gcd，快速幂，逆元，欧拉函数，分解因数balabala一通乱搞。 POJ1845 Sumdiv （数论：算数基本定理+数论基本操作） 题目： Sumdiv Time Limit: 1000MS Memory Limit: 30000K Total Submissions: 2901 ...
分类：
其他好文 时间：
2018-11-07 15:27:50
阅读次数：
145
前言 逆元是一直以来想补的坑了 求逆元 即求 x 满足 a/b % mod =a*x%mod =1%mod 费马小定理 限制：mod为质数 exgcd 限制：b和mod互质 数学变形 限制：b能被a整除 公式递推 限制：mod是质数 逆元不存在的时候会输出0 ...
分类：
其他好文 时间：
2018-11-05 19:20:57
阅读次数：
786
//法1：费马小定理求逆元，p需为质数 : typedef long long ll; const int MOD = 9973; ll PowerMod(ll a, ll b, ll p) { //(a^b)%p ll ans = 1; a = a%p; while (b>0) { if (b &... ...
分类：
其他好文 时间：
2018-11-01 15:05:17
阅读次数：
169
