![](https://img-blog.csdnimg.cn/20201014180756780.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数论
文章平均质量分 80
夨落旳尐孩
我要这天,再遮不住我眼,要这地,再埋不了我心,要这众生,都明白我意,要那诸佛,都烟消云散。
展开
-
poj 数论 Semi-prime H-numbers
题意:H-number是4*n+1这样的数,如1,5,9,13... 。H-primes是这样一个H-number:它只能唯一分解成1*它本身,而不能表示为其他两个H-number的乘积。一个H-semi-prime是一个这样的H-number:它正好能表示成两个H-primes的乘积(除了1*它本身),这种表示法可以不唯一,但它不能表示为3个或者以上H-primes的乘积。现在给出一个数n,要求原创 2013-09-11 14:08:56 · 845 阅读 · 0 评论 -
poj 2635 The Embarrassed Crypto…
The Embarrassed Cryptographer Time Limit: 2000MS Memory Limit: 65536K Total Submissions: 9405 Accepted: 2433 Description 2635 The Embarrassed Cryptographer" TITLE="poj 2635 The Embarras原创 2013-09-11 14:08:53 · 965 阅读 · 0 评论 -
poj 数论 1845
思路: 将a进行质因数分解,则a^b=(p1^k1)*(p2^k2)*...*(pn^kn) 根据公式,所有因子之和为 (1+p1+p1^2+...p1^k1)*(1+p2+p2^2+...p2^k2)*...*(1+pn+pn^2+...+pn^kn) 计算1+p+p^2+...p^n可以利用二分进行加速 当n为奇数时,例如n=5 则1+p+p^2+p^3+p^4+p^5=(1+p+p原创 2013-09-11 14:08:58 · 883 阅读 · 0 评论 -
poj 1061 青蛙的约会
此题其实就是扩展欧几里德算法-求解不定方程,线性同余方程。 设过s步后两青蛙相遇,则必满足以下等式: (x+m*s)-(y+n*s)=k*l(k=0,1,2....) 稍微变一下形得: (n-m)*s+k*l=x-y 令n-m=a,k=b,x-y=c,即 a*s+b*l=c 只要上式存在整数解,则两青蛙能相遇,否则不能。 求a * x + b * y = n的整原创 2013-09-11 14:09:15 · 869 阅读 · 0 评论 -
poj 2891
假设C ≡ A1 (mod B1),C ≡ A2 (mod B2)。令C = A1 + X1B,那么X1B1 ≡ A2 − A1 (mod B2)。用扩展欧几里德算法求出X1,也就求出C。令B = lcm(B1, B2),那么上面两条方程就可以被C’ ≡ C (mod B)代替。迭代直到只剩下一条方程。 #include using namespace std; __int64 egcd(__i原创 2013-09-11 14:09:17 · 868 阅读 · 0 评论