仍然是学弟出的题目的原题@lher 学弟将题目改成了多组数据,n在ll范围内,所以我就只讲提高版的做法。
链接:https://www.luogu.org/problem/show?pid=2233
题意分析:看题目:)
解题思路:显然对于n为奇数的情况不存在任意路线。接下来我们进行观察数据,显然这题是要递推的。接下来通过暴力打表加手算,我们推出了这个公式:
f[i]=4*f[i-1]-2*f[i-2],f[1]=2,f[2]=8,然后构造对应矩阵进行矩阵快速幂即可得到答案。时间效率\( O ( 2^3 \lg n) \)。
附AC代码:
#include<stdio.h>
#include<iostream>
#include<cmath>
#include<cstring>
#include<string>
#include<algorithm>
#define ll long long
#define mod 10007
using namespace std;
struct zxy{
ll a[][];
}jichu,a,b;
ll n;
inline ll in(){
ll x=,f=;
char ch=getchar();
while(ch<''||ch>'') {if(ch=='-') f=-; ch=getchar();}
while(ch>=''&&ch<='') x=x*+ch-'',ch=getchar();
return x*f;
}
inline zxy mult(zxy a,zxy b){
zxy c;
memset(c.a,,sizeof(c.a));
for (int i=; i<=; ++i)
for (int j=; j<=; ++j)
for (int k=; k<=; ++k)
c.a[i][j]+=(a.a[i][k]*b.a[k][j]+*mod)%mod,c.a[i][j]=(c.a[i][j]+*mod)%mod;
return c;
}
inline zxy ksm(zxy a,ll k){
if (k==) return jichu;
if (k==) return a;
zxy tt=ksm(a,k>>);
if (k&) return mult(a,mult(tt,tt));
return mult(tt,tt);
}
int main(){
int T=in();
memset(jichu.a,,sizeof(jichu.a));
memset(b.a,,sizeof(b.a));
memset(a.a,,sizeof(a.a));
jichu.a[][]=jichu.a[][]=;
b.a[][]=;b.a[][]=;
a.a[][]=;a.a[][]=-+mod,a.a[][]=;
while(T--){
n=in();
if (n&) printf("0\n");
else{
n=n/-;
if (n<) printf("%lld\n",b.a[n][]);
else{ zxy ans=mult(ksm(a,n-),b);printf("%lld\n",ans.a[][]);}
}
}
return ;
}
本文由Melacau编写,Melacau代表M星向您问好,如果您不是在我的博客http://www.cnblogs.com/Melacau上看到本文,请您向我联系,email:13960948839@163.com.