Codeforces 575A

题面传送门

题意：
给出 (s_0,s_1,s_2,dots,s_{n-1})，对于 (igeq n)，有 (m) 个 (s_i) 满足 (s_i eq s_{imod n})，这 (m) 个 (s_i) 特殊给定，其它的 (s_i=s_{imod n})。
另有序列 (f_i) 满足 (f_0=0,f_1=1)，(f_i=f_{i-1}s_{i-1}+f_{i-2}s_{i-2})。
求 (f_xmod P)。
(1leq n,mleq 5 imes 10^4,0leq xleq 10^{18})

Fibo纳粹
粉兔：矩阵快速幂 sb 题。
我：矩阵快速幂神仙题。
这好像是上赛季的一个 jxd 作业。
首先写出矩阵递推式：

([f_{i-2},f_{i-1}] imesegin{bmatrix}0&s_{i-2}\1&s_{i-1}end{bmatrix}=[f_{i-1},f_i])

于是 ([f_{x-1},f_x]=[0,1] imesprodlimits_{i=1}^{x-1}egin{bmatrix}0&s_{i-1}\1&s_{i}end{bmatrix})
首先考虑 (m=0) 的情况。后面 (prod) 的矩阵成周期性分部。所以可以求出每个周期内所有矩阵的乘积，跑个矩阵快速幂，最后的零头暴力算一下就可以了。
接下来考虑 (m eq 0) 的情况。记 (P=prodlimits_{i=0}^{n-1}egin{bmatrix}0&s_i\1&s_{(i+1)mod n}end{bmatrix})，不难发现，对于大多数周期，其中所有矩阵的乘积都为 (P)，最多只有 (2m) 个周期的乘积 ( eq P)。
故可以对这些特殊的周期特殊计算，其它周期跑矩阵快速幂。
考虑一个特殊点 (x_i) 会造成什么影响，会使 (x_i) 和 (x_i-1) 的转移矩阵发生变化。
所以我们对于所有 (x_i) 和 (x_i-1) 所在的周期单独计算它的贡献。
我们建立 (n) 个矩阵 (A_0,A_1,A_2,dots,A_{n-1})，初始 (A_i=egin{bmatrix}0&s_i\1&s_{(i+1)mod n}end{bmatrix})
对于某个周期中特殊处理的点 (p)，我们修改 (pmod n) 转移矩阵的值，然后求一遍所有 (A_i) 的乘积。处理完这个周期后再改回去。
于是此题变为：修改某个矩阵的值，求所有矩阵的乘积。
注意到矩阵是不支持除法操作的，所以这玩意儿只能用线段树维护。
时间复杂度 (2^3nlog n)。
细节恶心得要死。。。。

