UVA10655-Contemplation! Algebra

题目分析：有种错误的方法是这样的：利用已知的两个方程联立，求解出a和b，进而求出答案。这种方法之所以错，是因为这种方法有局限性。联立之后会得到一个二元一次方程，只有当该方程有实数解确切的说是当某个数据满足该方程有实数解时，这种方法得到的结果才有可能正确。显然，题中数据不可能这么片面。正确的方法是这样的：

令a+b=A，ab=B，S(n)=aⁿ+bⁿ。则S(n)=aⁿ+bⁿ=(a+b)(a^n-1+b^n-1)-ab^n-1-a^n-1b=(a+b)(a^n-1+b^n-1)-ab(a^n-2+b^n-2)=A*S(n-1)-B*S(n-2)（n≥2)。

n=1或n=0时，答案显然。

要注意：数据中可能会出现A=B=0的情况，此时a=b=0，所以不能简单的判定当A=B=0时就结束输入数据。

所以输入要用c++，或者while(scanf("%d %d %d",&x,&y,&n)==3)

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <cstring>

using namespace std;

#define ll long long

struct Matrix
{
    ll v[3][3];
    Matrix()
    {
        memset(v,0,sizeof v);
    }
}m,k,t;

Matrix mul(Matrix a,Matrix b)
{
    Matrix sum;
    for(int i=1; i<3; i++)
        for(int j=1; j<3; j++)
            for(int k=1; k<3; k++)
                sum.v[i][j]+=a.v[i][k]*b.v[k][j];
    return sum;
}

Matrix pow(Matrix a,ll n)
{
    Matrix sum=k;
    while(n)
    {
        if(n&1) sum=mul(sum,a);
        n>>=1;
        a=mul(a,a);
    }
    return sum;
}

int main()
{
    ll x,y,n;
    while(cin>>x>>y>>n)
    {
        if(n==0) {printf("2\n");continue;}
        if(n==1) {printf("%lld\n",x);continue;}
        k.v[1][1]=k.v[2][2]=1;
        m.v[1][1]=x;m.v[1][2]=1;m.v[2][1]=-y;
        t=pow(m,n-1);
        printf("%lld\n",x*t.v[1][1]+2*t.v[2][1]);
    }
    return 0;
}