随笔分类 -  数论+大整数

hdu 2669 扩展欧几里得（裸）
摘要：#include#include#define ll __int64ll gcd(ll a,ll b,ll &x,ll &y) { ll k; if(b==0) { x=1;y=0; return a; } ll d=gcd(b,a%b,x,y); k=x; x=... 阅读全文
posted @ 2015-03-31 19:20 HYDhyd 阅读(151) 评论(0) 推荐(0)
利用牛顿迭代法求平方根
摘要：求n的平方根，先假设一猜测值X0 = 1，然后根据以下公式求出X1，再将X1代入公式右边，继续求出X2…通过有效次迭代后即可求出n的平方根，Xk+1先让我们来验证下这个巧妙的方法准确性，来算下2的平方根 (Computed by Mathomatic)1-> x_new = ( x_old + y/... 阅读全文
posted @ 2014-12-27 21:46 HYDhyd 阅读(944) 评论(0) 推荐(0)
hdu 5073 推公式相邻质心转换
摘要：#include#include#include#define N 51000int cmp(const void *a,const void *b) {return *(int *)a-*(int *)b;}int a[N];int sum[N];double all[N];double ss(... 阅读全文
posted @ 2014-10-22 18:52 HYDhyd 阅读(105) 评论(0) 推荐(0)
高斯消元法模板
摘要：转载：http://yzmduncan.iteye.com/blog/1740520例：ZOJ3645题意：高斯消元模板题(浮点型)/**高斯消元求解线性方程组.*/#include #include #include #include #include using namespace std;//... 阅读全文
posted @ 2014-10-17 11:16 HYDhyd 阅读(301) 评论(0) 推荐(0)
高斯消元中用到的相关概念和模板
摘要：线性独立：一组向量，如果其中的一个向量能用别的向量的线性组合的形式表示出来，则这组向量是线性相关的；否则向量组就是线性无关的。秩：在线性代数中，一个矩阵A的列秩是A的线性獨立的纵列的极大数目。类似地，行秩是A的线性獨立的横行的极大数目。矩阵的列秩和行秩总是相等的，因此它们可以简单地称作矩阵A的秩。通... 阅读全文
posted @ 2014-10-17 11:09 HYDhyd 阅读(168) 评论(0) 推荐(0)
hdu 2604 矩阵快速幂模板题
摘要：/*矩阵快速幂： 第n个人如果是m，有f(n-1)种合法结果 第n个人如果是f,对于第n-1和n-2个人有四种ff,fm,mf,mm其中合法的只有fm和mm 对于ffm第n-3个人只能是m那么有f(n-4)种 对于fmm那么对于第n... 阅读全文
posted @ 2014-09-23 15:12 HYDhyd 阅读(136) 评论(0) 推荐(0)
hdu 5035 概率论
摘要：HDU 5035 Delivery 题解题解数学概率期望表示今年刚学概率论//https://www.zybuluo.com/rihkddd/note/34286题目大意：Matt一个发快递，快递点有N个职员，职员处理一个客户的时间服从指数分布：f(t)=λe−λt,其中的参数λ为职员的效率，现在给... 阅读全文
posted @ 2014-09-22 11:31 HYDhyd 阅读(199) 评论(0) 推荐(0)
hdu 5015 矩阵快速幂（可用作模板）
摘要：转载：http://blog.csdn.net/wdcjdtc/article/details/39318847之前各种犯傻 推了好久这个东西。。后来灵关一闪 就搞定了。。矩阵的题目，就是构造矩阵比较难想！题意：给出一个矩阵的第一列和第一行（下标从0开始），（0,0）位置为0，第一行为，233,23... 阅读全文
posted @ 2014-09-19 19:36 HYDhyd 阅读(128) 评论(0) 推荐(0)
poj 3006水题打素数表
摘要：#include#include#define N 1100000int isprim[N],prime[N];void isprime() { int i,j;memset(isprim,-1,sizeof(isprim)); isprim[1]=0; for(i=2;i<=1000;... 阅读全文
posted @ 2014-07-19 10:47 HYDhyd 阅读(134) 评论(0) 推荐(0)
hdu 1576扩展欧几里得算法
摘要：#include#define ll long long/* 2.那么x，y的一组解就是x1*m1，y1*m1，但是由于满足方程的解无穷多个，在实际的解题中一般都会去求解x或是y的最小正数的值。以求x为例，又该如何求解呢？还是从方程入手，现在的x，y已经满足a*x+b*y=m，那么a*(x+n*b)... 阅读全文
posted @ 2014-07-07 17:20 HYDhyd 阅读(166) 评论(0) 推荐(0)
欧几里德算法的扩展-求解不定方程
摘要：对于不完全为 0 的非负整数 a，b，gcd（a，b）表示 a，b 的最大公约数，必然存在整数对 x，y ，使得 gcd（a，b）=ax+by。c++语言实现#include#includeusingnamespacestd;intx,y,q;voidextend_Eulid(inta,intb){... 阅读全文
posted @ 2014-07-07 16:14 HYDhyd 阅读(299) 评论(0) 推荐(0)
hdu 1250 简单大整数加法
摘要：#include#include#define N 3100int a[N],b[N],c[N],d[N],e[N];int main() { int n,i,j,k,t; while(scanf("%d",&n)!=EOF) { memset(a,0,sizeof(a));memset... 阅读全文
posted @ 2014-07-05 11:41 HYDhyd 阅读(153) 评论(0) 推荐(0)
hdu 1568关于斐波那契数列的公式及其思维技巧
摘要：先看对数的性质,loga(b^c)=c*loga(b),loga(b*c)=loga(b)+loga(c);假设给出一个数10234432,那么log10(10234432)=log10(1.0234432*10^7)=log10(1.0234432)+7;log10(1.0234432)就是log... 阅读全文
posted @ 2014-07-05 09:10 HYDhyd 阅读(258) 评论(0) 推荐(0)
hdu 1713求分数的最小公倍数
摘要：题意中的圈数和天数说反了#include__int64 gcd(__int64 a,__int64 b) {/*比如4/3 3/5通分20/15 9/15所以这两个分数的最小公倍数为 180/15 （20，9的最小公倍数为180）..然后约分下就好了。。所以答案就是12*/if(b==0) ... 阅读全文
posted @ 2014-06-28 16:14 HYDhyd 阅读(274) 评论(0) 推荐(0)
nyoj 803 大数问题
摘要：#include#include#define ll long long#define N 110000int main() { ll i,sum,b,k,len,f[N]; char s[N],st[3]; while(scanf("%s%s%lld",s,st,&b)!=EOF) {... 阅读全文
posted @ 2014-05-22 09:49 HYDhyd 阅读(109) 评论(0) 推荐(0)
hdu 1163
摘要：别人的代码开始自己不知道什么数论解法： ab*ab=(a*10+b)(a*10+b)=a^2*100+2ab*10+b^2 所以the root digital=(a+b)*(a+b)； 而数论中的定理：两数之积对9取余数等于两数对9的余数的乘积。 ... 阅读全文
posted @ 2014-05-20 09:46 HYDhyd 阅读(150) 评论(0) 推荐(0)
1066N ！最右边非零数
摘要：http://hi.baidu.com/nicker2010/item/4fa83c4c5050b3e5a4c066ec另一个Lastnon-zeroDigitinN!TimeLimit:2000/1000MS(Java/Others)MemoryLimit:65536/32768K(Java/Ot... 阅读全文
posted @ 2014-04-14 15:33 HYDhyd 阅读(393) 评论(0) 推荐(0)
hdu 1787
摘要：#include#includeint huzhi(int n) {//欧拉函数公式int i,j=n;for(i=2;i1)//当n为2时用到此步j=j/n*(n-1);return j;}int main() {int n;while(scanf("%d",&n),n) {printf("%d\... 阅读全文
posted @ 2013-11-20 15:31 HYDhyd 阅读(107) 评论(0) 推荐(0)
同余定理
摘要：http://www.cnblogs.com/zhixingqiezhixing/archive/2012/04/03/2430676.htmlhttp://blog.sina.com.cn/s/blog_6e223f1d0101bkry.html 阅读全文
posted @ 2013-10-03 18:38 HYDhyd 阅读(182) 评论(0) 推荐(0)