![](https://img-blog.csdnimg.cn/20201014180756757.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数论
文章平均质量分 78
wxfwxf328
这个作者很懒,什么都没留下…
展开
-
hdu 3306 Another kind of Fibonacci
As we all known , the Fibonacci series : F(0) = 1, F(1) = 1, F(N) = F(N - 1) + F(N - 2) (N >= 2).Now we define another kind of Fibonacci : A(0) = 1 , A(1) = 1 , A(N) = X * A(N - 1) + Y * A(N - 2) (N >原创 2012-01-21 19:20:32 · 561 阅读 · 0 评论 -
poj 1365 Prime Land 【简单数论】
#include#include#include#includeusing namespace std;bool not_prime[40000];vectorprime,ans;int main(){ for(int i=2;i<40000;i++) if(!not_prime[i]) { prime.push_ba原创 2012-02-24 20:22:25 · 850 阅读 · 0 评论 -
求质数表 素数筛求素数 效率比较 数论相关
第一个是平时最常用的基本写法(还有一种写法是for(j=i*2;j第二个是以i*i为起点开始筛的,要注意浮点溢出!第三个不太好简洁地说清楚程序中途分别输出每次计算了的prime[500000],从而为了验证计算结果的正确性的#include#includeusing namespace std;const int N= 40000000;int prime[N],cnt原创 2012-01-21 04:17:15 · 684 阅读 · 0 评论 -
HDU 3589 Jacobi symbol 【二次剩余】
#include#includeusing namespace std;const int maxn=1000000+10;typedef long long ll;bool not_pri[maxn];int cnt=0,pri[maxn>>3];void init(){ for(int i=2;i<maxn;i++){ if(not_pri[i])原创 2012-03-02 18:07:51 · 909 阅读 · 0 评论 -
HDU 3579 Hello Kiki【中国剩余】
#includeusing namespace std;typedef long long ll;ll ext_gcd(ll a,ll b,ll &x,ll &y){ if(b==0){x=1;y=0;return a;} ll d=ext_gcd(b,a%b,x,y),t; t=x;x=y;y=t-a/b*y; return d;}inline ll原创 2012-02-24 20:25:01 · 926 阅读 · 0 评论 -
hdu 3524 Perfect Squares【打表、除法取余、快速幂】
先打表,找规律#include#includeusing namespace std;int main(){ setmyset; for(long long i=1;i<20;i++,cout<<endl) { myset.clear(); for(long long j=0;j<1000000;j++)原创 2012-02-10 17:24:36 · 1134 阅读 · 0 评论 -
hdu 3547 DIY Cube【polya】
ans=17*x^2+6*x^4+x^8#include#include#includeusing namespace std;const int base=10000;const int width=4;const int N=1300;char s[1000];struct bint{ int ln,v[N]; bint(int r=0) {原创 2012-02-02 20:01:00 · 862 阅读 · 0 评论 -
POJ 2154 Color 【polya+dfs】
看来有必要变得更加强大啊......... 141MS#include#include#includeusing namespace std;const int maxn=1<<15;vector >v;int vis[maxn],prime[maxn];int n,mod,cnt=0;void init(){ for(int i=2;i<maxn;i++)原创 2012-02-07 01:51:47 · 598 阅读 · 0 评论 -
HDU 1695 GCD 数论好题!
http://acm.hdu.edu.cn/showproblem.php?pid=1695此题思路基本分析出来了,可是写错来WA到不认识家了,果断删了自己代码,参考了别人的代码写的:参考自:http://blog.csdn.net/shiren_bod/article/details/5787722#quote(这个初始化写得好比较好了,一开始自己写用了写了两个函数一个求euler+原创 2012-01-25 17:16:16 · 1136 阅读 · 0 评论 -
hdu 3398 String 数论,catalon的一般形式(两次改进,比最初AC快了两倍)【完整版】
好难的题目啊,是cantalon数的一般情况啊,公式为c(a+b,b)-c(a+b,b-1),一开始想从cantalon推这一题,白浪费了一个小时,又是打了m,n的表和组合数的表,才推了这公式写完提交后无限WA,直到练习赛时间过了才想起来是我求的com(a+b, b)-com(a+b,b-1)可能为负值,因为a>b不代表a%p>b%p,即(a-b)%p = ((a%p-b%p)+p)%p而原创 2012-01-21 20:14:56 · 3001 阅读 · 1 评论 -
HDU 3307 Description has only two Sentences 数论
Description has only two SentencesTime Limit: 3000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 199 Accepted Submission(s): 9Problem Descrip原创 2012-01-21 19:59:02 · 1107 阅读 · 0 评论 -
Spoj数论专场解题报告
Spoj数论专题http://acm.hust.edu.cn:8080/judge/contest/view.action?cid=8035#overview Problem ASPOJ ARCTANUse of Function ArctanProblem BSPOJ CUTSQRSCutting off原创 2012-05-29 11:03:52 · 6139 阅读 · 5 评论