[agc003F]Fraction of Fractal

Description

传送门

Solution

本篇博客思路来自大佬的博客（侵删）。

我们定义如果网格的第一行和最后一行的第i列都为黑色，则它是一个上下界接口。左右界接口定义同上。

如果上下界接口和左右界接口都为0个，则答案为节点数^(k-1)。或者上下界接口和左右界接口都存在，则答案为1（就只有一个联通块了啊）。

将以上两种情况特判后，我们目前的网格要么只有上下界接口，要么只有左右界接口。我们只探讨前者。

对于k=2的分形（如果k<2特判就好），我们将1级分形看作一个节点（以下的点都是指1级分形，边同理），如果它与其他节点连通则视为将两个节点之间加一条边。

因为我们的网格只有上下界接口，所以2级分形事实上是由若干条链构成。以此类推，k级分形也是如此。

将链看成树，我们就可以运用森林的性质：树的个数等于总点数减去总边数。（这个性质是真的牛啊我自己完全想不到），即最后联通块的个数为总点数减去总边数。

考虑求总边数：

　　记原图中节点总数为all。

　　定义若a[i][j]与a[i+1][j]都为黑点则这是一个上下接口，（注意不是上下界接口），记其个数为usb。

　　记上下界接口个数为u。

　　设E_k为k级分形的边数。

　　　　边界E₂=usb   递推式为 Ek=E_k-1*all+usb*u^k-1

至于总点数（这里的点还是指1级分形）。。快速幂就好了。

Code

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
const int mod=1e9+;
typedef long long ll;
char ch[];
int a[][],h,w,up_down,left_right,cnt,usb;//cnt:黑点总数，usb：原图内总接口数
long long _pow;
ll ksm(ll x,long long k)
{
    ll ans=;
    while (k)
    {
        if (k&)ans=ans*x%mod;
        k>>=;
        x=x*x%mod;
    }
    return ans;
}
struct Matrix
{
    ll t[][];
    friend Matrix operator *(Matrix x,Matrix y)
    {
        Matrix z;
        memset(z.t,,sizeof(z.t));
        for (int i=;i<;i++)
            for (int j=;j<;j++)
                for (int k=;k<;k++)
                    z.t[i][j]=(z.t[i][j]+x.t[i][k]*y.t[k][j])%mod;
        return z;
    }
    friend Matrix operator ^(Matrix x,long long k)
    {
        Matrix ans;
        memset(ans.t,,sizeof(ans.t));
        ans.t[][]=ans.t[][]=;
        while (k)
        {
            if (k&) ans=ans*x;
            k>>=;
            x=x*x;
        }
        return ans;
    }
}A,B;
int main()
{
    scanf("%d%d%lld",&h,&w,&_pow);
    if (_pow<=) {printf("");return ;}
    for (int i=;i<=h;i++)
    {
        scanf("%s",ch+);
        for (int j=;j<=w;j++)
        {
            a[i][j]=(ch[j]=='#');
            if (a[i][j]) cnt++;
        }
    }
    for (int i=;i<=w;i++)
        if (a[][i]&&a[h][i]) up_down++;
    for (int i=;i<=h;i++)
        if (a[i][]&&a[i][w]) left_right++;
    if (!up_down&&!left_right)
    {
        printf("%lld",ksm(cnt,_pow-));return ;
    }
    if (up_down&&left_right)
    {
        printf("");return ;
    }
    for (int i=;i<=h;i++)
        for (int j=;j<=w;j++)
        {
            if (a[i][j]&&a[i][j+]&&left_right) usb++;
            if (a[i][j]&&a[i+][j]&&up_down) usb++;
        }
    A.t[][]=cnt;A.t[][]=usb;A.t[][]=max(up_down,left_right);
    A=A^(_pow-);
    ll ans;ans=A.t[][]*usb%mod+A.t[][]*max(up_down,left_right)%mod;
    ans%=mod;
    printf("%lld",(ksm(cnt,_pow-)-ans+mod)%mod);
}