HDU6050 Funny Function[矩阵快速幂]

Funny Function

 HDU - 6050 






 

题意:

给出一个函数由题目给的公式可以求出其他项,给出一个n和m,根据这个求出。结果取模1e9+7


题解:

根据前几项

1  1  3  5  11  21  43  85  171  341  683  ......

都是由第二条公式推导出来的

那么当n等于3的时候,第二行 就变成  5  9  19  37  75  149   ......

当n等于4的时候,第二行就变成  10  20  40  80  160  320  640  ......

手动自己写一下几项  会发现同样还是符合第二条 

那么就会有一条

T_T 观察一下吧

然后就会发现一个规律


最后根据这三个式子

我们可以得出

当n为偶数的时候




当n为奇数的时候






然后我们就得到我们要的矩阵了


用矩阵快速幂就可以直接求出答案了



#pragma comment(linker, "/STACK:102400000,102400000")
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<string>
#include<algorithm>
#include<queue>
#include<stack>
#include<set>
#include<map>
#include<vector>
using namespace std;
typedef long long ll;
const int mod=1e9+7;
const int N=2;
ll a[N][N],b[N][N],ans[N][N],tmp[N][N];
ll Mul(ll a[][N],ll b[][N])
{
    memset(tmp,0,sizeof(tmp));
    for (int i=0 ; i<N ; ++i)
        for (int j=0 ; j<N ; ++j)
            if (a[i][j])
                for (int k=0 ; k<N ; ++k)
                    tmp[i][k]=(tmp[i][k]+a[i][j]*b[j][k]%mod+mod)%mod;
    for (int i=0 ; i<N ; ++i)
        for (int j=0 ; j<N ; ++j)
            a[i][j]=tmp[i][j];
}
ll Pow(ll ans[][N],ll a[][N],ll m)
{
    while (m)
    {
        if (m&1)
            Mul(ans,a);
        Mul(a,a);
        m>>=1;
    }
}
ll Add(ll a[][N],ll b[][N],int flag)
{
    for (int i=0 ; i<N ; ++i)
        for (int j=0 ; j<N ; ++j)
            a[i][j]=(a[i][j]+flag*b[i][j]+mod)%mod;
}
void initial(ll n,int flag)
{
    memset(b,0,sizeof(b));
    memset(ans,0,sizeof(ans));
    memset(a,0,sizeof(a));
    for (int i=0 ; i<N ; ++i)
        b[i][i]=ans[i][i]=1;
    a[0][1]=a[1][1]=1;
    a[1][0]=2;
    Pow(ans,a,n);
    for (int i=0 ; i<N ; ++i)
        for (int j=0 ; j<N ; ++j)
            a[i][j]=ans[i][j];
    memset(ans,0,sizeof(ans));
    for (int i=0 ; i<N ; ++i)
        ans[i][i]=1;
    if (flag==1)
    {
        b[0][0]=1;
        b[0][1]=-1;
        b[1][0]=-2;
        b[1][1]=0;
    }
    Add(a,b,flag);
}
int main()
{
    int T;
    scanf("%d",&T);
    while (T--)
    {
        ll n,m;
        scanf("%lld%lld",&n,&m);
        int flag=n&1?1:-1;
        initial(n,flag);
        Pow(ans,a,m-1);
        ll res=0;
        for (int i=0 ; i<N ; ++i)
            res=(res+ans[0][i])%mod;
        printf("%lld\n",res);
    }
    return 0;
}

阅读更多
版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/ControlBear/article/details/76229036
文章标签: HDU 矩阵快速幂
上一篇ZOJ3626 Treasure Hunt Ⅰ[树形DP]
下一篇HDU6030 Happy Necklace[矩阵快速幂]
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页

关闭
关闭
关闭