若$\min(n,m)<\min(n+1,m-1)$,则考虑计算左边与右边不连通的概率,然后用$1$减去它得到答案。

若$\min(n,m)\geq \min(n+1,m-1)$,则考虑计算对偶图中上边与下边不连通的概率。

同时当$n<m$时还可以旋转$90$°来交换$n$和$m$,使得$m\leq 9$。

那么问题转化为,给定$n\times m$的网格图,每个点颜色为$0,1,2$,每条边有存在的概率,求每个连通块不同时含有颜色$1$和$2$的概率。

连通性状压DP,设$f[i][j][S]$表示考虑到了$(i,j)$,轮廓线状态为$S$的概率。

其中$S$记录轮廓线上$m$个点所属连通块的最小表示,以及每个连通块的颜色。

为了避免map的使用,可以考虑扩展出所有状态后再进行排序,将相同状态的DP值合并,减小常数。

同时可以将$1$和$2$内部的边存在的概率设置为$1$,减小连通块的数量已达到优化状态数的效果。

#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=115,M=1000000,P=1000000007;
int n,m,o,i,j,k,_id[N][N],_wr[N][N],_wd[N][N],id[N][N],wr[N][N],wd[N][N];
struct E{
ll x;int y,v;
E(){}
E(ll _x,int _y,int _v){x=_x,y=_y,v=_v;}
}e[2][M];
int cnt[2];
inline bool cmp(const E&a,const E&b){
if(a.x!=b.x)return a.x<b.x;
return a.y<b.y;
}
inline void up(int&a,int b){a=a+b<P?a+b:a+b-P;}
inline int po(int a,int b){int t=1;for(;b;b>>=1,a=1LL*a*a%P)if(b&1)t=1LL*t*a%P;return t;}
inline int read(){
int x;
scanf("%d",&x);
return 1LL*x*po(100,P-2)%P;
}
inline int read2(){
int x=read();
x=(P+1-x)%P;
return x;
}
inline void fix(){
int m=cnt[o],t=0,i,j,k;
sort(e[o],e[o]+m,cmp);
for(i=0;i<m;i=j){
for(k=0,j=i;j<m&&e[o][i].x==e[o][j].x&&e[o][i].y==e[o][j].y;j++)up(k,e[o][j].v);
if(k)e[o^1][t++]=E(e[o][i].x,e[o][i].y,k);
}
cnt[o^=1]=t;
}
inline ll get(ll x,int y){return x>>(y<<2)&15;}
inline void cal(ll&A,int&B){
static int v[N],f[N];
int i,x,t=0,C=0;
for(i=0;i<=m;i++)v[i]=-1;
for(t=i=0;i<m;i++){
x=get(A,i);
if(v[x]==-1){
C|=(B>>(x<<1)&3)<<(t<<1);
v[x]=t++;
}
f[i]=v[x];
}
for(A=0,B=C,i=m-1;~i;i--)A=A<<4|f[i];
}
inline void up0(ll A,int B,int C,int x,int y){
A^=(get(A,x)^m)<<(x<<2);
B|=y<<(m<<1);
cal(A,B);
e[o^1][cnt[o^1]++]=E(A,B,C);
}
inline void upl(ll A,int B,int C,int x,int y){
int t=get(A,x-1);
if(y+(B>>(t<<1)&3)==3)return;
A^=(get(A,x)^t)<<(x<<2);
B|=y<<(t<<1);
cal(A,B);
e[o^1][cnt[o^1]++]=E(A,B,C);
}
inline void upu(ll A,int B,int C,int x,int y){
int t=get(A,x);
if(y+(B>>(t<<1)&3)==3)return;
B|=y<<(t<<1);
e[o^1][cnt[o^1]++]=E(A,B,C);
}
inline void upul(ll A,int B,int C,int x,int y){
int t=get(A,x),z=get(A,x-1);
if(y+(B>>(t<<1)&3)==3)return;
B|=y<<(t<<1);
if((B>>(z<<1)&3)+(B>>(t<<1)&3)==3)return;
B|=(B>>(z<<1)&3)<<(t<<1);
ll S=z^t;
for(int i=0;i<m;i++)if(get(A,i)==z)A^=S<<(i<<2);
cal(A,B);
e[o^1][cnt[o^1]++]=E(A,B,C);
}
int solve(){
if(n<m){
for(i=1;i<=n;i++)for(j=1;j<m;j++)wd[j][i]=_wr[i][j];
for(i=1;i<n;i++)for(j=1;j<=m;j++)wr[j][i]=_wd[i][j];
for(i=1;i<=n;i++)for(j=1;j<=m;j++)id[j][i]=_id[i][j];
swap(n,m);
}else{
for(i=1;i<=n;i++)for(j=1;j<=m;j++)wd[i][j]=_wd[i][j],wr[i][j]=_wr[i][j],id[i][j]=_id[i][j];
}
cnt[o=0]=1;
e[0][0]=E(0,0,1);
for(i=m-1;~i;i--)e[0][0].x=e[0][0].x<<4|i;
for(i=1;i<=n;i++)for(j=1;j<=m;j++){
int pl=j>1?wr[i][j-1]:0,npl=(P+1-pl)%P,
pu=i>1?wd[i-1][j]:0,npu=(P+1-pu)%P,
pul=1LL*pl*pu%P,p0=1LL*npl*npu%P,t=id[i][j];
pl=1LL*pl*npu%P;
pu=1LL*pu*npl%P;
cnt[o^1]=0;
for(k=0;k<cnt[o];k++){
if(pl)upl(e[o][k].x,e[o][k].y,1LL*e[o][k].v*pl%P,j-1,t);
if(pu)upu(e[o][k].x,e[o][k].y,1LL*e[o][k].v*pu%P,j-1,t);
if(pul)upul(e[o][k].x,e[o][k].y,1LL*e[o][k].v*pul%P,j-1,t);
if(p0)up0(e[o][k].x,e[o][k].y,1LL*e[o][k].v*p0%P,j-1,t);
}
o^=1;
fix();
}
int ret=0;
for(i=0;i<cnt[o];i++)up(ret,e[o][i].v);
return ret;
}
int main(){
int r,c,flag=0;
scanf("%d%d",&r,&c);
if(c<=1)return puts("1"),0;
if(min(r,c)<min(r+1,c-1)){
flag=1;
for(i=1;i<=r;i++)for(j=1;j<c;j++)_wr[i][j]=read();
for(i=1;i<r;i++)for(j=2;j<c;j++)_wd[i][j]=read();
for(i=1;i<r;i++)_wd[i][1]=_wd[i][c]=1;
for(i=1;i<=r;i++)_id[i][1]=1,_id[i][c]=2;
n=r,m=c;
}else{
for(i=1;i<=r;i++)for(j=1;j<c;j++)_wd[i][j]=read2();
for(i=1;i<r;i++)for(j=2;j<c;j++)_wr[i+1][j-1]=read2();
for(i=1;i<c-1;i++)_wr[1][i]=_wr[r+1][i]=1;
for(i=1;i<c;i++)_id[1][i]=1,_id[r+1][i]=2;
n=r+1,m=c-1;
}
int ans=solve();
if(flag)ans=(P+1-ans)%P;
return printf("%d",ans),0;
}

  

