bzoj1311: 最优压缩

Description

 其中: Auv是与Aij相邻的像素(为了简化,认为(i-1,j),(i+1,j,(i,j-1),(i,j+1)为相邻元素); Wij取值0或者1,表示Aij修改后取V0或者V1. E的定义直观上的理解是,当修改了A之后,各像素上的值与原来的值相差了多少,以及相邻的像素对比程度的变化.为了图像的保真度,我们希望E的值越小越好.

Input

第一行二个整数N,M(1<=N,M<=35) 第二行二个整数V0,V1 接下来N行M列,对应矩阵元素Aij

Output

一个数E,表示最小可能的估价值

Sample Input

1 2
0 255
10 20

Sample Output

30

 

题解：

二元组建图：http://www.cnblogs.com/chenyushuo/p/5146626.html 

直接二元组建图即可，没有坑。。。。。

code：

 1 #include<cstdio>
 2 #include<iostream>
 3 #include<cmath>
 4 #include<cstring>
 5 #include<algorithm>
 6 #define maxn 1600
 7 #define maxm 20000
 8 #define inf 1061109567
 9 using namespace std;
10 char ch;
11 bool ok;
12 void read(int &x){
13     for (ok=0,ch=getchar();!isdigit(ch);ch=getchar()) if (ch=='-') ok=1;
14     for (x=0;isdigit(ch);x=x*10+ch-'0',ch=getchar());
15     if (ok) x=-x;
16 }
17 int n,m,v0,v1,a[42][42],idx,pos[42][42];
18 struct flow{
19     int s,t,tot,now[maxn],son[maxm],pre[maxm],val[maxm];
20     int dis[maxn],head,tail,list[maxn];
21     bool bo[maxn];
22     void init(){s=0,t=n*m+1,tot=1,memset(now,0,sizeof(now));}
23     void put(int a,int b,int c){pre[++tot]=now[a],now[a]=tot,son[tot]=b,val[tot]=c;}
24     void add(int a,int b,int c){put(a,b,c),put(b,a,0);}
25     bool bfs(){
26         memset(bo,0,sizeof(bo));
27         head=0,tail=1,list[1]=s,dis[s]=0,bo[s]=1;
28         while (head<tail){
29             int u=list[++head];
30             for (int p=now[u],v=son[p];p;p=pre[p],v=son[p])
31                 if (val[p]&&!bo[v]) bo[v]=1,dis[v]=dis[u]+1,list[++tail]=v;
32         }
33         return bo[t];
34     }
35     int dfs(int u,int rest){
36         if (u==t) return rest;
37         int ans=0;
38         for (int p=now[u],v=son[p];p&&rest;p=pre[p],v=son[p])
39             if (val[p]&&dis[v]==dis[u]+1){
40                 int d=dfs(v,min(rest,val[p]));
41                 val[p]-=d,val[p^1]+=d,ans+=d,rest-=d;
42             }
43         if (!ans) dis[u]=-1;
44         return ans;
45     }
46     int dinic(){
47         int ans=0;
48         while (bfs()) ans+=dfs(s,inf);
49         return ans;
50     }
51 }f;
52 const int dx[2]={0,1};
53 const int dy[2]={1,0};
54 int main(){
55     read(n),read(m),read(v0),read(v1),f.init();
56     for (int i=1;i<=n;i++) for (int j=1;j<=m;j++) read(a[i][j]),pos[i][j]=++idx;
57     for (int i=1;i<=n;i++) for (int j=1;j<=m;j++){
58         f.add(f.s,pos[i][j],abs(a[i][j]-v0)),f.add(pos[i][j],f.t,abs(a[i][j]-v1));
59         for (int k=0;k<2;k++){
60             int x=i+dx[k],y=j+dy[k];
61             if (x<=0||x>n||y<=0||y>m) continue;
62             f.add(pos[i][j],pos[x][y],abs(a[i][j]-a[x][y])),f.add(pos[x][y],pos[i][j],abs(a[i][j]-a[x][y]));
63         }
64     }
65     printf("%d\n",f.dinic());
66     return 0;
67 }

 

posted @ 2016-01-21 11:05  chenyushuo  阅读(285)  评论(0)  编辑  收藏  举报