Description

​ SnukeSnuke从他的母亲那里得到了生日礼物——一个网格。网格有HH行WW列。每个单元格都是黑色或白色。所有黑色单元格都是四联通的,也就是说,只做水平或垂直移动且只经过黑色单元格即可从任何黑色单元格移动到任何其他黑色单元格。

第ii行第jj列的单元格的颜色由字符si,jsi,j表示。如果si,jsi,j是 #,该单元格为黑色;如果si,jsi,j是 .,该单元格为白色。至少一个单元格是黑色的。

我们定义「分形」如下:00级分形是一个 1×11×1的黑色单元格.kk级分形由HH行WW列较小一级的分形按照 SnukeSnuke 的网格的样式拼成:与SnukeSnuke 网格中的黑色单元格对应的位置是一个kk级分形;与SnukeSnuke 网格中的白色单元格对应的位置是一个单元格全部为白色,尺寸与kk级分形相同的网格。

您将得到 SnukeSnuke 的网格的描述和整数 KK。请求出KK级分形中黑色单元格组成的连通分量数,模109+7109+7。

Input

​ 第一行三个整数H,W,K如题目描述

​ 接下来H行,每行W个字符

Output

​ 输出KK级分形中黑色单元格组成的连通分量数,模109+7109+7。

Sample Input

3 3 3
.#.
###
#.#

Sample Output

20

HINT

本题采用subtask。

真的抱歉部分分一拖再拖……

  • 存在5%5%的数据满足黑色方块与边界不相连。
  • 存在10%10%的数据满足KK级分形的边长≤2000≤2000。
  • 存在50%50%的数据满足黑色方块与四边界相连。

对于所有数据:

  • 1≤H,W≤10001≤H,W≤1000
  • 0≤K≤10180≤K≤1018
  • 每个 si,jsi,j是 # 或 .
  • 网格中所有黑色单元格四联通
  • 网格中至少有一个黑色单元格
.............#.............
............###............
............#.#............
..........#..#..#..........
.........#########.........
.........#.##.##.#.........
..........#.....#..........
.........###...###.........
.........#.#...#.#.........
....#........#........#....
...###......###......###...
...#.#......#.#......#.#...
.#..#..#..#..#..#..#..#..#.
###########################
#.##.##.##.##.##.##.##.##.#
.#.....#..#.....#..#.....#.
###...######...######...###
#.#...#.##.#...#.##.#...#.#
....#.................#....
...###...............###...
...#.#...............#.#...
.#..#..#...........#..#..#.
#########.........#########
#.##.##.#.........#.##.##.#
.#.....#...........#.....#.
###...###.........###...###
#.#...#.#.........#.#...#.#

这个是第一个样例的样例解释

第二个样例:

3 3 3
###
#.#
### OUTPUT:1

第三个样例:

11 15 1000000000000000000
.....#.........
....###........
....####.......
...######......
...#######.....
..##.###.##....
..##########...
.###.....####..
.####...######.
###############
#.##..##..##..# OUTPUT:301811921

Sol

%%%DTZ

其实这题不是很难想,但是看起来太神了觉得不是很可做所以当时就索性放弃了QWQ

Code

#include <bits/stdc++.h>
using namespace std;
#define P 1000000007
int n,m,x,y,z;char s[1005][1005],t[1005][1005];long long k;
struct M
{
long long m[5][5];
M(){memset(m,0,sizeof(m));}
M operator*(const M &b)const
{M c;for(int i=1,j,k;i<=2;i++) for(j=1;j<=2;j++) for(k=1;k<=2;k++) (c.m[i][j]+=m[i][k]*b.m[k][j])%=P;return c;}
M operator^(long long x)const
{M a(*this),ans;for(int i=1;i<=2;i++)ans.m[i][i]=1;for(;x;x>>=1,a=a*a) if(x&1)ans=ans*a;return ans;}
}ans;
int ksm(int x,long long y){int ans=1;for(;y;y>>=1,x=1ll*x*x%P) if(y&1) ans=1ll*ans*x%P;return ans;}
bool chk(){for(int i=1;i<=n;i++) if(s[i][1]=='#'&&s[i][m]=='#') return 1;return 0;}
void turn90(){for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) t[m-j+1][i]=s[i][j];swap(n,m);memcpy(s,t,sizeof(s));}
int main()
{
scanf("%d%d%lld",&n,&m,&k);
for(int i=1,j;i<=n;i++) for(scanf("%s",s[i]+1),j=1;j<=m;j++) x+=s[i][j]=='#';
bool a=chk(),b=(turn90(),chk());turn90();
if(a&&b) return printf("1"),0;
else if(!a&&!b) return printf("%d",ksm(x,k-1)),0;
else if(b) turn90();
for(int i=1;i<=n;i++) z+=s[i][1]=='#'&&s[i][m]=='#';
for(int i=1;i<=n;i++) for(int j=1;j<m;j++) y+=s[i][j]=='#'&&s[i][j+1]=='#';
ans.m[1][1]=x,ans.m[1][2]=y,ans.m[2][2]=z;ans=ans^(k-1);
printf("%lld",(ans.m[1][1]-ans.m[1][2]+P)%P);
}

