Hdu 2256 Problem of Precision

Problem of Precision

Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 614    Accepted Submission(s): 344


Problem Description

 

Input
The first line of input gives the number of cases, T. T test cases follow, each on a separate line. Each test case contains one positive integer n. (1 <= n <= 10^9)
 

Output
For each input case, you should output the answer in one line.
 

Sample Input
  
  
3 1 2 5
 

Sample Output
  
  
9 97 841
 

Source
 

Recommend
lcy


参考别人的解题思路,写的代码,感觉他的解题思路写的很清晰。下面便是他的解题思路:
解题思路原文
这个题目用快速幂+fmod是AC不了的,应该是浮点数精度的原因,推导的过程都体现在图上了,剩下的工作就是二分矩阵求解了。

  此外,在计算完x[n]和y[n]之后不能直接用(x[n]+(int)(y[n]*sqrt(6.0)))24来得到最后的结果的,先取整再模和先模再取整的结果是不一样的这一点举个例子就比较容易看出来了。

   比如(2000*1.372)00,这样先乘再取模得到的结果就是744,如果先把2000模1000,显然最后结果就是0了。

#include<stdio.h>

int Matrix[3][3];
int tmp1[3][3],tmp2[3][3];

void fun(int n)                         //求矩阵的快速幂
{
    int i,j,k;
    n-=2;
    Matrix[1][1]=tmp1[1][1]=5;
    Matrix[1][2]=tmp1[1][2]=12;
    Matrix[2][1]=tmp1[2][1]=2;
    Matrix[2][2]=tmp1[2][2]=5;
    while(n)
    {
        if(n%2)
        {
            for(i=1;i<=2;i++)
                for(j=1;j<=2;j++)
                    tmp2[i][j]=Matrix[i][j];
            for(i=1;i<=2;i++)
                for(j=1;j<=2;j++)
                {
                    Matrix[i][j]=0;
                    for(k=1;k<=2;k++)
                        Matrix[i][j]=(Matrix[i][j]+tmp1[i][k]*tmp2[k][j])24;
                }
        }
        n/=2;
        for(i=1;i<=2;i++)
            for(j=1;j<=2;j++)
                tmp2[i][j]=tmp1[i][j];
        for(i=1;i<=2;i++)
            for(j=1;j<=2;j++)
            {
                tmp1[i][j]=0;
                for(k=1;k<=2;k++)
                    tmp1[i][j]=(tmp1[i][j]+tmp2[i][k]*tmp2[k][j])24;
            }
    }
}

int main(void)
{
    int T,n;
//    freopen("d:\\in.txt","r",stdin);
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n);
        fun(n);
        if(n==1)                                         //n==1属于特殊情况
            printf("9\n");
        else
            printf("%d\n",(2 * (Matrix[1][1]*5+Matrix[1][2]*2 )-1)24);

    }
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值