#include <bits/stdc++.h>
using namespace std;
#define fi first
#define se second
#define fz(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
#define ffe(it,v) for(__typeof(v.begin()) it=v.begin();it!=v.end();it++)
#define fill0(a) memset(a,0,sizeof(a))
#define fill1(a) memset(a,-1,sizeof(a))
#define fillbig(a) memset(a,63,sizeof(a))
#define pb push_back
#define ppb pop_back
#define mp make_pair
template<typename T1,typename T2> void chkmin(T1 &x,T2 y){if(x>y) x=y;}
template<typename T1,typename T2> void chkmax(T1 &x,T2 y){if(x<y) x=y;}
typedef pair<int,int> pii;
typedef long long ll;
template<typename T> void read(T &x){
	char c=getchar();T neg=1;
	while(!isdigit(c)){if(c=='-') neg=-1;c=getchar();}
	while(isdigit(c)) x=x*10+c-'0',c=getchar();
	x*=neg;
}
const int MAXN=5e4;
const int SIZE=2;
ll x;int mod,n,m,k,b[MAXN+5];
pair<ll,int> p[MAXN+5];
ll tk[MAXN*2+5];
struct mat{
	int a[SIZE][SIZE];
	mat(){memset(a,0,sizeof(a));}
	friend mat operator *(mat x,mat y){
		mat ret;
		for(int i=0;i<2;i++) for(int j=0;j<2;j++) for(int k=0;k<2;k++)
			ret.a[i][j]=(ret.a[i][j]+1ll*x.a[i][k]*y.a[k][j]%mod)%mod;
		return ret;
	}
};
struct node{int l,r;mat val;} s[MAXN*4+5];
void build(int k,int l,int r){
	s[k].l=l;s[k].r=r;if(l==r) return;
	int mid=(l+r)>>1;build(k<<1,l,mid);build(k<<1|1,mid+1,r);
}
void modify(int k,int p,mat x){
	if(s[k].l==s[k].r){s[k].val=x;return;}
	int mid=(s[k].l+s[k].r)>>1;
	if(p<=mid) modify(k<<1,p,x);else modify(k<<1|1,p,x);
	s[k].val=s[k<<1].val*s[k<<1|1].val;
}
mat query(int k,int l,int r){
	if(l<=s[k].l&&s[k].r<=r) return s[k].val;
	int mid=(s[k].l+s[k].r)>>1;
	if(r<=mid) return query(k<<1,l,r);
	else if(l>mid) return query(k<<1|1,l,r);
	else return query(k<<1,l,mid)*query(k<<1|1,mid+1,r);
}
mat qpow(mat x,ll e){
	mat ret;ret.a[0][0]=ret.a[1][1]=1;
	while(e){if(e&1) ret=ret*x;x=x*x;e>>=1;}
	return ret;
}
map<ll,int> mmp;
signed main(){
//	freopen("in.txt","r",stdin);
	scanf("%lld%d%d",&x,&mod,&n);x--;
	if(x==-1) return printf("0
"),0;
	for(int i=0;i<n;i++) scanf("%d",&b[i]);scanf("%d",&m);
	for(int i=1;i<=m;i++) scanf("%lld%d",&p[i].fi,&p[i].se);
	for(int i=1;i<=m;i++) mmp[p[i].fi]=p[i].se;
	for(int i=1;i<=m;i++){
		tk[++k]=p[i].fi;
		if(p[i].fi!=n) tk[++k]=p[i].fi-1;
	}
	sort(tk+1,tk+k+1);k=unique(tk+1,tk+k+1)-tk-1;build(1,0,n-1);
	for(int i=0;i<n;i++){
		mat t;t.a[0][0]=0;t.a[1][0]=1;t.a[0][1]=b[i];
		if(!mmp[i+1]) t.a[1][1]=b[(i+1)%n];else t.a[1][1]=mmp[i+1];
		modify(1,i,t);
	}
	mat ret;
	if(x<n) ret=query(1,0,x);
	else{
		ret=query(1,0,n-1);
		if(mmp[n]){
			mat t;t.a[0][0]=0;t.a[1][0]=1;t.a[0][1]=b[n-1];t.a[1][1]=b[0];
			modify(1,n-1,t);
		}
		int cur=1,pre=0;ll pn=0;bool flg=0;
		while(cur<=k){
//			puts("-1");
//			printf("%d %d %d
",cur,pn,pre);
			if(x/n<tk[pre+1]/n){
				ret=ret*qpow(query(1,0,n-1),(x/n)-pn-1);
				ret=ret*query(1,0,x%n);flg=1;break;
			}
			ret=ret*qpow(query(1,0,n-1),(tk[pre+1]/n)-pn-1);
			while(cur<=k&&tk[cur]/n==tk[pre+1]/n){
				ll id=tk[cur];int v1,v2;
				if(!mmp[id]) v1=b[id%n];else v1=mmp[id];
				if(!mmp[id+1]) v2=b[(id+1)%n];else v2=mmp[id+1];
//				printf("%lld %d %d %d
",id,id%n,v1,v2);
				mat t;t.a[0][0]=0;t.a[1][0]=1;t.a[0][1]=v1;t.a[1][1]=v2;
				modify(1,id%n,t);cur++;
			}
			if((x/n)==(tk[pre+1]/n)){
				ret=ret*query(1,0,x%n);flg=1;break;
			} else ret=ret*query(1,0,n-1);
			for(int i=pre+1;i<cur;i++){
				mat t;t.a[0][0]=0;t.a[1][0]=1;
				t.a[0][1]=b[tk[i]%n];t.a[1][1]=b[(tk[i]+1)%n];
				modify(1,tk[i]%n,t);
			}
			pn=tk[pre+1]/n;pre=cur-1;
		}
		if(!flg){
//			puts("-1");
			ret=ret*qpow(query(1,0,n-1),(x/n)-pn-1);
			ret=ret*query(1,0,x%n);
		}
	}
	mat f;f.a[0][1]=1;f=f*ret;
	printf("%d
",f.a[0][0]);
	return 0;
}
/*
3 998244353
3
1 2 1
3
3 3
4 6
10 8

20 998244353
3
100 370 250
5
7 230
23 500
5 480
15 530
19 570

100 998244353
3
1 2 1
2
3 2
5 1

1000000 998244353
3
100 370 250
5
7 230
23 500
5 480
15 530
19 570
*/