P3160:局部极小值(容斥、状压)

解析

又是一道我不会的容斥题
qwq
本题的一个关键性质:答案有解时,极小值不超过8个
所以可以对其进行状压

考虑从小到大填数
那么在极小值填完之前,它的八连通必然是不能填的

设计dpi,sdp_{i,s}dpi,s表示从小到大填了i个数,已经填完的极小值状态为s的方案数
不难作出转移

但是这样会统计一些不合法的方案!
有的非极小值可能由于周围全是非极小值,又随便填,导致成为了极小值
所以要扣去所有非极小值成为极小值的方案

方法就是dfs枚举哪些非极小值成为极小值dp再按这些非极小值的个数的奇偶性进行容斥

代码

//暴力
#include<bits/stdc++.h>
using namespace std;
const int mod=12345678;
#define ll long long
#define il inline
il ll read(){ll x=0,f=1;char c=getchar();while(!isdigit(c)){if(c=='-')f=-1;c=getchar();}while(isdigit(c)){x=x*10+c-'0';c=getchar();}return x*f;
}
int n,m;
bool jd[20][20];
int a[20][20];
int dx[9]={0,0,-1,-1,-1,0,1,1,1},dy[9]={0,-1,-1,0,1,1,1,0,-1};
int ans;
inline bool exi(int x,int y){return x>=1&&x<=n&&y>=1&&y<=m;
}
int x[50],y[50],tot,mi[50];
bool vis[12][12];
int dp[35][1050];
int num[1050];
int calc(){tot=0;for(int i=1;i<=n;i++){for(int j=1;j<=m;j++){if(jd[i][j]){++tot;x[tot]=i;y[tot]=j;}}}for(int s=0;s<mi[tot];s++){memset(vis,0,sizeof(vis));num[s]=n*m;for(int i=1;i<=tot;i++){if(s&mi[i-1]) continue;int xx=x[i],yy=y[i];if(!vis[xx][yy]){vis[xx][yy]=1;num[s]--;}for(int k=1;k<=8;k++){int nx=xx+dx[k],ny=yy+dy[k];if(exi(nx,ny)&&vis[nx][ny]==0){vis[nx][ny]=1;num[s]--;}}}//printf("s=%d num=%d\n",s,num[s]);}memset(dp,0,sizeof(dp));dp[0][0]=1;for(int i=1;i<=n*m;i++){for(int s=0;s<mi[tot];s++){dp[i][s]+=1ll*dp[i-1][s]*max(num[s]-i+1,0)%mod;if(dp[i][s]>=mod) dp[i][s]-=mod;for(int k=1;k<=tot;k++){if((s&mi[k-1])==0) continue;dp[i][s]+=dp[i-1][s-mi[k-1]];if(dp[i][s]>=mod) dp[i][s]-=mod;}}}return dp[n*m][mi[tot]-1];
}
void dfs(int x,int y,int o){if(x>n){/*for(int i=1;i<=n;i++){for(int j=1;j<=m;j++) printf("%d ",jd[i][j]);putchar('\n');}*/if(o&1){ans-=calc();if(ans<0) ans+=mod;}else{ans+=calc();if(ans>=mod) ans-=mod;}//printf("ans=%d\n\n",ans);return;}if(y>m){dfs(x+1,1,o);return;}dfs(x,y+1,o);if(!jd[x][y]){for(int i=1;i<=8;i++){int nx=x+dx[i],ny=y+dy[i];if(exi(nx,ny)&&jd[nx][ny]) return;}jd[x][y]=1;dfs(x,y+1,o+1);jd[x][y]=0;}
}
int main(){
#ifndef ONLINE_JUDGEfreopen("a.in","r",stdin);freopen("a.out","w",stdout);#endifmi[0]=1;for(int i=1;i<=28;i++) mi[i]=mi[i-1]<<1;n=read();m=read();char c;for(int i=1;i<=n;i++){for(int j=1;j<=m;j++){scanf(" %c",&c);jd[i][j]=c=='X';}}for(int i=1;i<=n;i++){for(int j=1;j<=m;j++){if(!jd[i][j]) continue;for(int k=1;k<=8;k++){int nx=i+dx[k],ny=j+dy[k];if(exi(nx,ny)&&jd[nx][ny]){printf("0");return 0;}}}}dfs(1,1,0);printf("%d\n",ans);
}
/**/

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/317679.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

[SNOI2017]遗失的答案 (FWT)

description 小皮球在计算出答案之后&#xff0c;买了一堆皮肤&#xff0c;他心里很开心&#xff0c;但是一不小心&#xff0c;就忘记自己买了哪些皮肤了。 ||| 万幸的是&#xff0c;他还记得他把所有皮肤按照 1∼N 来编号&#xff0c;他买来的那些皮肤的编号&#xff08;他至…

Abp中使用可视化的日志面板

如果你还不了解LogDashboard请看这里 使用logdashboard查看可视化日志。ABP的相关知识不做介绍如果有需要请阅读ABP官方文档ABP是Net下非常优秀的开发框架,在中国很多的项目都正在使用它。现在我们可以使用LogDashboard增强在使用ABP开发中的查看日志能力。下载ABP模板项目打开…

J - Just Multiplicative Inverse Gym - 102875J

J - Just Multiplicative Inverse Gym - 102875J 题目&#xff1a; 题解&#xff1a; 给定一个x&#xff0c;求出F(1,x)F(2,x)…F(x-1,x) 的和除以&#xff08;x-1&#xff09; F(x,p)题目已经给出 我们观察F()含义&#xff0c;再结合本题含义&#xff0c;本题并不是要求F(x,…

P7137-[THUPC2021 初赛]切切糕【dp】

正题 题目链接:https://www.luogu.com.cn/problem/P7137 题目大意 有两个人&#xff0c;有nnn个蛋糕&#xff0c;第iii个蛋糕大小为aia_iai​。 每一次第一个人可以选择一个蛋糕把它切成任意大小的两份&#xff08;一份可以为空&#xff09;。 然后第二个人有mmm次机会优先…

利用Topshelf把.NET Core Generic Host管理的应用程序部署为Windows服务

2019第一篇文章。此文源于前公司在迁移项目到.NET Core的过程中&#xff0c;希望使用Generic Host来管理定时任务程序时&#xff0c;没法部署到Windows服务的问题&#xff0c;而且官方也没给出解决方案&#xff0c;只能关注一下官方issue #809 等他们方解决了。官方文档只提供了…

洛谷P3270:成绩比较(容斥、组合数学)

解析 依然不会亚qwq 但这次至少有点上道了 (指推出了一个会导致重复计数的错误式子) 首先&#xff0c;我们要选出碾压那些人&#xff0c;方案数就是Cn−1kC_{n-1}^kCn−1k​ 然后&#xff0c;我们要统计每门学科的排名情况 考虑比B神分数高的人一定是从没被碾压的人里选。所…

游戏 (匈牙利)

description 在 2016 年&#xff0c;佳媛姐姐喜欢上了一款游戏&#xff0c;叫做泡泡堂。 简单的说&#xff0c;这个游戏就是在一张地图上放上若干个炸弹&#xff0c;看是否能炸到对手&#xff0c;或者躲开对手的炸弹。在玩游戏的过程中&#xff0c;小 H 想到了这样一个问题&a…

D - Delete Prime Gym - 102875D

D - Delete Prime Gym - 102875D 题意&#xff1a; 长度为n的序列&#xff0c;每次从中取出编号为1或质数的数&#xff0c;组成新的序列d 现在给一个k&#xff0c;问序列中哪一位的值是k以及序列中第k位是多少 共T个询问 T< 2 * 105 n&#xff0c;k< 105 题解&#…

CF1628A-Meximum Array【二分】

正题 题目链接:http://codeforces.com/contest/1628/problem/A 题目大意 给出一个长度为nnn的序列aaa和一个空序列bbb&#xff0c;你每次可以选择aaa的一个前缀&#xff0c;将它的mexmexmex加入序列bbb的末尾&#xff0c;然后将aaa的这个前缀删除。 求bbb的最大字典序。 1≤…

gym 102875 H. Happy Morse Code

gym 102875 H. Happy Morse Code 题意&#xff1a; 一个长度为n的字符串&#xff0c;现在给你m个小字符串&#xff0c;问小字符串拼成大字符串有多少种方法&#xff1f; 答案mod128 题解&#xff1a; 其实也不难&#xff0c;但是本人对于dp及其不敏感&#xff0c;我看出来是…

开源项目商业模式分析(2) - 持续维护的重要性 - Selenium和WatiN

该系列第一篇发布后收到不少反馈&#xff0c;包括&#xff1a;第一篇里说的MonicaHQ不一定盈利没错&#xff0c;但是问题在于绝大多数开源项目商业数据并没有公开&#xff0c;从而无法判断其具体是否盈利。难得MonicaHQ是公开的&#xff0c;所以才用来做这系列文章的开篇。很多…

[SCOI2015]小凸玩矩阵 (匈牙利+二分)

description 题目描述 小凸和小方是好朋友&#xff0c;小方给小凸一个 NM&#xff08;N≤M&#xff09;的矩阵 A&#xff0c;要求小凸从其中选出 N 个数&#xff0c;其中任意两个数字不能在同一行或同一列&#xff0c;现小凸想知道选出来的 N 个数中第 K大的数字的最小值是多…

洛谷P4971:断罪者(左偏树)

解析 看起来就是左偏树的基本操作啊… 然而就是调不过去 吐了qwq 参考了望月大神的实现 感觉清晰的多 就定义并查集维护的是每个点所在的堆的根节点 一下子少了很多恶心的套娃 代码 #include<bits/stdc.h> using namespace std; #define ll long long #define il inl…

CF1267G-Game Relics【数学期望,dp】

正题 题目链接:https://www.luogu.com.cn/problem/CF1267G 题目大意 给出nnn个物品&#xff0c;你可以进行如下操作 花费xxx获得随机一个物品。花费cic_ici​获得第iii个物品。 1≤n≤100,1≤x≤10000,∑ai≤1041\leq n\leq 100,1\leq x\leq 10000,\sum a_i\leq 10^41≤n≤1…

深入业务成为更好的软件架构师——信息化建设图鉴一二例

软件开发实际上跟英语比较类似&#xff0c;都是一项工具&#xff0c;服务于各行各业。从程序员的个人修养上来讲&#xff0c;一是要研习好软件开发这门技艺&#xff0c;二是要深入到所服务的行业。说到底&#xff0c;软件的终极目标是模拟业务&#xff0c;在此期间常常会有一个…

2020牛客暑期多校训练营(第三场)

比赛链接 题号题目通过率AClam and Fish1453/4902BClassical String Problem1547/6542COperation Love1431/5631DPoints Construction Problem548/2238ETwo Matchings943/4141FFraction Construction Problem1034/6445GOperating on a Graph972/4259HSort the Strings Revisio…

文理分科 (最小割问题)

Description 文理分科是一件很纠结的事情&#xff01;&#xff08;虽然看到这个题目的人肯定都没有纠结过&#xff09; 小P所在的班级要进行文理分科。他的班级可以用一个n*m的矩阵进行描述&#xff0c;每个格子代表一个同学的座位。每位同学必须从文科和理科中选择一科。同学…

NOIP2014洛谷P2296:寻找道路(bfs)

解析 杀鸡焉用AS47 做的巨差的一道题 WA3遍写的巨麻烦复杂度被吊打 qwq 说起来很玄学&#xff0c;但第一交之前就有一种强烈的预感觉得自己会WA 一开始像个傻子一样对有环图上记搜 WA了还不知道为啥… 可能真是傻了吧qwq 过了一会看着lojn50,m800干瞪了会眼后终终终于发现了bu…

P7097-[yLOI2020]牵丝戏【dp】

正题 题目链接:https://www.luogu.com.cn/problem/P7097 题目大意 有两个人在打游戏&#xff0c;开始时每个人有个初始伤害xxx和delay值ddd。 有mmm种道具&#xff0c;每个道具一个回合只能用一次&#xff0c;可以使当前回合玩家伤害初始伤害的ki10000\frac{k_i}{10000}1000…

恭贺微软技术俱乐部苏州站正式成立

今天去苏州微软中国&#xff0c;参加了微软技术俱乐部苏州站的成立大会。大会的历程悉数经历&#xff0c;这会儿仍在为各位大佬、社区领袖的奉献精神所感动。在通过我们公司同事分享的链接报名时就了解到&#xff0c;大会上有苏震巍老师和蒋金楠老师的分享&#xff0c;便笃定了…