BZOJ5207 : [Jsoi2017]隧道的更多相关文章

  1. SSH隧道应用, 突破网络限制

    文/怡文圣美 这篇文章可以帮你解决下面三个问题: 不能直连服务器, 要先登陆跳板机, 造成客户端工具无法连接服务器. 服务器没有公网IP, 且只允许公司IP访问, 要在家里操作要先远程桌面到工作机再登 ...

  2. 前端学HTTP之网关、隧道和中继

    前面的话 Web是一种强大的内容发布工具.人们已经从只在网上发送静态的在线文档,发展到共享更复杂的资源,比如数据库内容或动态生成的HTML页面.Web浏览器为用户提供了一种统一的方式来访问因特网上的内 ...

  3. ngrok反向隧道--获取内网IP

    ngrok反向隧道 前情提要:小明与小白各有一台主机,两台主机在同一内网,小明想直接通过内网ssh到小白的主机上.但是小白的ip地址会不断的变化,而小明不想每次都要麻烦小白查看ip.于是小明催生了一个 ...

  4. SSH正反向隧道

    正向隧道 拓扑如下: 说明: CLIENT不能直接访问WEB服务器,AGENT可访问WEB服务器: 在AGENT上通过创建ssh正向隧道,使CLIENT可以通过AGENT间接访问WEB服务器: AGE ...

  5. Windows下使用Xshell建立反向隧道

    反向隧道是一个进行内网穿透的简单而有用的方法.在Linux下通过OpenSSH和AutoSSH可以很容易地建立稳定的反向隧道.但是在Windows下,还能看到有人特意装个Cygwin来运行这些工具…… ...

  6. WPF的路由事件、冒泡事件、隧道事件(预览事件)

    本文摘要: 1:什么是路由事件: 2:中断事件路由: 3:自定义路由事件: 4:为什么需要自定义路由事件: 5:什么是冒泡事件和预览事件(隧道事件): 1:什么是路由事件 WPF中的事件为路由事件,所 ...

  7. XSHELL使用隧道

    线上系统中,搭建了一个elasticsearch环境,想要访问页面,发现环境的内网中没有windows机器,无法使用浏览器来直接进行web页面的访问,于是直接使用了XSELL中强大的功能"隧 ...

  8. HTTP权威协议笔记-7.集成点:网关、隧道及中继

    .8.1 网关 定义:网关类似与翻译器,它抽象出了一种能够到达资源的方法. 实用:网关可以自动将HTTP流量转换为其他协议,这样使用HTTP协议的一方就不需要了解其他协议,也可实现与其他程序或设备交互 ...

  9. Android终端配置isatap隧道使用IPV6的方法

    使用isatap隧道可以在手机有IPV4网络的情况下访问IPv6网络资源.关于isatap隧道的配置方法,清华.上交两所学校都给出了相应的在windows xp/2003,windows 2000,以 ...

