[AGC003F] Fraction of Fractal

题目链接：

AGC003F

题意：

给你一个n*m的网格，一个k，每个格子要么就是黑的要么就是白的，要让你求如果用这个图形构成k阶分形，求联通块数量。

题解：

一看这个k的范围，，就肯定是矩阵乘法的题了，可惜场上没有推出来。

首先如果这个图形与四周不连通，那么答案就一定是，其中cnt为黑点个数。

如果这个图形上接下是联通的，左接右也是联通的，那么答案就一定是1，因为最后的分形一定联通。

那么剩下的情况就只有上边和下边联通，左边和右边联通两种情况之一了。

那么我们就只需求出cnt，a，b，其中a表示多少个黑点满足在这行与它相邻的下一个点也是黑的，b表示有多少行是连通的。

构造的矩阵长这样：

用快速幂求，答案就是。

代码：

#include <bits/stdc++.h>
using namespace std;

#define mod 1000000007

int n,m,cnt,flag1,flag2,s1,s2;
long long k;
char s[2001][2001];

long long quickpow(long long base,long long to)
{
	if(to<=0)return 1;
	if(to==1)return base;
	long long mid=quickpow(base,to>>1);
	if(to&1)return mid*mid%mod*base%mod;
	else return mid*mid%mod;
}

struct Matrix{
	long long a[3][3];
	Matrix(){memset(a,0,sizeof(a));}
	Matrix operator*(Matrix &m)
	{
		Matrix ans;
		for(int i=1;i<=2;++i)
			for(int j=1;j<=2;++j)
				for(int k=1;k<=2;++k)
					ans.a[i][j]=(ans.a[i][j]+a[i][k]*m.a[k][j]%mod)%mod;
		return ans;
	}
};

int main()
{
	scanf("%d%d%lld",&n,&m,&k);
	for(int i=1; i<=n; ++i)
		scanf("%s",s[i]+1);
	for(int i=1; i<=n; ++i)
	{
		for(int j=1; j<=m; ++j)
			if(s[i][j]=='#')
			{
				++cnt;
				s1+=s[i][j-1]=='#';
				s2+=s[i-1][j]=='#';
			}
		if(s[i][1]=='#' && s[i][m]=='#')++flag1;
	}
	for(int i=1;i<=m;++i)if(s[1][i]=='#' && s[n][i]=='#')++flag2;
	if(flag1 && flag2)
	{
		puts("1");
		return 0;
	}
	if(!flag1 && !flag2)
	{
		printf("%lld\n",quickpow(cnt,k-1));
		return 0;
	}
	if(!flag1)swap(flag1,flag2),swap(s1,s2);
	Matrix A;
	A.a[1][1]=cnt;
	A.a[2][1]=s1;
	A.a[2][2]=flag1;
	--k;
	Matrix ANS;
	ANS.a[1][1]=ANS.a[2][2]=1;
	while(k)
	{
		if(k&1)ANS=ANS*A;
		A=A*A;
		k>>=1;
	}
	printf("%lld\n",(ANS.a[1][1]-ANS.a[2][1]+mod)%mod);
}