Bubble Cup 8 finals A. Fibonotci (575A)

题意：

定义类循环序列为 长度无限，且除了有限个元素外，均满足s[i] ≡ s[i mod N] (i≥N)。

现在有数列F，定义为 F[i] = s[i-2]*F[i-1] + s[i-1]*F[i-1]，特别的，F[0] = 0, F[1] = 1。

给定正整数K,P,N代表要求输出的答案为F[k] mod P，类循环序列s的长度为N。

接下来给出s[0]..s[n-1]。

然后是一个正整数M，代表不满足循环条件的元素个数。

接下来M行每行两个正整数j,v表示s[j] = v，保证所有j不同。

1<=N,M<=10^5,

1<=P,s[i],v<=10^9,

1<=K,j<=10^18

 

题解：

会矩阵乘法快速转移的一看就知道大概的做法…所以思考难度不算高。

但是实现起来比较麻烦。

可以明显地看出，没有被特殊位置影响到递推式的部分可以快速转移。

所以我们排序并求出所有被特殊位置影响到递推式的位置（即特殊位置的后两个位置）。

然后按照排序相邻的两个受影响位置的位置关系分类来做：

1、两个受影响位置分离，中间包含若干个循环

那么我们可以在预处理时求出一个循环内所有转移矩阵的乘积，

然后找出当前位置所处的循环节末尾，下一个位置所处的循环节开头，

把这中间的部分用快速幂处理出来。

而当前位置到循环节末的部分，可以用反向的前缀和或线段树处理。

循环节头到下一个位置的部分，可以用前缀和或线段树处理。

（这里用前缀和是O(1)，但是因为还有快速幂，总的复杂度没有降低）

2、两个受影响位置分离，在同一循环节之中

这种情况可以直接用线段树求处理出区间的乘积。

3、两个受影响位置相邻

直接暴力计算。

这种做法不用考虑连续的受影响位置跨越循环节的情况，写起来比较方便，不那么容易写挂。

（然而我因为诸如没开long long这类的低级错误挂了很多次）

时间复杂度O( (2^3) * M log N )

 

  1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
 
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
using namespace std;
const int N = 100010, NS = 131072+10;
typedef long long lint;
typedef pair<lint,int> pli;
inline int read()
{
	int s = 0; char c; while((c=getchar())<'0'||c>'9');
	do{s=s*10+c-'0';}while((c=getchar())>='0'&&c<='9');
	return s;
}
inline lint readll()
{
	lint s = 0;char c; while((c=getchar())<'0'||c>'9');
	do{s=s*10+c-'0';}while((c=getchar())>='0'&&c<='9');
	return s;
}
int n,m,p,s[N],br[N+N],cr[N+N],cur,S,tot;
lint k,ar[N+N];
pli q[N];
struct mrx
{
	int n[2][2];
	void clr(){ memset(n,0,sizeof n); }
	friend mrx operator * (const mrx &a,const mrx &b)
	{
		mrx c; c.clr();
		for(int i=0;i<2;i++) for(int k=0;k<2;k++) if(a.n[i][k]) for(int j=0;j<2;j++) 
			c.n[i][j] = ((lint)c.n[i][j]+(lint)a.n[i][k]*b.n[k][j])%p;
		return c;
	}
}ym,tm[N],stm,I,tmpm,tr[NS],psm[N],ssm[N];
inline lint mo(lint a){ return (a+n-1)%n+1; }
mrx powmod(mrx a,lint b)
{
	mrx ans = I;
	for(;b;b>>=1)
	{
		if(b&1) ans = ans*a;
		a = a*a;
	}
	return ans;
}
mrx query(int l,int r)
{
	if(l>r) return I;
	mrx ls = I, rs = I;
	for(l=l+S-1,r=r+S+1;l^r^1;l>>=1,r>>=1)
	{
		if(~l&1) ls = ls*tr[l^1];
		if( r&1) rs = tr[r^1]*rs;
	}
	return ls*rs;
}
int calc()
{
	int i;
	for(i=1;i<=m;i++)
	{
		if(q[i].first>=k) break;
		ar[++tot] = q[i].first+1;
		br[tot] = (q[i].first==q[i-1].first+1)?q[i-1].second:s[(q[i].first+n-1)%n+1];
		cr[tot] = q[i].second;
		if(q[i].first+1!=q[i+1].first&&q[i].first+2<=k)
			ar[++tot] = q[i].first+2, br[tot] = cr[tot-1], cr[tot] = s[(q[i].first+1)%n+1];
	}
	tmpm = tm[1];
	cur = 1; while(cur<=tot&&ar[cur]<2) cur++;
	ar[--cur] = 1;
	ar[++tot] = k+1;
	lint ql,qr;
	for(;cur<tot;cur++)
	{
		if(cur+1!=tot&&ar[cur+1]==ar[cur]+1)
		{
			tmpm.n[0][1] = br[cur+1], tmpm.n[1][1] = cr[cur+1];
			ym = ym*tmpm;
			continue;
		}
		ql = (lint)((ar[cur]-1)/n)*n+n;
		qr = (lint)((ar[cur+1]-1)/n)*n+1;
		if(ql<qr)
			ym = ((ym * ssm[(ar[cur]-1)%n+2]) * powmod(psm[n],(qr-ql-1)/n)) * psm[ar[cur+1]-qr];
		else
			ym = ym * query(mo(ar[cur]+1),mo(ar[cur+1]-1));
		if(cur+1==tot) break;
		tmpm.n[0][1] = br[cur+1], tmpm.n[1][1] = cr[cur+1];
		ym = ym*tmpm;
	}
	return ym.n[0][1];
}
int main()
{
	int i;
	k = readll(), p = read();
	if(k<2){ printf("%d\n",(int)k%p); return 0; }
	for(i=1,n=read();i<=n;i++) s[i] = read();
	for(S=1;S<=n+2;S<<=1);
	s[n+1] = s[1];
	ym.n[0][1] = 1; I.n[0][0] = I.n[1][1] = 1; psm[0] = ssm[n+1] = I;
	for(i=1;i<=n;i++)
	{
		tm[i].n[1][0] = 1;
		tm[i].n[0][1] = s[(i+n*3-2)%n+1];
		tm[i].n[1][1] = s[(i+n*3-1)%n+1];
		tr[S+i] = tm[i];
		psm[i] = psm[i-1] * tm[i];
	}
	for(i=n;i>=1;i--) ssm[i] = tm[i] * ssm[i+1];
	for(i=S;i>=1;i--) tr[i] = tr[i+i] * tr[i+i+1];
	for(i=1,m=read();i<=m;i++) q[i].first = readll(), q[i].second = read();
	q[0].first = q[m+1].first = -1;
	sort(q+1,q+1+m);
	printf("%d\n",calc());
	return 0;
}
posted @ 2015-11-23 19:35  MoebiusMeow  阅读(380)  评论(0)  编辑  收藏  举报