BZOJ 2669([cqoi2012]局部极小值-状态压缩+dp)
来源:互联网 发布:centos安装无线网卡 编辑:程序博客网 时间:2024/05/16 09:33
有一个n行m列的整数矩阵,其中1到nm之间的每个整数恰好出现一次。如果一个格子比所有相邻格子(相邻是指有公共边或公共顶点)都小,我们说这个格子是局部极小值。
给出所有局部极小值的位置,你的任务是判断有多少个可能的矩阵。(1<=n<=4, 1<=m<=7)
状压dp+容斥 容斥哪些X是此次DP的内容 dp[i][2^9] 表示该填i了 2^9表示哪些X被填了
首先dfs所有的X,看哪些X是这次DP要考虑的
#include<cstdio>#include<cstring>#include<cstdlib>#include<algorithm>#include<functional>#include<iostream>#include<cmath>#include<cctype>#include<ctime>#include<iomanip> #include<vector>#include<string>#include<queue>#include<stack>#include<map>#include<sstream>using namespace std;#define For(i,n) for(int i=1;i<=n;i++)#define Fork(i,k,n) for(int i=k;i<=n;i++)#define Rep(i,n) for(int i=0;i<n;i++)#define ForD(i,n) for(int i=n;i;i--)#define ForkD(i,k,n) for(int i=n;i>=k;i--)#define RepD(i,n) for(int i=n;i>=0;i--)#define Forp(x) for(int p=Pre[x];p;p=Next[p])#define Forpiter(x) for(int &p=iter[x];p;p=Next[p]) #define Lson (o<<1)#define Rson ((o<<1)+1)#define MEM(a) memset(a,0,sizeof(a));#define MEMI(a) memset(a,127,sizeof(a));#define MEMi(a) memset(a,128,sizeof(a));#define INF (2139062143)#define F (12345678)#define pb push_back#define mp make_pair #define fi first#define se second#define vi vector<int> #define pi pair<int,int>#define SI(a) ((a).size())#define Pr(kcase) printf("Case #%d: ",kcase);#define PRi(a,n) For(i,n-1) cout<<a[i]<<' '; cout<<a[n]<<endl;#define PRi2D(a,n,m) For(i,n) { \ For(j,m-1) cout<<a[i][j]<<' ';\ cout<<a[i][m]<<endl; \ } #pragma comment(linker, "/STACK:102400000,102400000")typedef long long ll;typedef long double ld;typedef unsigned long long ull;ll mul(ll a,ll b){return (a*b)%F;}ll add(ll a,ll b){return (a+b)%F;}ll sub(ll a,ll b){return ((a-b)%F+F)%F;}void upd(ll &a,ll b){a=(a%F+b%F)%F;}int read(){ int x=0,f=1; char ch=getchar(); while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();} while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();} return x*f;} #define MAXN (100)int n,m;int p2[100];ll f[MAXN][1<<9];int g[MAXN],g2[MAXN],dir[8][2]={{1,0},{-1,0},{0,1},{0,-1},{1,1},{1,-1},{-1,1},{-1,-1}};bool inside(int x,int y) { return 1<=x&&1<=y&&x<=n&&y<=m;}int id[MAXN][MAXN];int idx(int i,int j) { return (i-1)*m+j-1;}int x[MAXN],y[MAXN],siz=0;int b[50];int calc() { if (siz>9) return 0; memset(b,-1,sizeof(b)); For(i,n) For(j,m) id[i][j]=idx(i,j); MEM(g) MEM(g2) Rep(i,siz) { int X=x[i],Y=y[i]; g2[i]=1<<idx(X,Y); b[id[X][Y]]=i; Rep(di,8) { int a=X+dir[di][0],b=Y+dir[di][1]; if (inside(a,b)) { int p=1<<idx(a,b); g[i]|=p; } } } MEM(f) f[0][0]=1; int S=1<<siz; int all=(1<<(n*m))-1; Rep(i,n*m) { Rep(st,S) if (f[i][st]) { int now=all; Rep(j,siz) if (st&p2[j]) { now-=now&g2[j]; } else { now-=now&g[j]; } int tot=0; int tot1=0,tot2=0; Rep(p,n*m) if (now&p2[p]) { tot++; if (-1 == b[p]) tot2++; else tot1++; } int lef=n*m-i; int choseX=siz-tot1; int useBlank=i-choseX; tot2-=useBlank; tot2=min(tot2,lef-tot1); if (tot2>0) upd(f[i+1][st],mul(tot2,f[i][st])); Rep(p,n*m) { if (now&p2[p]) { if (b[p]==-1) { } else { upd(f[i+1][st|p2[b[p]]],f[i][st]); } } } } } return f[n*m][S-1];} bool ma[10][10];ll ans=0;void dfs(int dep,int l) { if (dep==1) upd(ans,calc()); else ans=sub(ans,calc()); For(i,n) For(j,m) if (idx(i,j)>=l && !ma[i][j]){ bool fl=0; Rep(di,8) { int a=i+dir[di][0],b=j+dir[di][1]; if (!inside(a,b)) continue; if (ma[a][b]) { fl=1; break; } } if (!fl) { ma[i][j]=1; x[siz]=i,y[siz]=j; ++siz; dfs(-dep,idx(i,j)+1); ma[i][j]=0; --siz; } }}int main(){ p2[0]=1; For(i,30) p2[i]=p2[i-1]*2; int kcase=1; while(cin>>n>>m) { siz=0; For(i,n) { char s[10]; scanf("%s",s); Rep(j,m) if (s[j]=='X') { x[siz]=i; y[siz]=j+1; ++siz; ma[i][j+1]=1; } else ma[i][j+1]=0; } ans=0; dfs(1,0); cout<<ans<<endl; } return 0;}
0 0
- BZOJ 2669([cqoi2012]局部极小值-状态压缩+dp)
- [DP 容斥] BZOJ 2669 [cqoi2012]局部极小值
- bzoj 2669: [cqoi2012]局部极小值
- BZOJ 2669 局部极小值 CQOI2012
- BZOJ 2669- [cqoi2012]局部极小值
- BZOJ 2669: [cqoi2012]局部极小值
- bzoj 2669: [cqoi2012]局部极小值
- BZOJ 2669 cqoi2012 局部极小值 状压DP+容斥原理
- BZOJ 2669 cqoi2012 局部极小值 状压DP+容斥原理
- bzoj 2669: [cqoi2012]局部极小值(dp+容斥原理)
- 【bzoj 2669】[cqoi2012]局部极小值(状压dp+容斥原理)
- BZOJ 2669: [cqoi2012]局部极小值 状压dp 容斥原理
- 【bzoj2669】【cqoi2012】【局部极小值】【状压dp】
- 【CQOI2012】局部极小值
- bzoj2669: [cqoi2012]局部极小值
- BZOJ2669: [cqoi2012]局部极小值
- CQOI2012局部极小值 BZOJ2669 状压DP+容斥
- [BZOJ2669][cqoi2012][状压DP][容斥原理]局部极小值
- session
- 批处理操作系统
- POJ3258【二分】
- html+css3
- Linux Netlink 基本使用
- BZOJ 2669([cqoi2012]局部极小值-状态压缩+dp)
- css3动画
- 关于运算符号 ++ -- 的一个最简单明了的说明
- 页面布局
- jquery+php实现用户输入搜索内容时自动提示
- NFC相关
- 响应式布局
- 1002 1002 A + B Problem II大数相加
- 求数组中第k大的数的位置