bzoj 2669 状压DP
因为最多有8个'X',所以我们可以用w[i][s]来表示现在我们填了前i个数,填的X的为S,因为每次新加进来的数都不影响前面的最小值,所以我们可以随便添加,这样就有了剩下所有位置的方案,每次都这样转移。
但是这样会造成不是规定的地方出现局部最小值的情况,对于这样的情况,我们只需要枚举所有可能成为局部最小值的不合法状态来做容斥就可以了。
反思:这道题的容斥开始写错了,本来应该是判奇偶来判断正负,写成了全是负的,还是A了,应该是后面的容斥没有合法方案所以符号无所谓的关系,真是rp++。
/**************************************************************
Problem: 2669
User: BLADEVIL
Language: C++
Result: Accepted
Time:0 ms
Memory:928 kb
****************************************************************/
//By BLADEVIL
#include <cstdio>
#include <cstring>
#define d39 12345678
#define get(x) ((x&1)?-1:1)
using namespace std;
char s[][];
int map[][],flag[][],w[][],X[],Y[],tot[],tmp[][];
int ans,sum,n,m;
const int go[][]={{-,-},{-,},{-,},{,},{,},{,},{,-},{,-}};
int solve() {
int cnt=,x,y,flag;
for (int i=;i<=n;i++)
for (int j=;j<=m;j++) if (map[i][j]) X[cnt]=i,Y[cnt++]=j;
for (int p=;p<(<<cnt);p++) {
tot[p]=;
memset(tmp,,sizeof tmp);
for (int i=;i<cnt;i++) if (p&(<<i)) tmp[X[i]][Y[i]]=;
for (int i=;i<=n;i++)
for (int j=;j<=m;j++) if (!tmp[i][j]) {
flag=;
for (int k=;k<&&flag;k++) {
x=i+go[k][],y=j+go[k][];
if (x<||x>n||y<||y>m) continue;
if (tmp[x][y]) flag=;
}
if (flag) tot[p]++;
}
}
memset(w,,sizeof w);
w[][]=;
int opp;
for (int i=;i<=sum;i++) {
for (int p=;p<(<<cnt);p++) if (w[i-][p]) {
opp=;
for (int j=;j<cnt;j++) if (!(p&(<<j)))
opp|=(<<j),(w[i][p|(<<j)]+=w[i-][p])%=d39;
(w[i][p]+=w[i-][p]*(tot[opp]-(i-))%d39)%=d39;
}
}
return w[sum][(<<cnt)-];
}
void ie(int x,int y,int t){
map[x][y]=;
int i,j,ret;
for (int k=;k<;k++) {
i=x+go[k][],j=y+go[k][];
if (i<||i>n||j<||j>m) continue;
if (flag[i][j]) continue;
flag[i][j]=t;
}
ret=solve();
ans=(ans+get(t)*ret)%d39;
for (j=y+;j<=m;j++) if (flag[x][j]==) ie(x,j,t+);
for (i=x+;i<=n;i++)
for (j=;j<=m;j++) if (flag[i][j]==) ie(i,j,t+);
map[x][y]=;
for (int k=;k<;k++) {
i=x+go[k][],j=y+go[k][];
if (i<||i>n||j<||j>m) continue;
if (flag[i][j]==t) flag[i][j]=;
}
}
int main(){
scanf("%d%d",&n,&m);
int cur=,x,y;
sum=n*m;
for (int i = ;i <= n;i ++) {
scanf("%s",s[i]);
for (int j=;j<m;j++)map[i][j+]=s[i][j]=='X'?:;
}
for (int i=;i<=n&&cur;i++) {
for (int j=;j<=m&&cur;j++) if (map[i][j]) {
if (flag[i][j]==-) cur=;
flag[i][j]=-;
for (int k=;k<;k++) {
x=i+go[k][],y=j+go[k][];
if (x<||x>n||y<||y>m) continue;
flag[x][y]=-;
}
}
}
if (!flag) {
printf("0\n");
return ;
}
ans=solve();
for (int i=;i<=n;i++)
for (int j=;j<=m;j++) if (flag[i][j]==) ie(i,j,);
(ans+=d39)%=d39;
printf("%d\n",ans);
return ;
}
bzoj 2669 状压DP的更多相关文章
- bzoj 1879 状压dp
879: [Sdoi2009]Bill的挑战 Time Limit: 4 Sec Memory Limit: 64 MBSubmit: 852 Solved: 435[Submit][Status ...
- bzoj 1087 状压dp
1087: [SCOI2005]互不侵犯King Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 4130 Solved: 2390[Submit][ ...
- BZOJ 2064 - 状压DP
传送门 题目大意: 给两个数组, 数组中的两个元素可以合并成两元素之和,每个元素都可以分裂成相应的大小,问从数组1变化到数组2至少需要多少步? 题目分析: 看到数据范围\(n<=10\), 显然 ...
- BZOJ 4057 状压DP
思路: 状压一下 就完了... f[i]表示选了的集合为i 转移的时候判一判就好了.. //By SiriusRen #include <cstdio> #include <cstr ...
- BZOJ 4565 状压DP
思路: f[i][j][S]表示从i到j压成S状态 j-m是k-1的倍数 $f[i][j][S<<1]=max(f[i][j][S<<1],f[i][m-1][S]+f[m][ ...
- bzoj 1072状压DP
1072: [SCOI2007]排列perm Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 2293 Solved: 1448[Submit][St ...
- bzoj 1072 状压DP
我们用w[i][j]来表示,i是一个二进制表示我们选取了s中的某些位,j表示这些位%d为j,w[i][j]则表示这样情况下的方案数,那么我们可以得到转移.w[i|(1<<k)][(j*10 ...
- bzoj 1076 状压DP
我们设w[i][s]为当前到第i关,手中的物品为s的时候,期望得分为多少,其中s为二进制表示每种物品是否存在. 那么就比较容易转移了w[i][s]=(w[i-1][s']+v[j]) *(1/k),其 ...
- BZOJ 1231 状压DP
思路: f[i][j] i表示集合的组成 j表示选最后一个数 f[i][j]表示能选的方案数 f[i|(1<< k)][k]+=f[i][j]; k不属于i j属于i且符合题意 最后Σf[ ...
随机推荐
- C/S结构 B/S结构
[1]C/S 结构,即大家熟知的客户机和服务器结构.它是软件系统体系结构,通过它可以充分利用两端硬件环境的优势,将任务合理分配到Client端和Server端来实现,降低了系统的通讯开销.目前大多数应 ...
- java 基础--继承--007
1,子类只能继承父类所有非私有成员 2,子类不能继承父类的构造方法,但可以通过super去访问父类构造方法 3,子类成员变量和父类成员变量名称不一样,如果一样类似于重写,按子类处理,如果一样,就近原则 ...
- bpf移植到3.10
bpf_common.h中显示的是/usr/src/linux-headersXXXX/include/uapi/linux 竟然会识别系统的挂载选项:
- HttpServletRequestWrapper 是HttpServletRequest的包装类 ·关系相当于 int 与integer的关系
HttpServletRequestWrapper 是HttpServletRequest的包装类 ·关系相当于 int 与integer的关系
- 同步锁(synchronized)使用三要素
1.代码被多个线程访问 2.代码中有共享的数据 3.共享数据被多个语句操作
- Android开发中如何解析Json
解析Json 有了请求,自然会有响应,那我们该如何解析服务端响应返回给我们的Json数据呢? 了解什么是Json JSON(JavaScript object notation)是一种轻量级的数据交换 ...
- poj2018——Best Cow Fences
Description Farmer John's farm consists of a long row of N (1 <= N <= 100,000)fields. Each fie ...
- 斜率优化第一题! HDU3507 | 单调队列优化DP
放一手原题 题解: 第一次写(抄)斜率优化,心里还是有点小激动的.讲一下怎么实现的! 首先我们可以考虑一个朴素的dp:DP[i]表示前i个数字的最少花费,显然我们有一个转移方程 DP[i]=min{D ...
- 【树形DP】【P1351】 【NOIP2014D1T2】联合权值
传送门 Description 无向连通图 \(G\) 有 \(n\) 个点, \(n-1\) 条边.点从 \(1\) 到 \(n\) 依次编号,编号为 \(i\) 的点的权值为 \(W_i\) ,每 ...
- spring全局异常处理 自定义返回数据结构
在写api接口中,正常返回和异常错误返回我们都希望很清楚的将这些信息清楚的返回给用户,出现异常情况下需要清楚的知道是参数异常还是未知异常,而不是返回一个不正确的数据结构. 所以此处只针对写api接口时 ...