网站流量统计软件,价目表海报app制作,一个简单的网站怎么做的,广州市天河区住房和建设局网站应队友要求#xff0c;开始学线性代数#xff0c;具体路线是矩阵 → \rightarrow →高斯消元 → \rightarrow →线性基。为多项式做个准备
P3390 【模板】矩阵快速幂
题面 板子#xff0c;用结构体写的#xff0c;感觉有点丑#xff0c;一会儿看看题解有没有写得好看的 …应队友要求开始学线性代数具体路线是矩阵 → \rightarrow →高斯消元 → \rightarrow →线性基。为多项式做个准备
P3390 【模板】矩阵快速幂
题面 板子用结构体写的感觉有点丑一会儿看看题解有没有写得好看的
#includebits/stdc.h
using namespace std;
typedef long long ll;
const int N 110;
const ll mod1e97;
struct node{ll a[N][N];int len;}sqr;
void sqr0(node x){memset(x.a,0,sizeof x.a);x.lensqr.len;
}
void sqr1(node x){memset(x.a,0,sizeof x.a);x.lensqr.len;for(int i1;ix.len;i)x.a[i][i]1;
}
node operator*(node x, node b){node c;sqr0(c);for(int i1;ix.len;i){for(int j1;jx.len;j){for(int k1;kx.len;k)(c.a[i][j]x.a[i][k]*b.a[k][j]%mod)%mod;}}return c;
}void qpow(node x, ll y){node re;sqr1(re);while(y){if(y1)rere*x;xx*x;y1;}xre;
}
ll k;
int main(){scanf(%d%lld,sqr.len,k);for(int i1;isqr.len;i){for(int j1;jsqr.len;j)scanf(%lld,sqr.a[i][j]);}qpow(sqr,k);for(int i1;isqr.len;i){for(int j1;jsqr.len;j)printf(%lld ,sqr.a[i][j]);puts();}
}
P1939 【模板】矩阵加速数列
题面 搞个方阵 A 3 [ a 3 a 2 a 1 0 0 0 0 0 0 ] , X [ 1 1 0 0 0 1 1 0 0 ] , A_3\left [ \begin{matrix} a_3 a_2 a_1 \\ 0 0 0 \\ 0 0 0 \\ \end{matrix} \right] ,X\left [ \begin{matrix} 1 1 0 \\ 0 0 1 \\ 1 0 0 \\ \end{matrix} \right], A3 a300a200a100 ,X 101100010 , 则 A 3 X [ a 4 a 3 a 2 0 0 0 0 0 0 ] A 4 , A_3X\left [ \begin{matrix} a_4 a_3 a_2 \\ 0 0 0 \\ 0 0 0 \\ \end{matrix} \right]A_4, A3X a400a300a200 A4, 因此对 X X X进行矩阵快速幂即可。
#includebits/stdc.h
using namespace std;
typedef long long ll;
const int N 5;
const ll mod1e97;
struct node{ll a[N][N];}sqr,A;
void sqr0(node x){memset(x.a,0,sizeof x.a);
}
void sqr1(node x){memset(x.a,0,sizeof x.a);for(int i1;i3;i)x.a[i][i]1;
}
node operator*(node x, node b){node c;sqr0(c);for(int i1;i3;i){for(int j1;j3;j){for(int k1;k3;k)(c.a[i][j]x.a[i][k]*b.a[k][j]%mod)%mod;}}return c;
}void qpow(node x, ll y){node re;sqr1(re);while(y){if(y1)rere*x;xx*x;y1;}xre;
}
ll n,T;
int main(){cinT;while(T--){cinn;if(n3){puts(1);continue;}sqr0(sqr);sqr.a[1][1]sqr.a[1][2]sqr.a[2][3]sqr.a[3][1]1;sqr0(A);A.a[1][1]A.a[1][2]A.a[1][3]1;qpow(sqr,n-3);AA*sqr;coutA.a[1][1]endl;}}
P4783 【模板】矩阵求逆
题面 把一个矩阵通过行变换变为单位矩阵所需要的行变换操作操作给一个单位矩阵就可以得到其逆矩阵。故应用高斯消元即可。
#includebits/stdc.h
#define N 1000
using namespace std;
const int mod1e97;
inline void read(int x){int s0,w1;char chgetchar();while(ch0||ch9){if(ch-)w-1;chgetchar();}while(ch0ch9){s(s3)(s1)(ch15);chgetchar();}xs*w;
}
int n,a[N][N];
int qpow(int x, int y){int re1;while(y){if(y1)re1LL*re*x%mod;x1LL*x*x%mod,y1;}return re;
}
int main(){read(n);for(int i1;in;i){for(int j1;jn;j)read(a[i][j]);a[i][ni]1;}for(int i1;in;i){int nowi;for(int ji;jn;j)if(a[now][i]a[j][i])nowj;if(a[now][i]0){puts(No Solution);return 0;}if(now!i)swap(a[now],a[i]);for(int ji1;jn1;j)a[i][j]1LL*a[i][j]*qpow(a[i][i],mod-2)%mod;a[i][i]1;for(int j1;jn;j){if(ji)continue;int div1LL*a[j][i]*qpow(a[i][i],mod-2)%mod;for(int ki;kn1;k)a[j][k](a[j][k]-1LL*a[i][k]*div%modmod)%mod;}}for(int i1;in;i,puts())for(int j1;jn;j)printf(%d ,a[i][nj]);}P1962 斐波那契数列
题面 构造矩阵 A 2 [ f 2 f 1 0 0 ] , X [ 1 1 1 0 ] , A_2\left [ \begin{matrix} f_2 f_1 \\ 0 0 \\ \end{matrix} \right] ,X\left [ \begin{matrix} 1 1 \\ 1 0 \\ \end{matrix} \right], A2[f20f10],X[1110], 则 A 2 X [ f 3 f 2 0 0 ] A 3 , A_2X\left [ \begin{matrix} f_3 f_2 \\ 0 0 \\ \end{matrix} \right]A_3, A2X[f30f20]A3,
#includecstdio
typedef long long ll;
const ll modll(1e97);
struct node
{ll sqr[5][5];
}a;
node operator*(node a, node b)
{node c;c.sqr[1][1](a.sqr[1][1]*b.sqr[1][1]%moda.sqr[1][2]*b.sqr[2][1]%mod)%mod;c.sqr[1][2](a.sqr[1][1]*b.sqr[1][2]%moda.sqr[1][2]*b.sqr[2][2]%mod)%mod;c.sqr[2][1](a.sqr[2][1]*b.sqr[1][1]%moda.sqr[2][2]*b.sqr[2][1]%mod)%mod;c.sqr[2][2](a.sqr[2][1]*b.sqr[1][2]%moda.sqr[2][2]*b.sqr[2][2]%mod)%mod;return c;
}
ll n;
void quickpow(node x, ll y)
{node rec;rec.sqr[1][1]rec.sqr[2][2]1,rec.sqr[1][2]rec.sqr[2][1]0;while(y){if(y1)recrec*x;xx*x,y1;}xrec;
}
int main()
{scanf(%lld,n);if(n0)return puts(0);a.sqr[1][1]a.sqr[1][2]a.sqr[2][1]1,a.sqr[2][2]0;quickpow(a,n-1);printf(%lld\n,a.sqr[1][1]);
}P1349 广义斐波那契数列
题面 构造矩阵 A 2 [ f 2 f 1 0 0 ] , X [ P 1 Q 0 ] , A_2\left [ \begin{matrix} f_2 f_1 \\ 0 0 \\ \end{matrix} \right] ,X\left [ \begin{matrix} P 1 \\ Q 0 \\ \end{matrix} \right], A2[f20f10],X[PQ10], 则 A 2 X [ f 3 f 2 0 0 ] A 3 , A_2X\left [ \begin{matrix} f_3 f_2 \\ 0 0 \\ \end{matrix} \right]A_3, A2X[f30f20]A3,
#includebits/stdc.h
using namespace std;
typedef long long ll;
ll mod;
struct node
{ll sqr[5][5];node(){memset(sqr,0,sizeof sqr);}
}a,b;
node operator*(node a, node b)
{node c;c.sqr[1][1](a.sqr[1][1]*b.sqr[1][1]%moda.sqr[1][2]*b.sqr[2][1]%mod)%mod;c.sqr[1][2](a.sqr[1][1]*b.sqr[1][2]%moda.sqr[1][2]*b.sqr[2][2]%mod)%mod;c.sqr[2][1](a.sqr[2][1]*b.sqr[1][1]%moda.sqr[2][2]*b.sqr[2][1]%mod)%mod;c.sqr[2][2](a.sqr[2][1]*b.sqr[1][2]%moda.sqr[2][2]*b.sqr[2][2]%mod)%mod;return c;
}
ll n;
void quickpow(node x, ll y)
{node rec;rec.sqr[1][1]rec.sqr[2][2]1,rec.sqr[1][2]rec.sqr[2][1]0;while(y){if(y1)recrec*x;xx*x,y1;}xrec;
}
int main()
{scanf(%lld%lld%lld%lld%lld%lld,a.sqr[1][1],a.sqr[2][1],b.sqr[1][2],b.sqr[1][1],n,mod);if(n2)return printf(%lld\n,b.sqr[1][3-n]);a.sqr[1][2]1,a.sqr[2][2]0;quickpow(a,n-2);bb*a;printf(%lld\n,b.sqr[1][1]);
}P4000 斐波那契数列
题面 不是这什么题都往题单里放啊这是我 18 18 18年外出集训堵了个论文费了两三天时间才切了的人生中第一道黑题现在变成紫题了。 有一个性质是 f n m o d p f_n\mod p fnmodp有循环节且循环节长度不会超过 6 p 6p 6p这还有个名叫皮萨诺定理。所以我们考虑求出循环节的长度然后用矩阵乘法求出结果。 引理对于 f n m o d p f_n\mod p fnmodp的循环节 g ( p ) g(p) g(p)有如下性质 1. p p i α i 1.pp_i^{\alpha_i} 1.ppiαi即 p p p为质数的幂时 g ( p ) g ( p i ) × p i α i − 1 g(p)g(p_i)\times p_i^{\alpha_i-1} g(p)g(pi)×piαi−1 2. p ∏ p i α i 2.p\prod p_i^{\alpha_i} 2.p∏piαi即 p p p为合数时 g ( p ) l c m ( g ( p i α i ) g(p)lcm(g(p_i^{\alpha_i}) g(p)lcm(g(piαi) 对于 g ( p ) g(p) g(p)这么算如果 5 5 5是模 p p p的二次剩余那么循环节为 p − 1 p-1 p−1的因子否则为 2 p 2 2p2 2p2的因子。 因为 p p p不是特别大直接取 p − 1 p-1 p−1和 2 p 2 2p2 2p2即可。 对于 p ≤ 5 p\le 5 p≤5就暴力算即可 g ( 2 ) 3 , g ( 3 ) 5 , g ( 5 ) 20 g(2)3,g(3)5,g(5)20 g(2)3,g(3)5,g(5)20
// luogu-judger-enable-o2
#includecstdio
#includecmath
#includecstring
using namespace std;
#define rg register
typedef long long ll;
char str[30000000];
ll n,p,mod,len,fac[100000],power[100000],faccnt,s;
struct node
{ll sqr[5][5];
}b;
node operator *(node a, node b)
{node xx;xx.sqr[1][1](a.sqr[1][1]*b.sqr[1][1]%pa.sqr[1][2]*b.sqr[2][1]%p)%p;xx.sqr[1][2](a.sqr[1][1]*b.sqr[1][2]%pa.sqr[1][2]*b.sqr[2][2]%p)%p;xx.sqr[2][2](a.sqr[2][1]*b.sqr[1][2]%pa.sqr[2][2]*b.sqr[2][2]%p)%p;xx.sqr[2][1](a.sqr[2][1]*b.sqr[1][1]%pa.sqr[2][2]*b.sqr[2][1]%p)%p;return xx;
}
node quickpow(node x, ll y)
{node rec;rec.sqr[1][1]rec.sqr[2][2]1;rec.sqr[1][2]rec.sqr[2][1]0;while(y){if(y%21)recrec*x;xx*x;y/2;}return rec;
}
ll gcd(ll a, ll b)
{if(b0)return a;else return gcd(b,a%b);
}
ll lcm(ll a, ll b)
{return a*b/gcd(a,b);
}
ll get(ll k)
{ll nowk;for(rg ll i2;i*inow;i){if(now%i0){faccnt;fac[faccnt]i;power[faccnt]1;while(now%i0){now/i;power[faccnt]*i;}}}for(rg ll i1;ifaccnt;i)power[i]/fac[i];if(now!1){fac[faccnt]now;power[faccnt]1;}for(rg ll i1;ifaccnt;i){if(fac[i]2)power[i]*3;else if(fac[i]3)power[i]*5;else if(fac[i]5)power[i]*20;else if(fac[i]%51||fac[i]%54)power[i]*fac[i]-1;else power[i]*(fac[i]1)1;}ll anspower[1];for(rg ll i1;ifaccnt;i)anslcm(ans,power[i]);return ans;
}
int main()
{scanf(%s%lld,str,p);if(p1){printf(0\n);return 0;}modget(p);lenstrlen(str);for(rg ll i0;ilen;i)n((n3)(n1)(str[i]15))%mod;if(n0){printf(0\n);return 0;}if(n1||n2){printf(1\n);return 0;}b.sqr[2][2]0;b.sqr[1][1]b.sqr[1][2]b.sqr[2][1]1;bquickpow(b,n-1);printf(%lld\n,b.sqr[1][1]);
}P3758 [TJOI2017] 可乐
题面
#includebits/stdc.h
#define N 50
using namespace std;
const int mod2017;
int t,n,m;
struct node{int a[N][N];node(){memset(a,0,sizeof a);}
}sqr;
node operator*(node x, node b){node c;for(int i0;in;i){for(int j0;jn;j){for(int k0;kn;k)(c.a[i][j]x.a[i][k]*b.a[k][j]%mod)%mod;}}return c;
}
void qpow(node x, int y){node re;for(int i0;in;i)re.a[i][i]1;while(y){if(y1)rere*x;xx*x;y1;}xre;
}
int main(){cinnm;for(int i1,u,v;im;i)cinuv,sqr.a[u][v]sqr.a[v][u]1;cint;for(int i1;in;i)sqr.a[i][0]1;for(int i0;in;i)sqr.a[i][i]1;qpow(sqr,t);int ans0;for(int i0;in;i)(anssqr.a[1][i])%mod;coutansendl;}P5343 【XR-1】分块
题面 方程很简单 d p [ i ] ∑ j ∈ b l o c k , j ≤ i d p [ i − j ] dp[i]\sum_{j\in block,j\le i}dp[i-j] dp[i]∑j∈block,j≤idp[i−j]现在考虑如何矩阵优化。 由于块的大小不会超过 100 100 100所以我们开一个 100 × 100 100\times 100 100×100的矩阵首先预处理出 d p [ 1 ] − d p [ 100 ] dp[1]-dp[100] dp[1]−dp[100]将其填入 A A A矩阵第一行中再考虑所有 j ∈ b l o c k j\in block j∈block设 X [ j ] [ 1 ] 1 X[j][1]1 X[j][1]1对于后99列设 X [ j − 1 ] [ j ] 1 X[j-1][j]1 X[j−1][j]1则这样可以转移 A A A矩阵应用矩阵快速幂即可。
#includebits/stdc.h
#define N 120
using namespace std;
const int mod1e97;
const int len100;
inline void read(int x){int s0,w1;char chgetchar();while(ch0||ch9){if(ch-)w-1;chgetchar();}while(ch0ch9){s(s3)(s1)(ch15);chgetchar();}xs*w;
}
long long n;
int p,q,cnt,a[N],f[N],vis[N];
setint s;
struct node{int m[N][N];node(){memset(m,0,sizeof m);}
}sqr,A;
node operator*(node a, node b){node c;for(int i0;ilen;i)for(int j0;jlen;j)for(int k0;klen;k)c.m[i][j]int((1LL*c.m[i][j]1LL*a.m[i][k]*b.m[k][j]%mod)%mod);return c;
}
void qpow(node x, long long y){node re;for(int i0;ilen;i)re.m[i][i]1;while(y){if(y1)rere*x;xx*x;y1;}xre;
}
int main(){cinn;read(p);for(int i1,x;ip;i){read(x);if(s.find(x)s.end())s.insert(x);}read(q);for(int i1,x;iq;i){read(x);if(s.find(x)!s.end()!vis[x])a[cnt]x,vis[x]1;}f[0]1;for(int i1;ilen;i)for(int j1;jcnt;j)if(a[j]i)f[i](1LL*f[i]f[i-a[j]])%mod;if(n100){printf(%d\n,f[n]);return 0;}for(int i0;ilen;i)sqr.m[0][len-i]f[i];for(int i1;icnt;i)A.m[a[i]-1][0]1;for(int i1;ilen;i)A.m[i-1][i]1;qpow(A,n-100);sqrsqr*A;printf(%d\n,sqr.m[0][0]);return 0;
}
文章转载自： http://www.morning.lgnbr.cn.gov.cn.lgnbr.cn http://www.morning.qlkjh.cn.gov.cn.qlkjh.cn http://www.morning.gnmhy.cn.gov.cn.gnmhy.cn http://www.morning.wjyyg.cn.gov.cn.wjyyg.cn http://www.morning.fgxnb.cn.gov.cn.fgxnb.cn http://www.morning.jcfg.cn.gov.cn.jcfg.cn http://www.morning.wcrcy.cn.gov.cn.wcrcy.cn http://www.morning.hrpbq.cn.gov.cn.hrpbq.cn http://www.morning.hdrrk.cn.gov.cn.hdrrk.cn http://www.morning.gwjqq.cn.gov.cn.gwjqq.cn http://www.morning.coffeedelsol.com.gov.cn.coffeedelsol.com http://www.morning.rykw.cn.gov.cn.rykw.cn http://www.morning.ljfjm.cn.gov.cn.ljfjm.cn http://www.morning.yhljc.cn.gov.cn.yhljc.cn http://www.morning.wgbmj.cn.gov.cn.wgbmj.cn http://www.morning.trrpb.cn.gov.cn.trrpb.cn http://www.morning.srmpc.cn.gov.cn.srmpc.cn http://www.morning.kstgt.cn.gov.cn.kstgt.cn http://www.morning.jxtbr.cn.gov.cn.jxtbr.cn http://www.morning.xwgbr.cn.gov.cn.xwgbr.cn http://www.morning.hwhnx.cn.gov.cn.hwhnx.cn http://www.morning.smrkf.cn.gov.cn.smrkf.cn http://www.morning.kzpy.cn.gov.cn.kzpy.cn http://www.morning.hjrjr.cn.gov.cn.hjrjr.cn http://www.morning.kyflr.cn.gov.cn.kyflr.cn http://www.morning.xzlp.cn.gov.cn.xzlp.cn http://www.morning.zkzjm.cn.gov.cn.zkzjm.cn http://www.morning.rkxdp.cn.gov.cn.rkxdp.cn http://www.morning.kndst.cn.gov.cn.kndst.cn http://www.morning.pqcbx.cn.gov.cn.pqcbx.cn http://www.morning.tntqr.cn.gov.cn.tntqr.cn http://www.morning.sjsks.cn.gov.cn.sjsks.cn http://www.morning.tbrnl.cn.gov.cn.tbrnl.cn http://www.morning.ymjrg.cn.gov.cn.ymjrg.cn http://www.morning.pkpqh.cn.gov.cn.pkpqh.cn http://www.morning.ywxln.cn.gov.cn.ywxln.cn http://www.morning.jkpnm.cn.gov.cn.jkpnm.cn http://www.morning.bpmdn.cn.gov.cn.bpmdn.cn http://www.morning.skbkq.cn.gov.cn.skbkq.cn http://www.morning.crrmg.cn.gov.cn.crrmg.cn http://www.morning.kjmcq.cn.gov.cn.kjmcq.cn http://www.morning.xrpwk.cn.gov.cn.xrpwk.cn http://www.morning.psqs.cn.gov.cn.psqs.cn http://www.morning.plhyc.cn.gov.cn.plhyc.cn http://www.morning.mtbth.cn.gov.cn.mtbth.cn http://www.morning.bsjxh.cn.gov.cn.bsjxh.cn http://www.morning.mqdr.cn.gov.cn.mqdr.cn http://www.morning.rbbzn.cn.gov.cn.rbbzn.cn http://www.morning.ycnqk.cn.gov.cn.ycnqk.cn http://www.morning.qnxzx.cn.gov.cn.qnxzx.cn http://www.morning.rwzkp.cn.gov.cn.rwzkp.cn http://www.morning.rckmz.cn.gov.cn.rckmz.cn http://www.morning.jzkqg.cn.gov.cn.jzkqg.cn http://www.morning.tongweishi.cn.gov.cn.tongweishi.cn http://www.morning.tpnxj.cn.gov.cn.tpnxj.cn http://www.morning.trzmb.cn.gov.cn.trzmb.cn http://www.morning.drnfc.cn.gov.cn.drnfc.cn http://www.morning.bynf.cn.gov.cn.bynf.cn http://www.morning.bdsyu.cn.gov.cn.bdsyu.cn http://www.morning.lnsnyc.com.gov.cn.lnsnyc.com http://www.morning.frpfk.cn.gov.cn.frpfk.cn http://www.morning.thmlt.cn.gov.cn.thmlt.cn http://www.morning.xgzwj.cn.gov.cn.xgzwj.cn http://www.morning.pthmn.cn.gov.cn.pthmn.cn http://www.morning.nqbs.cn.gov.cn.nqbs.cn http://www.morning.kjyfq.cn.gov.cn.kjyfq.cn http://www.morning.rfbt.cn.gov.cn.rfbt.cn http://www.morning.nfpct.cn.gov.cn.nfpct.cn http://www.morning.qrqg.cn.gov.cn.qrqg.cn http://www.morning.kpxky.cn.gov.cn.kpxky.cn http://www.morning.cqyhdy.cn.gov.cn.cqyhdy.cn http://www.morning.lcmhq.cn.gov.cn.lcmhq.cn http://www.morning.yrskc.cn.gov.cn.yrskc.cn http://www.morning.gpxbc.cn.gov.cn.gpxbc.cn http://www.morning.zfhwm.cn.gov.cn.zfhwm.cn http://www.morning.dmlsk.cn.gov.cn.dmlsk.cn http://www.morning.gcrlb.cn.gov.cn.gcrlb.cn http://www.morning.brscd.cn.gov.cn.brscd.cn http://www.morning.rfhm.cn.gov.cn.rfhm.cn http://www.morning.hhqjf.cn.gov.cn.hhqjf.cn
