UVA11653_Buses

这个题目很有意思，一不小心就会让人坑在里面。

题意是这样的，给你n，k，l。分别表示总共的长度，长度为5和10的车的不同颜色数量现在问你要把n的填满有多少种方案。

很多人一开始都会脑子一根筋地想用排列组合去搞这个题目。然而实际上不是这样的。因为排列组合计算量巨大，而且这个题目的数据范围是10^15，绝对无法承受。

其实我们可以先把n/5，这样相当于是放长度为1和长度为2的方案了。

我们加入一个状态量f[i]，其意义为长度为i的排列方案有多少种？

那么我们可以迅速地得出这个状态转移的递推式：f[i]=k*f[i-1]+l*f[i-2]。（分别表示放长度为1和2的情况嘛）

这样你是否有些眼熟了。。。。 对没有错，就是它——矩阵快速幂。

这样我们要求解的范围十分之小，瞬间变为了log级别，所以答案最终就可以轻松飘过啦。

#include <iostream>
#include <cstdio>
#include <cstring>
#define M 1000000
#define ll long long
using namespace std;

struct mat{
    ll a[][];
    void init(ll k,ll l)
    {
        a[][]=k,a[][]=l,a[][]=,a[][]=;
    }
    void E()
    {
        a[][]=,a[][]=,a[][]=,a[][]=;
    }
}tep,ans;

mat mul(mat a1,mat a2)
{
    mat f;
    memset(f.a,,sizeof f.a);
    for (ll i=; i<; i++)
        for (ll j=; j<; j++)
            for (ll k=; k<; k++)
                f.a[i][j]=(f.a[i][j]+a1.a[i][k]*a2.a[k][j])%M;
    return f;
}

mat power(mat cur,ll y)
{
    mat now;
    now.E();
    while (y)
    {
        if (y&) now=mul(now,cur);
        y>>=;
        cur=mul(cur,cur);
    }
    return now;
}

int main()
{
    ll n,k,l,f1,f2;
    while (scanf("%lld%lld%lld",&n,&k,&l)!=EOF)
    {
        n/=;
        k%=M,l%=M;
        f1=k,f2=(k*k+l)%M;
        if (n==)
        {
            printf("%06lld\n",f1);
            continue;
        }
        if (n==)
        {
            printf("%06lld\n",f2);
            continue;
        }
        tep.init(k,l);
        ans=power(tep,n-);
        printf("%06lld\n",(ans.a[][]*f2+ans.a[][]*f1)%M);
    }
    return ;
}