序文
(この記事を見て、最初にマトリックスとは何かを理解してください)
マトリックスアクセラレーションは非常に魔法のアルゴリズムであり、極端なデータの下で複雑な問題を短時間で解決でき、非常にシンプルです。例から始めましょう:
例1:フィボナッチアイテムn
フィボナッチ数列は、を満たします。(N <= 1e9)
このデータは大きく、再帰によって日常的に解決することはできません。次に、迂回せずに、式の導出を開始します。
(キーの繰り返し)のため、行列の乗算(n * mおよびm * q行列)に入れることができます。2項×2項であるため、m = 2であり、通常の累積乗算を実現するには、結果もn * 2の行列、つまりq = 2であることを確認する必要があります。他に必要な数がないため、n = 1です。行列は、乗算後、(1,1)の位置はです。ただし、左行列はである必要があるため、(2,2)の位置は、と他の2つの数値の「加算と乗算」演算によって取得する必要があります。明らかに、したがって、全体の乗算は、を取得でき、次に1を乗算してを取得できます。したがって、f2とf1は既知であり、行列を使用してそれらをすばやく指数化します。
(コード自身の水)
例2:TRのシーケンス
TRは数学が大好きで、奇妙な数学の問題を研究するためにドラフトペーパーをよく取り出します。最近、彼は突然数字のシーケンスに興味を持ち、フィボナッチに似た数字のシーケンス、つまりTn = 1 * f1 + 2 *を見つけました。 f2 + 3 * f3 +……+ n * fn(fnはフィボナッチのn番目の項の値)、TRはTn%m(1≤n、m≤2^ 31)の値を見つけるのを手伝ってくれるように頼みます-1)。
とすると、
したがって(キー再帰)、左側の行列はに設定され、次に
だから。
わかりますか?最後に、行列の高速べき乗を使用します。
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
long long n,mod;
struct matrix{
int n,m;
long long c[105][105];
matrix(){memset(c,0,sizeof(c));}
matrix operator*(const matrix&a){
matrix r;r.n=n,r.m=a.m;
for(int i=1;i<=r.n;i++)
for(int j=1;j<=r.m;j++)
for(int k=1;k<=m;k++)r.c[i][j]=(r.c[i][j]+c[i][k]*a.c[k][j])%mod;
return r;
}
}a,b;
matrix mpow(matrix a,long long b){
matrix res;res.n=res.m=a.n;
for(int i=1;i<=res.n;i++)res.c[i][i]=1;
for(;b;b>>=1){
if(b&1)res=res*a;
a=a*a;
}
return res;
}
int main()
{
scanf("%lld%lld",&n,&mod);
if(n==1){
printf("%lld",1%mod);
putchar('\n');
return 0;
}
a.n=1,a.m=5;
a.c[1][1]=3,a.c[1][2]=a.c[1][4]=a.c[1][5]=1,a.c[1][3]=2;
b.n=b.m=5;
b.c[1][1]=1;
b.c[2][1]=b.c[2][3]=b.c[3][1]=b.c[3][2]=b.c[3][3]=1;
b.c[4][1]=b.c[4][3]=2,b.c[4][5]=1;
b.c[5][1]=b.c[5][3]=b.c[5][4]=b.c[5][5]=1;
a=a*mpow(b,n-2);
printf("%lld",a.c[1][1]%mod);
putchar('\n');
return 0;
}
総括する
(たとえば、質問1と2、1つの解決策を理解するだけです)行列加速法は、最初に主要な漸化式をリストし(右側は1回のm項でなければなりません)、次に1 * mとm *をリストします。式に従ってm行列、そしてそれはとても簡単です。
(update2021.3.13)拡張:一般化モーメント乗算
上記の方法は、プラス記号とマイナス記号で構成される再帰式のみを解くようです。実際、最大記号と最小記号を使用して接続された式は、行列で最適化することもできます。このとき、一般化された行列乗算が使用されます。
、(これはそのうちの1つにすぎません)
一般化された行列の乗算は結合法則を満たしますが、分布率を満たさないことが証明できますが、行列の加速には結合法則のみが必要なため、これで十分です。
プラス記号とマイナス記号の制約を取り除くため、行列アクセラレーションは、一見最適化されていない多くの問題を最適化できます。動的DPは、一般化モーメント乗算に基づいて確立されたアルゴリズムです。