zoukankan      html  css  js  c++  java
  • BZOJ 2127 happiness

    传送门

    题解:和BZOJ3894一个套路,差不多的题意,给你3894的题解传送门

    #include <iostream>
    #include <cstdio>
    #include <cmath>
    #include <algorithm>
    #include <map>
    #include <queue>
    #include <vector>
    #include <cstring>
    #include <iomanip>
    #include <set>
    #include<ctime>
    //CLOCKS_PER_SEC
    #define se second
    #define fi first
    #define ll long long
    #define lson l,m,rt<<1
    #define rson m+1,r,rt<<1|1
    #define Pii pair<int,int>
    #define Pli pair<ll,int>
    #define ull unsigned long long
    #define pb push_back
    #define fio ios::sync_with_stdio(false);cin.tie(0)
    const double Pi=3.14159265;
    const int N=3e6+5;
    const ull base=163;
    const int INF=0x3f3f3f3f;
    using namespace std;
    int head[N],nx[10000010],tot=0,to[1000010];
    ll cap[1000010];
    int cur[N];
    void add(int u,int v,ll c){
        nx[tot]=head[u];
        to[tot]=v;
        cap[tot]=c;
        head[u]=tot++;
         
        nx[tot]=head[v];
        to[tot]=u;
        cap[tot]=0;
        head[v]=tot++;
    }
    int s,t,d[N];
    bool bfs(){
        memset(d,-1,sizeof(d));
        queue<int>q;
        d[s]=1;
        q.push(s);
        while(!q.empty()){
            int u=q.front();q.pop();
            for(int i=head[u];~i;i=nx[i]){
                int v=to[i];
                if(d[v]==-1&&cap[i]>0){
                    d[v]=d[u]+1;
                    q.push(v);
                }
            }
        }
        return d[t]!=-1;
    }
    ll dfs(int s,ll a){
        if(s==t||a==0)return a;
        ll flow=0,f;
        for(int &i=cur[s];~i;i=nx[i]){
            int v=to[i];
            if(d[v]==d[s]+1&&cap[i]>0&&(f=dfs(v,min(a,cap[i])))>0){
                flow+=f;
                cap[i]-=f;
                cap[i^1]+=f;
                a-=f;
                if(a==0)break;
            }
        }
        return flow;
    }
    ll dinic(){
        ll ans=0;
        while(bfs()){
            for(int i=0;i<=t;i++)cur[i]=head[i];
            while(ll flow=dfs(s,INF)){
                ans+=flow;
            }
        }
        return ans;
    }
    int n,m;
    int getid(int x,int y){
        return (x-1)*m+y;
    }
    int gets1(int x,int y){
        return (x-1)*m+y+n*m;
    }
    int gets2(int x,int y){
        return (x-1)*m+y+n*m*2;
    }
    int gets3(int x,int y){
        return (x-1)*m+y+n*m*3;
    }
    int gets4(int x,int y){
        return (x-1)*m+y+n*m*4;
    }
    int main(){
        fio;cin>>n>>m;
        memset(head,-1,sizeof(head));
        s=0,t=n*m*6+10;
        ll sum=0;
        for(int i=1;i<=n;i++){
            for(int j=1;j<=m;j++){
                int x;cin>>x;
                sum+=x;
                add(s,getid(i,j),x);
            }
        }
        for(int i=1;i<=n;i++){
            for(int j=1;j<=m;j++){
                int x;cin>>x;
                sum+=x;
                add(getid(i,j),t,x);
            }
        }
        for(int i=1;i<=n-1;i++){
            for(int j=1;j<=m;j++){
                int x;cin>>x;
                sum+=x;
                add(s,gets1(i,j),x);
                add(gets1(i,j),getid(i,j),INF);
                add(gets1(i,j),getid(i+1,j),INF);
            }
        }
        for(int i=1;i<=n-1;i++){
            for(int j=1;j<=m;j++){
                int x;cin>>x;
                sum+=x;
                add(gets2(i,j),t,x);
                add(getid(i,j),gets2(i,j),INF);
                add(getid(i+1,j),gets2(i,j),INF);
            }
        }
        for(int i=1;i<=n;i++){
            for(int j=1;j<=m-1;j++){
                int x;cin>>x;
                sum+=x;
                add(s,gets3(i,j),x);
                add(gets3(i,j),getid(i,j),INF);
                add(gets3(i,j),getid(i,j+1),INF);
            }
        }
        for(int i=1;i<=n;i++){
            for(int j=1;j<=m-1;j++){
                int x;cin>>x;
                sum+=x;
                add(gets4(i,j),t,x);
                add(getid(i,j),gets4(i,j),INF);
                add(getid(i,j+1),gets4(i,j),INF);
            }
        }
        cout<<sum-dinic();
        return 0;
    }
  • 相关阅读:
    atitit.基于http json api 接口设计 最佳实践 总结o7
    atitit.表单验证 的dsl 本质跟 easyui ligerui比较
    Atitit. 提升软件开发效率and 开发质量---java 实现dsl 4gl 的本质and 精髓 O725
    atitit.设计模式(2) -----查表模式/ command 总结
    atitit.spring3 mvc url配置最佳实践
    atitit.java给属性赋值方法总结and BeanUtils 1.6.1 .copyProperty的bug
    atitit.RESTful服务的概览and框架选型
    atitit.查看预编译sql问号 本质and原理and查看原生sql语句
    atitit.atitit.hb many2one relate hibernate 多对一关联配置..
    atitit.hbnt orm db 新新增更新最佳实践o7
  • 原文地址:https://www.cnblogs.com/Mrleon/p/9136369.html
Copyright © 2011-2022 走看看