随机推荐

  1. Linux下安装软件命令详解

    ---------------------------------------------------------------- 或许你对于linux还不够了解,但是一旦你步入公司后,你就会发现lin ...

  2. HDU5950 Recursive sequence (矩阵快速幂)

    Recursive sequence Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Other ...

  3. python练习册 0002随机生成验证

    这个题需要用到random库的方法,不用就会忘,暂把random库的常用方法贴出来 import random import string # 随机整数 # randint(a, b),生成a~b之间 ...

  4. 使用android-ndk官方ndkbuild例子

    Why this blog 现在(2018年9月27日),Android Studio中新建ndk项目都使用cmake而不是Android.mk+Application.mk的方式.但老项目需要维护, ...

  5. error: not found: value sqlContext/import sqlContext.implicits._/error: not found: value sqlContext /import sqlContext.sql/Caused by: java.net.ConnectException: Connection refused

    1.今天启动启动spark的spark-shell命令的时候报下面的错误,百度了很多,也没解决问题,最后想着是不是没有启动hadoop集群的问题 ,可是之前启动spark-shell命令是不用启动ha ...

  6. JSP基础知识➣Cookie和Session(五)

    JSP Cookie 处理 Cookie是存储在客户机的文本文件,它们保存了大量轨迹信息.在servlet技术基础上,JSP显然能够提供对HTTP cookie的支持.JSP脚本通过request对象 ...

  7. bzoj2870

    题解: 边分治入门题 当然并查集+维护直径更加简单 就是两棵树二合一直径是两颗树上的4个直径两两组合的最大值 查询路径长度你搞个差分查个lca就好了 点分治并不能做这题 分成多个联通块就gg了(点分治 ...

  8. PHP Manager 安装失败的解决方法, PHP Manager 1.4 for IIS 10,经验证支持windows server 2016版本

    // 另有无需进行修改注册表的安装包,经测试最高支持Windows Server 2016 版本,下载地址如下: https://github.com/EnhWeb/PHPManager/tree/m ...

  9. Flink在流处理上常见的Source和sink操作

    flink在流处理上的source和在批处理上的source基本一致.大致有4大类 1.基于本地集合的source(Collection-based-source) 2.基于文件的source(Fil ...

  10. 2018牛客网暑假ACM多校训练赛(第五场)H subseq 树状数组

    原文链接https://www.cnblogs.com/zhouzhendong/p/NowCoder-2018-Summer-Round5-H.html 题目传送门 - https://www.no ...