数学
文章平均质量分 68
INCINCIBLE
这个作者很懒,什么都没留下…
展开
-
扩欧——NKOJ P3677 观光车
问题描述何老板带领n名游客来到一景区大门口,需要乘坐观光车游览景区。景区提供两种观光车,一种是每辆车可以坐a名游客,包一辆车费用是p1块钱;另一种每辆车可以坐b名游客,包一辆车费用是p2块钱。何老板想让这n名游客都坐上观光车,且每辆车都坐满。问何老板至少要花费多少钱?输入格式第一行,一个整数n,表示游客的总数。第二行,两个空格间隔的整数,表示p1和a第三行,两个原创 2016-05-12 18:37:21 · 876 阅读 · 0 评论 -
大组合数取模
考虑C(n,m)%P 情况一:n,m很大,P为素数 处理小范围的阶乘和阶乘的逆元 用卢卡斯定理即可。卢卡斯定理: 情况二: 当P= p1∗p2∗p3∗...∗pn% MathType!MTEF!2!1!+-% feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9w原创 2017-03-01 18:48:28 · 1000 阅读 · 0 评论 -
Miller_Robin素数判定和Pollard_rho质因数分解模板
#include<cstdio>#include<iostream>#include<cstdlib>#include<cmath>#include<cstring>#include<queue>#include<vector>#include<algorithm>#define ULL unsigned long long using namespace std;const int原创 2017-02-18 22:39:02 · 706 阅读 · 0 评论 -
常用数论算法
扩展欧几里得:解不定方程,求逆元int extended_gcd(int a, int b, int &x, int &y) { int r, tmp; if (b==0) { x = 1; y = 0; return a; } r = extended_gcd(b, a % b, x, y); tmp = x; x = y; y = tmp原创 2016-10-22 23:14:53 · 553 阅读 · 0 评论 -
高斯消元(xor方程组)模板
#include#include#include#include#includeusing namespace std;const int inf=0x3f3f3f3f;template inline void _read(T& x){ char ch=getchar();bool sign=true; while(!isdigit(ch)){if(ch=='-')原创 2016-12-14 19:34:38 · 406 阅读 · 0 评论 -
UVa 1642 魔法GCD + 区间问题总结
题目和题解请见刘汝佳紫皮书340页。代码中用了STL中的map。也是现在才知道到map在系统中是以pair的形式存储的,first是键值,second是对应的映射值。map的遍历仍然是使用类的迭代器。代码如下:#include#include#include#include#include#include#define LL long longusing namesp原创 2016-08-10 21:50:55 · 845 阅读 · 0 评论 -
POJ 3150 Cellular Automaton
Cellular AutomatonTime Limit: 12000MS Memory Limit: 65536KTotal Submissions: 3553 Accepted: 1434Case Time Limit: 2000MSDescriptionA cellular automaton is原创 2016-07-21 22:41:39 · 269 阅读 · 0 评论 -
坑爹水题——数列
数列A的满足An=x*An-1+y,给出n,A0,x,y,求模100000007后的结果。输入格式第1行包含4个整数n,A0,x,y输出格式只有1行包含1个整数,An模100000007后的结果。样例输入样例输入1:1000 333 1 233样例输入2:100 2 3 3样例输出样例输出1:233333样原创 2016-06-06 18:56:27 · 1076 阅读 · 0 评论 -
【欧拉】nkoj3685--8数
P3685 8数时间限制 : - MS 空间限制 : 65536 KB 评测说明 : 时限:1000ms问题描述给一个正整数N,问是否存在N的倍数M,且M的各个位全部由数字8组成,如果存在多个取最小的M并输出M由几个8组成。输入格式一行,一个整数N输出格式一行一个整数,表示所求的结果。如果无解,输出0原创 2016-05-22 17:14:34 · 838 阅读 · 0 评论 -
高精度模板
by spark:#include#include#include#include#define LL long long using namespace std;struct sparkint{ static const int BASE=100000000; static const int width=8; vector s; sparkint (LL num=0原创 2016-05-19 18:43:05 · 367 阅读 · 0 评论 -
欧拉 uva 11426 - GCD - Extreme (II)
题意:(详见蓝皮书P125)给定正整数n,求下列表达式的值:G =ij∑ ∑ GCD(i, j)i=1 j=i+1InputThe input file contains at most 100 lines of inputs. Each line contains an integer N (1 < N < 4000001).The meaning of原创 2016-05-21 00:03:38 · 383 阅读 · 0 评论 -
高精度catalan数模板
//高精度的Catalan数模板 公式: h ( i ) = h ( i - 1 ) * ( 4 * i - 2 ) / (i + 1)const int MAX=100;const int BASE=10000;int h[101][MAX],i,j,n;void multiply(int *h,int Max,int b) //大数乘法 求”高精度数a“乘以”低精度数b”原创 2016-05-16 19:27:55 · 476 阅读 · 0 评论 -
回文--nkoj3681
P3681 回文时间限制 : - MS 空间限制 : 65536 KB 评测说明 : 时限1000ms问题描述 给你一个由N个小写字母构成的字符串,请你将它改成“回文串”。你可以任意调整串中字符的位置,但不可以删掉其中的字符。问,最多能得到多少个不同的回文串? 结果可能很大,mod 1,000,000,007 后再输出!原创 2016-05-15 16:30:15 · 506 阅读 · 0 评论 -
费马小定理降幂--nkoj3687 整数拆分
P3687 整数拆分时间限制 : - MS 空间限制 : 65536 KB 评测说明 : 时限1000ms问题描述给你一个正整数N,F(x)表示把N拆分成x个正整数之和的方案数。例如,当n=5时:F(1)=1,方案为:{5}F(2)=4,方案为:{1+4} {4+1} {2+3} {3+2}F(3)=6,方案为:{1+1+3} {1+3+原创 2016-05-26 18:51:29 · 655 阅读 · 0 评论 -
[ZJOI2015] 地震后的幻想乡
P3255【ZJOI2015 Day1】地震后的幻想乡时间限制 : 20000 MS 空间限制 : 262144 KB问题描述 傲娇少女幽香是一个很萌很萌的妹子,而且她非常非常地有爱心,很喜欢为幻想乡的人们做一些自己力所能及的事情来帮助他们。 这不,幻想乡突然发生了地震,所有的道路都崩塌原创 2017-04-25 10:27:44 · 744 阅读 · 0 评论