斐波那契数列（升级版）

#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<algorithm>
using namespace std;
int read()
{
    int x=0,f=1;
    char ch=getchar();
    while(ch<‘0‘||ch>‘9‘)
    {
        if(ch==‘-‘) f=-1;
        ch=getchar();
    }
    while(ch>=‘0‘&&ch<=‘9‘)
    {
        x=x*10+ch-‘0‘;
        ch=getchar();
    }
    return f*x;
}//读入优化 
int main()
{
    int p=pow(2,31);
    int n=read();
    double x=sqrt(5.0);
    long long s=(pow(((1+x)/2),n)/x-pow(((1-x)/2),n)/x);//斐波那契数列通项公式 
    s=s%p;//取模运算 
    printf("%lld=",s);
    long long sum=0;
    int m=2;
    while(s!=1)
    {
        if(s%m)
         m++;
        else
        {
            sum++;
            if(sum==1)
             printf("%d",m);
            else 
              printf("*%d",m);
            s/=m;
        }
    }//分解质因数 
    return 0;
}