bzoj 2965 平面图+暴力+网络流

平面图转对偶图

暴力枚举保护那些点

然后求最小割即可

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
#define maxn 4000
#define maxp 12
#define ll long long
#define inf 0x3fffffff
struct P{
int a,b;
}a1[maxn],p1[12];
int n,m,p;
int nn=1;
int e[maxn],ne[maxn],v[maxn],w[maxn];
int ver[maxn],he[maxn],nex[maxn],u[maxn];
void add(int x,int y,int z){
ne[++nn]=e[x],e[x]=nn,v[nn]=y,w[nn]=z;
}
int nn1=1;
void add1(int x,int y,int z){
nex[++nn1]=he[x],he[x]=nn1,u[nn1]=z,ver[nn1]=y;
}
ll mul(P x,P y,P z){
return (ll)(y.a-x.a)*(z.b-x.b)-(ll)(y.b-x.b)*(z.a-x.a);
}
int sn,s[maxn],cnt,in[maxn];
void dfs(int ee){
if(in[ee])return;
s[++sn]=ee;
in[ee]=cnt;
int y=v[ee],x=v[ee^1],now=0;
ll temp,t1,t2;
for(int i=e[y];i;i=ne[i])if(i!=(ee^1)){
t1=mul(a1[y],a1[x],a1[v[i]]);
t2=mul(a1[y],a1[v[i]],a1[v[now]]);
if(!now||temp>=0&&t1<=0||((t1<0&&temp<0||t1>0&&temp>0)&&t2<0)){
temp=t1,now=i;
}
}
dfs(now);
}
bool cross(P x1,P y1,P x2,P y2){
ll t1=mul(x1,y1,x2),t2=mul(x1,y1,y2),t3=mul(x2,y2,x1),t4=mul(x2,y2,y1);
if((t1<0&&t2>0||t1>0&&t2<0)&&(t3>0&&t4<0||t3<0&&t4>0))return 1;
return 0;
}
bool ins(int x){
int j=0;
for(int i=1;i<=sn;i++)if(cross(a1[v[s[i]]],a1[v[s[i]^1]],p1[x],(P){919941229,919940311}))j++;
return j&1;
}
int cc[12];
int S,T;
ll solv(){
ll sum=0;
for(int i=1;i<=sn;i++)sum+=mul((P){0,0},a1[v[s[i]^1]],a1[v[s[i]]]);
return sum;
}
int ch[maxn];
int tot;
bool bfs(){
for(int i=1;i<=S;i++)ch[i]=-1;
queue<int> jj;
jj.push(S);
ch[S]=0;
while(!jj.empty()){
int x=jj.front();
jj.pop();
for(int i=he[x];i;i=nex[i]){
if(ch[ver[i]]==-1&&u[i]){
ch[ver[i]]=ch[x]+1;
jj.push(ver[i]);
}
}
}
if(ch[T]!=-1)return 1;
return 0;
}
int zeng(int x,int mm){
if(x==T)return mm;
int r=mm;
for(int i=he[x];i&&r;i=nex[i])if(u[i]&&ch[ver[i]]==ch[x]+1){
int k=zeng(ver[i],min(u[i],r));
u[i]-=k;
u[i^1]+=k;
r-=k;
}
if(r==mm)ch[x]=-1;
return mm-r;
}
int dinic(){
int r=0,tt;
while(bfs())while(tt=zeng(S,inf))r+=tt;
return r;
}
int ans[maxn];
int main(){
scanf("%d%d%d",&p,&n,&m);
for(int i=1;i<=p;i++)scanf("%d%d",&p1[i].a,&p1[i].b);
for(int i=1;i<=n;i++)scanf("%d%d",&a1[i].a,&a1[i].b);
for(int i=1;i<=m;i++){
int a,b,c;
scanf("%d%d%d",&a,&b,&c);add(a,b,c);add(b,a,c);
}
for(int i=2;i<=nn;i++)if(!in[i]){
cnt++;
sn=0;
dfs(i);
if(solv()<0)T=cnt;
else for(int j=1;j<=p;j++)if(ins(j))cc[j]=cnt;
}
S=cnt+1;
for(int i=1;i<=p;i++)ans[i]=inf;
for(int op=1;op<(1<<p);op++){
int kk=0;
nn1=1;
memset(he,0,sizeof(he));
for(int i=1;i<=p;i++)if((1<<i-1)&op){
kk++;
add1(S,cc[i],inf);
add1(cc[i],S,0);
}
for(int i=2;i<=nn;i++)if(i&1){
add1(in[i],in[i^1],w[i]);
add1(in[i^1],in[i],w[i]);
}
ans[kk]=min(ans[kk],dinic());
}
for(int i=1;i<=p;i++){
cout<<ans[i]<<endl;
}
}