ACM[数学基础]
sevenster
= = 水男
展开
-
POJ 2313 中位数的应用
我发誓再也不看Discuss了~~ #include #include #include using namespace std; int A[101]; int B[101]; int Mid( int a,int b,int c ) { int num[4];原创 2011-09-24 16:35:57 · 868 阅读 · 0 评论 -
POJ 3233 快速矩阵乘法
觉得专业课学习有点落下了,而且acm这方面也没有弄好= =蛋疼死了!!稍微写写博客就开始学习吧!! 矩阵的快速幂和a^b%n类似。运用二分的思想。在这题中需要计算的是A+A^2+A^3+...+A^k的和,由于矩阵相乘有结合律,所以!@#$%^&*..... 不多说了,代码很丑。还是拿来占个位置吧.... 第二次写递归啊!!我的神啊~~ #include #define MAXN 31 us原创 2011-10-20 08:08:41 · 820 阅读 · 1 评论 -
POJ 2118 Matrix线性递推
第一次做这种少人做的题目,感觉好爽啊~~~ 这个题目的解法借鉴了Matrix67大神的好文章: 10种经典的矩阵解法 这里讲得很好,我用的是矩阵的经典7,专门用来借线性递推的。 再结合矩阵的快速幂,这题无疑是小细节要注意注意了! #include #define MAXN 101 #define mod 10000 #define ull unsigned long long using原创 2011-10-20 18:10:28 · 872 阅读 · 0 评论 -
k次方的首末三位
K次方 Time Limit : 2000/1000ms (Java/Other) Memory Limit : 65535/32768K (Java/Other) Total Submission(s) : 10 Accepted Submission(s) : 3 Font: Times New Roman | Verdana | Georgia Font Si原创 2012-08-24 18:34:22 · 2503 阅读 · 0 评论 -
约瑟夫环 数学解
题目大意: 有给定的N个人,主角编号13,求最小的报数序号M,使得主角为赢家。 报数序列为1,M%N这样.. 1是第一个出去的 网上看的数学方法,各种感觉不爽! 觉得算法中最神奇的就是2分和递归了。美妙! 数学分析过程网上都有,不解释了。 #include using namespace std; int find( int n,int m ) { int s=0; for(原创 2012-08-25 12:52:49 · 944 阅读 · 0 评论 -
POJ 2689 数论
刚做素数部分,把以前不懂的题目搞定了。还不错。 #include #define N 50001 using namespace std; bool f[N],seg[1111111]; int prime[N],pcnt; __int64 rec[1111111]; void deal( __int64 L,__int64 U ) { memset( seg,true,sizeof(原创 2012-09-19 11:41:59 · 1167 阅读 · 0 评论 -
nefu 120 梅森素数 Lucas-Lehmer判定法
定义:如果m是一个正整数,且2^m-1是一个素数,则m必是素数.反之,如果m是一个正整数,素数且Mm=2^m-1成为第m个梅森数;如果p是一个素数,并且Mp=2^p-1也是素数,那么Mp=2^p-1也是素数,那么Mp就称为梅森素数. Lucas-Lehmer判定方法. 设p是素数,第p个梅森数为Mp=2^p-1,r1=4; 对于k>=2,利用rk=((r(k-1))^2-1)%Mp,0 可原创 2012-09-19 16:40:40 · 2384 阅读 · 0 评论 -
poj 3641 Robin-Miller 素性测试
题意: a^p=a%p;则称p为以a为基的素数。给定p与a,求判断在p不为素数的情况下,a^p=a%p是否成立。 #include #define ll long long using namespace std; ll mul( ll a,ll b,ll mod ) { ll ret=0; while( b>0 ) { if( b&1 )原创 2012-09-19 16:43:36 · 1024 阅读 · 0 评论