【bzoj5133】[CodePlus2017年12月]白金元首与独舞 并查集+矩阵树定理
题目描述
给定一个 $n\times m$ 的方格图,每个格子有 ↑、↓、←、→,表示从该格子能够走到相邻的哪个格子。
有一些格子是空着的,需要填上四者之一,需要满足:最终的方格图中,从任意一个位置出发都能够走出方格图。求方案数 mod 10^9+7。
$数据组数\le 10$ ,$n,m\le 300$ ,$空格子数k\le 200$
题解
并查集+矩阵树定理
由于k很小,又是计数问题,考虑矩阵树定理。
先使用并查集处理出从每个位置开始,最终会走到哪个位置。显然如果有环则答案为0,否则一定走到的是一个空格子或方格图外部。
这样就不用考虑已填好的格子的走法,只需要考虑空格子的走法即可。
每个空格子需要走到方格图外部,不能有环,相当于是一棵以方格图外部为根的内向树形图。
考虑每个空格子4个方向会走到哪个空格子(或外部),连边,矩阵树定理求解即可。
本题要求的是内向树,因此求 出度矩阵-邻接矩阵 删去根节点所在行列,得到的行列式的值 即可。
时间复杂度 $O(nm+k^3)$
#include <cstdio>
#include <cstring>
#include <algorithm>
#define mod 1000000007
using namespace std;
typedef long long ll;
int id[210][210] , f[40010] , flag , v[40010] , wx[310] , wy[310];
ll a[310][310];
char str[210];
inline ll pow(ll x , int y)
{
ll ans = 1;
while(y)
{
if(y & 1) ans = ans * x % mod;
x = x * x % mod , y >>= 1;
}
return ans;
}
int find(int x)
{
return x == f[x] ? x : f[x] = find(f[x]);
}
inline void link(int x , int y)
{
x = find(x) , y = find(y);
if(x == y) flag = 1;
f[x] = y;
}
int main()
{
int T;
scanf("%d" , &T);
while(T -- )
{
int n , m , p = 0 , i , j , k , d = 0;
ll t , ans = 1;
scanf("%d%d" , &n , &m);
memset(id , 0 , sizeof(id)) , flag = 0;
for(i = 1 ; i <= n ; i ++ )
for(j = 1 ; j <= m ; j ++ )
id[i][j] = (i - 1) * m + j;
for(i = 0 ; i <= n * m ; i ++ ) f[i] = i;
for(i = 1 ; i <= n ; i ++ )
{
scanf("%s" , str + 1);
for(j = 1 ; j <= m ; j ++ )
{
switch(str[j])
{
case 'L': link(id[i][j] , id[i][j - 1]); break;
case 'R': link(id[i][j] , id[i][j + 1]); break;
case 'U': link(id[i][j] , id[i - 1][j]); break;
case 'D': link(id[i][j] , id[i + 1][j]); break;
default: v[id[i][j]] = ++p , wx[p] = i , wy[p] = j;
}
}
}
if(flag) puts("0");
else
{
memset(a , 0 , sizeof(a));
for(i = 1 ; i <= p ; i ++ )
{
a[i][i] += 4;
a[i][v[find(id[wx[i]][wy[i] - 1])]] -- ;
a[i][v[find(id[wx[i]][wy[i] + 1])]] -- ;
a[i][v[find(id[wx[i] - 1][wy[i]])]] -- ;
a[i][v[find(id[wx[i] + 1][wy[i]])]] -- ;
}
for(i = 1 ; i <= p ; i ++ )
for(j = 1 ; j <= p ; j ++ )
a[i][j] = (a[i][j] + mod) % mod;
for(i = 1 ; i <= p ; i ++ )
{
for(j = i ; j <= p ; j ++ )
if(a[i][j])
break;
if(j > p) continue;
if(j != i)
{
d ^= 1;
for(k = i ; k <= p ; k ++ )
swap(a[i][k] , a[j][k]);
}
ans = ans * a[i][i] % mod;
t = pow(a[i][i] , mod - 2);
for(j = i ; j <= p ; j ++ ) a[i][j] = a[i][j] * t % mod;
for(j = i + 1 ; j <= p ; j ++ )
for(t = a[j][i] , k = i ; k <= p ; k ++ )
a[j][k] = (a[j][k] - a[i][k] * t % mod + mod) % mod;
}
for(i = 1 ; i <= p ; i ++ ) ans = ans * a[i][i] % mod;
if(d) ans = (mod - ans) % mod;
printf("%lld\n" , ans);
}
}
return 0;
}
【bzoj5133】[CodePlus2017年12月]白金元首与独舞 并查集+矩阵树定理的更多相关文章
- [BZOJ5133][CodePlus2017年12月]白金元首与独舞
bzoj luogu 题意 给你一个\(n*m\)的网格,每个位置上有一个箭头指向上或下或左或右.有些位置上还没有箭头,现在要求你在这些没有箭头的位置上填入箭头,使得从网格的任意一个位置开始,都可以沿 ...
- 【BZOJ5133】[CodePlus2017年12月]白金元首与独舞 矩阵树定理
[BZOJ5133][CodePlus2017年12月]白金元首与独舞 题面:www.lydsy.com/JudgeOnline/upload/201712/div1.pdf 题解:由于k很小,考虑用 ...
- [LOJ#6259]「CodePlus 2017 12 月赛」白金元首与独舞
[LOJ#6259]「CodePlus 2017 12 月赛」白金元首与独舞 试题描述 到河北省 见斯大林 / 在月光下 你的背影 / 让我们一起跳舞吧 うそだよ~ 河北省怎么可能有 Stalin. ...
- BZOJ5131: [CodePlus2017年12月]可做题2
BZOJ没有题面,差评 洛谷的题目链接 题解 其实这题很久之前就写了,也想写个题解但是太懒了,咕到了今天 在typora写完题解不想copy过来再改格式了,于是直接贴截图qwq #include &l ...
- 【LibreOJ】#6259. 「CodePlus 2017 12 月赛」白金元首与独舞
[题目]给定n行m列的矩阵,每个位置有一个指示方向(上下左右)或没有指示方向(任意选择),要求给未定格(没有指示方向的位置)确定方向,使得从任意一个开始走都可以都出矩阵,求方案数.n,m<=20 ...
- 「CodePlus 2017 12 月赛」白金元首与独舞
description 题面 data range \[ 1 \leq T \leq 10, 1 \leq n, m \leq 200 , 0 \leq k \leq \min(nm, 300)\] ...
- 走进矩阵树定理--「CodePlus 2017 12 月赛」白金元首与独舞
n,m<=200,n*m的方阵,有ULRD表示在这个格子时下一步要走到哪里,有一些待决策的格子用.表示,可以填ULRD任意一个,问有多少种填法使得从每个格子出发都能走出这个方阵,答案取模.保证未 ...
- loj6259「CodePlus 2017 12 月赛」白金元首与独舞
分析 我们将没连的点连向周围四个点 其余的按照给定的方向连 我们将所有连出去的位置统一连到0点上 再以0作为树根 于是就将问题转化为了有向图内向树计数 代码 #include<iostream& ...
- Solution -「Code+#2」「洛谷 P4033」白金元首与独舞
\(\mathcal{Description}\) link. 给定一个 \(n\times m\) 的网格图,一些格子指定了走出该格的方向(上下左右),而有 \(k\) 格可以任意指定走出方 ...
随机推荐
- 2016-2017-2015329 《Java程序设计》第4周学习总结
学号 2016-2017-2015329 <Java程序设计>第4周学习总结 教材学习内容总结 面向对象有三大特性:封装.继承.多态 封装 封装是指,一种将抽象性函式接口的实例细节部份包装 ...
- day 2 异常传递 ,抛出
1.异常的传递 def test1(): print("---test1--") print(num) print('---test1 over---') def test2(): ...
- Drupal学习(19) 使用jQuery
本节学习如果在Drupal里交互使用jQuery. jQuery在Drupal是内置支持的.存在根目录的misc目录中. 当调用drupal_add_js方法,会自动加载jQuery. 在Drupal ...
- bootstrap的Alerts中 可以放置p标签 设置 align="center" 用来设置文本居中
效果
- NOIP2018出征策
蒟蒻的风之旅人即将退役,现在分享一下退休前的故事 首先,经过这么多时间的划水训练,我成功从一个萌新变成了一个蒟蒻.我学会了各种奇怪玄学的算法,比如说昨天老师讲的NOIP第三题通用的算法,叫做XG算法, ...
- python的30个编程技巧
1.原地交换两个数字 x, y =10, 20 print(x, y) y, x = x, y print(x, y) 10 20 20 10 2.链状比较操作符 n = 10 print(1 &l ...
- 在Office 365 的如何给管理员赋予查看所有人邮箱的权限的Powershell
连接至Office365 的Powershell Get-MsolUser -UserPrincipalName admin@***.partner.onmschina.cn //Get-MsolUs ...
- rest_framework组件
认证组件 局部认证 在需要认证的视图类里加上authentication_classes = [认证组件1类名,认证组件2类名....] 示例如下: seralizers.py from rest_f ...
- Vue工作原理小结
本文能帮你做什么?1.了解vue的双向数据绑定原理以及核心代码模块2.缓解好奇心的同时了解如何实现双向绑定为了便于说明原理与实现,本文相关代码主要摘自vue源码, 并进行了简化改造,相对较简陋,并未考 ...
- Python中import的as语法
在Python中,如果import的语句比较长,导致后续引用不方便,可以使用as语法,比如: import dir1.dir2.mod # 那么,后续对mod的引用,都必须是dir1.dir2.mod ...