[AGC003F] Fraction of Fractal 矩阵快速幂的更多相关文章

  1. 题解-AtCoder-agc003F Fraction of Fractal(非矩阵快速幂解法)

    Problem AtCoder-agc003F 题意:给出\(n\)行\(m\)列的01矩阵,一开始所有 \(1\) 连通,称此为\(1\)级分形,定义\(i\)级分形为\(i-1\)级分形中每个标示 ...

  2. [AGC003F] Fraction of Fractal(矩阵乘法)

    Description Snuke从他的母亲那里得到了生日礼物--一个网格.网格有H行W列.每个单元格都是黑色或白色.所有黑色单元格都是四联通的,也就是说,只做水平或垂直移动且只经过黑色单元格即可从任 ...

  3. 矩阵快速幂 HDU 4565 So Easy!(简单?才怪!)

    题目链接 题意: 思路: 直接拿别人的图,自己写太麻烦了~ 然后就可以用矩阵快速幂套模板求递推式啦~ 另外: 这题想不到或者不会矩阵快速幂,根本没法做,还是2013年长沙邀请赛水题,也是2008年Go ...

  4. 51nod 算法马拉松18 B 非010串 矩阵快速幂

    非010串 基准时间限制:1 秒 空间限制:131072 KB 分值: 80 如果一个01字符串满足不存在010这样的子串,那么称它为非010串. 求长度为n的非010串的个数.(对1e9+7取模) ...

  5. 51nod 1113 矩阵快速幂

    题目链接:51nod 1113 矩阵快速幂 模板题,学习下. #include<cstdio> #include<cmath> #include<cstring> ...

  6. 【66测试20161115】【树】【DP_LIS】【SPFA】【同余最短路】【递推】【矩阵快速幂】

    还有3天,今天考试又崩了.状态还没有调整过来... 第一题:小L的二叉树 勤奋又善于思考的小L接触了信息学竞赛,开始的学习十分顺利.但是,小L对数据结构的掌握实在十分渣渣.所以,小L当时卡在了二叉树. ...

  7. HDU5950(矩阵快速幂)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5950 题意:f(n) = f(n-1) + 2*f(n-2) + n^4,f(1) = a , f(2 ...

  8. 51nod 1126 矩阵快速幂 水

    有一个序列是这样定义的:f(1) = 1, f(2) = 1, f(n) = (A * f(n - 1) + B * f(n - 2)) mod 7. 给出A,B和N,求f(n)的值. Input 输 ...

  9. hdu2604(递推,矩阵快速幂)

    题目链接:hdu2604 这题重要的递推公式,找到公式就很easy了(这道题和hdu1757(题解)类似,只是这道题需要自己推公式) 可以直接找规律,推出递推公式,也有另一种找递推公式的方法:(PS: ...

随机推荐

  1. SQLSERVER2012误删数据恢复过程

    由于长时间从事企业应用系统开发,前往用户现场升级.调试系统是比较常做的事情,但是就在周一,由于同事的失误在毫无知觉的情况下误删了生产数据库几乎所有的数据.当我发现的那一刻,感觉头发都立起来了,心想这他 ...

  2. Gym101128G:Game of Cards

    题意: 有P摞纸牌和一个数字k,每次可以从一摞中拿0-k张牌,拿完再剩下的牌中的第一张数字是几,就必须再拿几张,谁不能拿谁输. emmm感觉好像就是裸的SG游戏啊,数据不大,递推出每一摞牌的SG值,然 ...

  3. impdp报错ORA-39083 ORA-02304 Object type TYPE failed to create

    环境Red Hat Enterprise Linux Server release 5.8 (Tikanga)ORACLE Release 11.2.0.3.0 Production 我用expdp, ...

  4. RTX二次开发::检索 COM 类工厂中 CLSID 为 {79210E58-99EB-45D0-8890-763EFEAA143F} 的组件失败,

    解决方法: 1.Interop.RTXSAPILib.dll这个是32位的 将IIS 应用程序池 启用32位应用程序 设置为true就可以了, 2.把编译好的程序选择X86平台运行,就好了 程序员的基 ...

  5. Nginx概述、安装及配置详解

    nginx概述 nginx是一款自由的.开源的.高性能的HTTP服务器和反向代理服务器:同时也是一个IMAP.POP3.SMTP代理服务器:nginx可以作为一个HTTP服务器进行网站的发布处理,另外 ...

  6. js正则表达式语法[转]

    1. 正则表达式规则 1.1 普通字符 字母.数字.汉字.下划线.以及后边章节中没有特殊定义的标点符号,都是"普通字符".表达式中的普通字符,在匹配一个字符串的时候,匹配与之相同的 ...

  7. source insight 保存时删除多余空格,去除多余空格 space tab键

    source insight 保存时删除多余空格,去除多余空格 space tab键 摘自:https://blog.csdn.net/lanmanck/article/details/8638391 ...

  8. 3.3.5 高效读取:不变模式下的CopyOnWriteArrayList

    源码分析:读写(get,add) 一:get 方法 private E get(Object[] a, int index) { return (E) a[index];}可以看到读取数据的时候 没有 ...

  9. [GO]通过结构体生成json

    package main import ( "encoding/json" "fmt" ) type IT struct { //一定要注意这里的成员变量的名字 ...

  10. CentOS 7.2配置Apache服务httpd(上)

    http://www.jb51.net/article/97434.htm 二.安装Apache httpd 安装httpd以配置Web服务器, HTTP使用80 / TCP ? 1 2 3 4 5 ...