题目链接:
题意:
给你一个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);
}