【BZOJ3814】【清华集训2014】简单回路 状压DP
题目描述
给你一个\(n\times m\)的网格图和\(k\)个障碍,有\(q\)个询问,每次问你有多少个不同的不经过任何一个障碍点且经过\((x,y)\)与\((x+1,y)\)之间的简单回路
\(n\leq 1000,m\leq 6,q\leq 10000\)
题解
简单插头DP
先用DP求出前面\(i\)行的的方案数和后面\(i\)行的方案数,询问时暴力合并。
其实状态数很少的(不到\(200\))
时间复杂度:\(O(???)\)
代码
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<ctime>
#include<utility>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
ll p=1000000007;
int c[100010];
int d[1010];
int ld[1010][10];
int rd[1010][10];
int cnt;
int n,m,k,q;
int la[10];
int ra[10];
int set(int s,int x,int v)
{
s&=~(3<<(2*(x-1)));
s|=v<<(2*(x-1));
return s;
}
int get(int s,int x)
{
return (s>>(2*(x-1)))&3;
}
int st[10];
void dfs(int x,int s)
{
if(x>m+1)
{
int i;
int now=0;
for(i=1;i<=m+1;i++)
{
if(get(s,i)==1)
st[++now]=i;
else if(get(s,i)==2)
{
la[i]=st[now];
ra[st[now]]=i;
now--;
}
if(now<0)
return;
}
if(now>0)
return;
for(i=1;i<=m+1;i++)
fprintf(stderr,"%d",get(s,i));
fprintf(stderr,"\n");
d[++cnt]=s;
c[s]=cnt;
memcpy(ld[cnt],la,sizeof la);
memcpy(rd[cnt],ra,sizeof ra);
return;
}
int i;
for(i=0;i<=2;i++)
dfs(x+1,set(s,x,i));
}
void add(ll &x,ll y)
{
x=(x+y)%p;
}
struct dp
{
int a[1010][10];
ll f[1010][10][200];
void solve()
{
f[1][1][1]=1;
int i,j,k;
for(i=1;i<=n;i++)
{
for(j=1;j<=m;j++)
for(k=1;k<=cnt;k++)
{
if(!f[i][j][k])
continue;
int s=d[k];
int l=get(s,j);
int r=get(s,j+1);
if(a[i][j])
{
if(!l&&!r)
add(f[i][j+1][k],f[i][j][k]);
continue;
}
if(!l)
if(!r)
{
add(f[i][j+1][k],f[i][j][k]);
int v=set(s,j,1);
v=set(v,j+1,2);
add(f[i][j+1][c[v]],f[i][j][k]);
}
else if(r==1)
{
add(f[i][j+1][k],f[i][j][k]);
int v=set(s,j,1);
v=set(v,j+1,0);
add(f[i][j+1][c[v]],f[i][j][k]);
}
else
{
add(f[i][j+1][k],f[i][j][k]);
int v=set(s,j,2);
v=set(v,j+1,0);
add(f[i][j+1][c[v]],f[i][j][k]);
}
else if(l==1)
{
if(!r)
{
add(f[i][j+1][k],f[i][j][k]);
int v=set(s,j,0);
v=set(v,j+1,1);
add(f[i][j+1][c[v]],f[i][j][k]);
}
else if(r==1)
{
int v=set(s,j,0);
v=set(v,j+1,0);
v=set(v,rd[k][j+1],1);
add(f[i][j+1][c[v]],f[i][j][k]);
}
else
{
// fprintf(stderr,"orzzjt\n");
}
}
else
{
if(!r)
{
add(f[i][j+1][k],f[i][j][k]);
int v=set(s,j,0);
v=set(v,j+1,2);
add(f[i][j+1][c[v]],f[i][j][k]);
}
else if(r==1)
{
int v=set(s,j,0);
v=set(v,j+1,0);
add(f[i][j+1][c[v]],f[i][j][k]);
}
else
{
int v=set(s,j,0);
v=set(v,j+1,0);
v=set(v,ld[k][j],2);
add(f[i][j+1][c[v]],f[i][j][k]);
}
}
}
for(j=1;j<=cnt;j++)
if(!get(d[j],m+1))
add(f[i+1][1][c[d[j]<<2]],f[i][m+1][j]);
}
}
};
dp a,b;
int tot;
int e1[10];
int e2[10];
int op[1010];
int f[1010];
struct list
{
int h[210];
int v[100010];
int t[100010];
int n;
list()
{
n=0;
memset(h,0,sizeof h);;
}
void add(int x,int y)
{
n++;
v[n]=y;
t[n]=h[x];
h[x]=n;
}
};
list e;
int find(int x)
{
return f[x]==x?x:f[x]=find(f[x]);
}
int check(int x,int y)
{
if(x==53&&y==53)
int xxx=1;
int i;
for(i=1;i<=m+1;i++)
{
bool x1=get(d[x],i);
bool y1=get(d[y],i);
if(x1^y1)
return 0;
}
for(i=1;i<=m+1;i++)
f[i]=i;
int fa;
for(i=1;i<=m+1;i++)
if(get(d[x],i))
{
fa=find(i);
if(get(d[x],i)==1)
f[find(i)]=find(rd[x][i]);
if(get(d[y],i)==1)
f[find(i)]=find(rd[y][i]);
}
for(i=1;i<=m+1;i++)
if(get(d[x],i)&&find(i)!=fa)
return 0;
return 1;
}
int main()
{
freopen("d2t1.in","r",stdin);
freopen("d2t1.out","w",stdout);
scanf("%d%d%d",&n,&m,&k);
int j,l;
int i,x,y;
for(i=1;i<=k;i++)
{
scanf("%d%d",&x,&y);
a.a[x][y]=b.a[n-x+1][y]=1;
}
dfs(1,0);
for(i=1;i<=cnt;i++)
for(j=1;j<=cnt;j++)
if(check(i,j))
e.add(i,j);
a.solve();
b.solve();
scanf("%d",&q);
for(i=1;i<=q;i++)
{
scanf("%d%d",&x,&y);
ll ans=0;
for(j=1;j<=cnt;j++)
{
int s=d[j];
if(get(s,m+1))
continue;
if(!get(s,y))
continue;
int bo=1;
tot=0;
for(l=1;l<=m;l++)
if(get(s,l))
{
if((a.a[x][l]||a.a[x+1][l]))
{
bo=0;
break;
}
tot++;
e1[tot]=l;
e2[tot]=get(s,l);
}
if(!bo)
continue;
for(l=e.h[j];l;l=e.t[l])
ans=(ans+a.f[x+1][1][c[s<<2]]*b.f[n-x+1][1][c[d[e.v[l]]<<2]])%p;
}
printf("%lld\n",ans);
}
return 0;
}
【BZOJ3814】【清华集训2014】简单回路 状压DP的更多相关文章
- [清华集训2015 Day1]主旋律-[状压dp+容斥]
Description Solution f[i]表示状态i所代表的点构成的强连通图方案数. g[i]表示状态i所代表的的点形成奇数个强连通图的方案数-偶数个强连通图的方案数. g是用来容斥的. 先用 ...
- URAL 1326. Bottle Taps(简单的状压dp)
题目不太好读懂,就是先给你一个n代表要从n个物品中买东西,然后告诉你这n个东西的单位价格,在给你m个集合的情况.就是每一个结合中有x件物品.他们合起来买的价格是k.这x件物品依次是:p1--px.之后 ...
- POJ 3311 Hie with the Pie (状压DP)
dp[i][j][k] i代表此层用的状态序号 j上一层用的状态序号 k是层数&1(滚动数组) 标准流程 先预处理出所有合法数据存在status里 然后独立处理第一层 然后根据前一层的max推 ...
- 状压dp Codeforces Beta Round #8 C
http://codeforces.com/contest/8/problem/C 题目大意:给你一个坐标系,给你一个人的目前的坐标(该坐标也是垃圾桶的坐标),再给你n个垃圾的坐标,这个人要捡完所有的 ...
- HDU 5067 Harry And Dig Machine(状压DP)(TSP问题)
题目地址:pid=5067">HDU 5067 经典的TSP旅行商问题模型. 状压DP. 先分别预处理出来每两个石子堆的距离.然后将题目转化成10个城市每一个城市至少经过一次的最短时间 ...
- hdu 2809(状压dp)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=2809 思路:简单的状压dp,看代码会更明白. #include<iostream> #in ...
- BZOJ5299:[CQOI2018]解锁屏幕(状压DP)
Description 使用过Android手机的同学一定对手势解锁屏幕不陌生.Android的解锁屏幕由3x3个点组成,手指在屏幕上画一条 线将其中一些点连接起来,即可构成一个解锁图案.如下面三个例 ...
- 对状压dp的见解
看了好几篇博客,终于对一些简单的状压dp有了点了解.就像HDU1074. 有个博客:https://blog.csdn.net/bentutut/article/details/70147989 感觉 ...
- 状压dp之二之三 炮兵阵地/玉米田 By cellur925
一.简单的状压dp 玉米田 题目描述 Farmer John has purchased a lush new rectangular pasture composed of M by N (1 ≤ ...
随机推荐
- itoa()函数和atoi()函数详解
C语言提供了几个标准库函数,可以将任意类型(整型.长整型.浮点型等)的数字转换为字符串. 以下是用itoa()函数将整数转换为字符串的一个例子:# include <stdio.h># i ...
- Python_每日习题_0006_斐波那契数列
程序设计: 斐波那契数列(Fibonacci sequence),从1,1开始,后面的每一项等于前面两项之和. 图方便就递归实现,图性能就用循环. # for 循环 target = int(inpu ...
- SVN插件和Tomcat插件地址
SVN插件: http://subclipse.tigris.org/update_1.8.x Tomcat插件: http://tomcatplugin.sf.net/update 备注:如果svn ...
- 使用publisher模式控制频繁的UI输出,避免Winform界面假死
http://www.cnblogs.com/Charltsing/p/publisher.html 最近测试task并发任务的效率与线程池的区别,发现了另外一个问题.task建立任务的速度很快,输出 ...
- Django之ORM操作(聚合 分组、F Q)
Django之ORM操作(聚合 分组.F Q) 聚合 aggregate()是QuerySet的一个终止子句,也就是说,他返回一个包含一些键值对的字典,在它的后面不可以再进行点(.)操作. 键的名 ...
- ocrosoft 1015 习题1.22 求一元二次方程a*x^2 + b*x + c = 0的根
http://acm.ocrosoft.com/problem.php?id=1015 题目描述 求一元二次方程a*x2 + b*x + c = 0的根.系数a.b.c为浮点数,其值在运行时由键盘输入 ...
- mysql-SQL Error: 1205, SQLState: 41000
mysql-SQL Error: 1205, SQLState: 41000——CSDN问答频道https://ask.csdn.net/questions/176492 mysql-SQL Erro ...
- MySQL慢查询日志配置方式 slow_query_log
MySQL慢查询(一) - 开启慢查询 - 鲁玉成 - 博客园 https://www.cnblogs.com/luyucheng/p/6265594.html mysql开启慢查询方法 - lava ...
- 三、如何设置npm镜像
一.临时使用 npm --registry https://registry.npm.taobao.org install express 二.永久使用 npm config set registry ...
- 基于vue-cli,sass,vant的移动端项目
项目架构 开始 vue init webpack 项目名称 //新建项目,cd进入新项目 npm install axios //先安装! ...