10655 - Contemplation! Algebra

最新推荐文章于 2019-07-06 14:08:00 发布

*ACoder*

最新推荐文章于 2019-07-06 14:08:00 发布

阅读量218

点赞数

分类专栏： # 线性代数

本文链接：https://blog.csdn.net/FSAHFGSADHSAKNDAS/article/details/89007445

版权

线性代数专栏收录该内容

15 篇文章 0 订阅

订阅专栏

这篇博客介绍了如何解决UVA在线判题平台上的10655号问题——Contemplation!，通过解析数学公式fn=an+bn，并利用递推关系fn+2=−qfn+pfn+1，博主探讨了使用矩阵快速幂的方法来求解问题的策略，并提供了相应的代码实现。

摘要由CSDN通过智能技术生成

链接

https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=1596

题解

令 $f_n=a^n+b^n$
则 $f_{n+2}=-qf_{n}+pf_{n+1}$
矩阵快速幂

代码

#include <bits/stdc++.h>
#define cl(x) memset(x,0,sizeof(x))
using namespace std;
typedef long long ll;
struct matrix
{
    ll a[5][5], n, m;
    matrix(ll n, ll m):n(n),m(m){cl(a);}
    ll* operator[](int index)
    {
        return a[index];
    }
    matrix operator*(matrix b)
    {
        ll i, j, k;
        matrix c(n,b.m);
        for(i=1;i<=c.n;i++)for(k=1;k<=m;k++)for(j=1;j<=c.m;j++)c[i][j]+=a[i][k]*b[k][j];
        return c;
    }
    matrix fastpow(ll b)
    {
        matrix ans(n,n), t(*this);
        ll i;
        for(i=1;i<=n;i++)ans[i][i]=1;
        for(;b;b>>=1,t=t*t)if(b&1)ans=ans*t;
        return ans;
    }
};
int main()
{
    ll i, p, q, n;
    while(scanf("%lld %lld %lld",&p,&q,&n)==3)
    {
        matrix f(1,2), trans(2,2);
        f[1][1]=2, f[1][2]=p;
        trans[1][1]=0, trans[1][2]=-q;
        trans[2][1]=1, trans[2][2]=p;
        f=f*trans.fastpow(n);
        printf("%lld\n",f[1][1]);
    }
    return 0;
}