[NOI.AC省選模擬賽3.31] 附耳而至 [平面圖+最小割]

題面

傳送門

思路

其實就是很明顯的平面圖模型。

不咕咕咕的平面圖學習筆記

用最左轉線求出對偶圖的點，以及原圖中每一個邊兩側的點是誰

創建網絡流圖：

源點鏈接至每個對偶圖點，權值爲這個區域的光明能量

每個對偶圖點鏈接至匯點，權值爲這個區域的黑暗能量

對於每一條原圖中的邊，在它兩側的對偶圖點之間連一條雙向邊，權值爲這個邊的代價

用全部點的光明能量和黑暗能量之和，減去最小割，獲得的就是答案

Code

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cassert>
#include<cmath>
#include<map>
#include<vector>
#include<queue>
#define next DEEP_DARK_FANTASY
#define mp make_pair
#define ll long long
using namespace std;
inline int read(){
    int re=0,flag=1;char ch=getchar();
    while(!isdigit(ch)){
        if(ch=='-') flag=-1;
        ch=getchar();
    }
    while(isdigit(ch)) re=(re<<1)+(re<<3)+ch-'0',ch=getchar();
    return re*flag;
}
int n,m,T,x[100010],y[100010],light[100010],dark[100010],cnt;
vector<pair<int,int> >e[100010];
vector<pair<long double,int> >ee;
map<int,int>suf[100010],col[100010];
int suml[100010],sumd[100010];
namespace g{//當前弧dinic
    int first[100010],cnte=-1;
    struct edge{
        int to,next,w;
    }a[1000010];
    inline void init(){memset(first,-1,sizeof(first));}
    inline void add(int u,int v,int w1,int w2){
        a[++cnte]=(edge){v,first[u],w1};first[u]=cnte;
        a[++cnte]=(edge){u,first[v],w2};first[v]=cnte;
    }
    int dep[100010],cur[100010];queue<int>q;
    inline bool bfs(int s,int t){
        int i,u,v;
        for(i=s;i<=t;i++) dep[i]=-1,cur[i]=first[i];
        dep[s]=0;q.push(s);
        while(!q.empty()){
            u=q.front();q.pop();
            for(i=first[u];~i;i=a[i].next){
                v=a[i].to;if(~dep[v]||!a[i].w) continue;
                dep[v]=dep[u]+1;q.push(v);
            }
        }
        return ~dep[t];
    }
    inline int dfs(int u,int t,int lim){
        if(u==t||!lim) return lim;
        int i,v,f,flow=0;
        for(i=cur[u];~i;i=a[i].next){
            v=a[i].to;cur[u]=i;
            if(dep[v]==dep[u]+1&&(f=dfs(v,t,min(a[i].w,lim)))){
                a[i].w-=f;a[i^1].w+=f;
                flow+=f;lim-=f;
                if(!lim) return flow;
            }
        }
        return flow;
    }
    inline int dinic(int s,int t){
        int re=0;
        while(bfs(s,t)) re+=dfs(s,t,1e9);
        return re;
    }
}
int main(){
    T=read();n=read();m=read();int i,j,t1,t2,t3,tot,pos,from,next;
    for(i=1;i<=n;i++){
        x[i]=read();y[i]=read();
        light[i]=read();dark[i]=read();
    }
    for(i=1;i<=m;i++){
        t1=read();t2=read();t3=read();
        e[t1].push_back(mp(t2,t3));
        e[t2].push_back(mp(t1,t3));
    }
    for(i=1;i<=n;i++){
        tot=e[i].size();ee.clear();
        for(j=0;j<tot;j++){
            ee.push_back(mp(atan2(y[e[i][j].first]-y[i],x[e[i][j].first]-x[i]),e[i][j].first));
        }
        sort(ee.begin(),ee.end());
        for(j=0;j<tot-1;j++){//最左轉線預處理：標記每個點的後繼
            suf[i][ee[j].second]=ee[j+1].second;
        }
        if(tot) suf[i][ee[tot-1].second]=ee[0].second;
    }
    for(i=1;i<=n;i++){
        for(j=0;j<e[i].size();j++){
            pos=e[i][j].first;from=i;
            if(col[i][pos]) continue;
            cnt++;
            col[i][pos]=cnt;
            suml[cnt]+=light[pos];
            sumd[cnt]+=dark[pos];
            while(1){//求出一個區域
                next=suf[pos][from];
                if(col[pos][next]) break;
                from=pos;pos=next;
                col[from][pos]=cnt;
                suml[cnt]+=light[pos];
                sumd[cnt]+=dark[pos];
            }
        }
    }
    g::init();int ans=0;
    for(i=1;i<=cnt;i++){
        g::add(0,i,suml[i],0);
        g::add(i,n<<1,sumd[i],0);
        ans+=suml[i];ans+=sumd[i];
    }
    for(i=1;i<=n;i++){
        for(j=0;j<e[i].size();j++){
            t1=col[i][e[i][j].first];
            t2=col[e[i][j].first][i];
            if(i<e[i][j].first) g::add(t1,t2,e[i][j].second,e[i][j].second);
        }
    }
    cout<<ans-g::dinic(0,n<<1)<<'\n';
}