注册 登录
编程论坛 数据结构与算法

uva 10236 斐波拉契素数问题

brian1994 发布于 2011-05-16 13:55, 1681 次点击
斐波那契素数
时限:8s
【题目描述】
斐波那契数列为1,1,2,3,5,8,13……,你可以看得出,除了第一第二个数以外,每个数都是前两个数的和。现在定义,斐波那契素数是指与其他比它小的斐波那契数都互质的斐波那契数,为2、3、5、13、89、233、1597、4181、28657、514229……
现在让你求出第i个斐波那契素数的前九为,若不足九位就全部输出;
【输入】
    输入共n行,每行一个整数i(1≤i≤22000);
【输出】
    输出共n行,每行按要求输出第i个斐波那契素数;
【输入样例】
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
【输出样例】
2
3
5
13
89
233
1597
4181
28657
514229
1346269
24157817
165580141
433494437
297121507
533162911
956722026
250473078
449455702
308061521
//----------------------------------------------
个人想法:用高精度计算斐波拉契数列,素数位置的斐波拉契数就是结果,试过前面几项可以,后面不知道行不行。
          但这种方法无疑是非常繁的,各位大牛有什么想法????
17 回复
#2
buffer2011-05-16 20:54
这个对你会有帮助:http://en.
GCD(F_n, F_m) = F_GCD(n, m)
Fibonacci_Prime[k] = Fibonacci[Prime[k]], k >= 3
再根据Fibonacci的通项公式(就是用黄金分割率表示的那个,这儿贴公式不方便,不写了~),就可以求出前9位数字了。
#3
fangdong652011-05-16 20:56
可以试试用字符表示数字
#4
brian19942011-05-17 13:31
回复 2楼 buffer
为什么没有翻译... ...还有,这个贴为什么被设为高亮????
#5
buffer2011-05-17 14:11
回复 4楼 brian1994
玩ACM的还怕这点儿英文吗
#6
brian19942011-05-17 14:14
回复 楼主 brian1994
推敲过了,我的想法是错的。
那个“素数位的斐波那契数是素数”是错的
#7
brian19942011-05-17 14:15
回复 5楼 buffer
我只是个高中生
#8
brian19942011-05-17 16:00
要结贴了!
通过:GCD(F_n, F_m) = F_GCD(n, m)
可以证明:素数位的斐波那契数是斐波那契素数!
稍后间上传代码0...0
#9
buffer2011-05-17 17:49
回复 8楼 brian1994
嗯,精度是这题最容易出错的原因~
试试这几组数据:
233
6879
15086
15520
输出:
117851144
842347072
920759405
386616082
#10
寒风中的细雨2011-05-17 19:43
回复 4楼 brian1994
置顶有什么问题   要理解  难得使用下版主身份 ~~~~~
#11
brian19942011-05-18 13:22
回复 9楼 buffer
我第15个就精度问题了...
#12
brian19942011-05-18 13:47
OK!!!!!!!!
#include<cstdio>
#include<cstring>
#include<cmath>
const int num=249439;
int prime[22001];   //素数
bool boprime[num+1];
double f[num+1];       //斐波那契数列
double goldnum;     //goldnum黄金分割比
int n;   

void find()
{
     //fill
     memset(boprime,0,sizeof(boprime));
     //find prime
     prime[0]=2; prime[1]=3; prime[2]=4;
     int i=2,j;
     while (i<=num)
     {
           if (!boprime[i])
           {
              j=2*i;
              while (j<=num)
              {
                    boprime[j]=true;
                    j+=i;
              }
              if (i!=2 && i!=3)
              prime[++prime[0]]=i;
           }
           i++;
     }
}

void doit()
{
     goldnum=(sqrt(5)-1)/2;
     f[1]=f[2]=1;
     for (int i=3;i<=83;i++) f[i]=f[i-1]+f[i-2];
     for (int i=40;i<=83;i++) while (f[i]>=1000000000) f[i]/=10;
     for (int i=84;i<=num;i++)
     {
         f[i]=f[i-1]/goldnum;
         while (f[i]>=1000000000) f[i]/=10;
     }
}

void init()
{
     while (scanf("%d",&n)==1)
       printf("%.lf\n",floor(f[prime[n]]));
}

int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    find();
    doit();
    init();
    return 0;
}
#13
brian19942011-05-19 14:01
我的程序过不了!!!!
#14
brian19942011-05-19 16:42
回复 9楼 buffer
看看我的程序如何解决精度问题!!!
本以为可以AC,谁知是WA
#15
brian19942011-05-20 15:56
如果高精度怎么处理相加要错位这个问题??
#16
brian19942011-05-20 16:51
好惨啊!弄了我这么久...
还好最后AC了
#include<cstdio>
int n;
int tot=0;
int ans[25000];

int main()
{
    long long t0,t1=1,t2=1;
    for (int i=3;tot<22000;i++)
    {
        t0=t1+t2;
        t2=t1;
        t1=t0;
        if (t0>=1e18)
        {
           t0/=10;
           t1/=10;
           t2/=10;
        }
        bool p=true;
        for (int j=2;j*j<=i;j++)
            if (i%j==0)
            {
               p=false;
               break;
            }
        if (p || i==4)
        {
           long long p=t0;
           while (p>=1e9) p/=10;
           ans[++tot]=p;
        }
    }
    while (scanf("%d",&n)==1) printf("%d\n",ans[n]);
    return 0;
}
#17
jiava2011-05-26 21:14
java皇冠高级群:7156436
java皇冠高级群:93787804
(注:群名原创,有特定群标,其他有同名或类似名称如java静区,或相同群标的群一律为"盗版"群)
群简介:快乐学习,真诚交友!以JAVA为主,面向C/C++,C#,Android,iphone,jsp,Sql,Oracle,php,web,.net..
欢迎工作者,爱好者加入。
#18
schnee2011-09-20 20:03
楼主牛逼
我TLE了好久

[ 本帖最后由 schnee 于 2011-9-20 20:05 编辑 ]
1