平面图求域之后增加超级根,然后建图求最小树形图即可。
#include<cstdio>
#include<cmath>
#include<set>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=3010,M=10010,inf=0x7fffffff;
int n,m,q,cnt,i,x,y,z1,z2;
struct P{
int x,y;
P(){}
P(int _x,int _y){x=_x,y=_y;}
ll operator*(const P&b){return 1LL*x*b.y-1LL*y*b.x;}
}a[N],b[M];
struct E{
int x,y,z;double o;
E(){}
E(int _x,int _y,int _z){x=_x,y=_y,z=_z,o=atan2(a[y].x-a[x].x,a[y].y-a[x].y);}
}e[M];
bool del[M];int from[M];
namespace GetArea{
struct cmp{bool operator()(int a,int b){return e[a].o<e[b].o;}};
set<int,cmp>g[N];set<int,cmp>::iterator k;int i,j,q[M],t;
void work(){
for(i=0;i<m+m;i++)if(!del[i]){
for(q[t=1]=j=i;;q[++t]=j=*k){
k=g[e[j].y].find(j^1);k++;
if(k==g[e[j].y].end())k=g[e[j].y].begin();
if(*k==i)break;
}
ll s=0;
for(j=1;j<=t;j++)s+=a[e[q[j]].x]*a[e[q[j]].y],del[q[j]]=1;
if(s<0)continue;
for(cnt++,j=1;j<=t;j++)from[q[j]]=cnt;
}
}
}
namespace DMST{
int n,size,pre[N],id[N],vis[N],in[N];
struct EDGE{
int u,v,cost;
EDGE(){}
EDGE(int a,int b,int c):u(a),v(b),cost(c){}
}E[M];
void init(int _n){n=_n,size=0;}
void add(int u,int v,int w){E[size++]=EDGE(u,v,w);}
int dmst(int root){
int u,v,cnt,ret=0,i;
while(1){
for(i=0;i<n;i++)in[i]=inf;
for(i=0;i<size;i++){
u=E[i].u,v=E[i].v;
if(E[i].cost<in[v]&&u!=v)pre[v]=u,in[v]=E[i].cost;
}
for(cnt=in[root]=i=0;i<n;i++)id[i]=vis[i]=-1;
for(i=0;i<n;i++){
ret+=in[i],v=i;
while(vis[v]!=i&&id[v]==-1&&v!=root)vis[v]=i,v=pre[v];
if(v!=root&&id[v]==-1){
for(u=pre[v];u!=v;u=pre[u])id[u]=cnt;
id[v]=cnt++;
}
}
if(!cnt)break;
for(i=0;i<n;i++)if(id[i]==-1)id[i]=cnt++;
for(i=0;v=E[i].v,i<size;i++){
E[i].u=id[E[i].u],E[i].v=id[E[i].v];
if(E[i].u!=E[i].v)E[i].cost-=in[v];
}
n=cnt,root=id[root];
}
return ret;
}
void work(int _n){
n=_n;
int i,tot=1;
for(i=0;i<m+m;i++)if(e[i].z)add(from[i],from[i^1],e[i].z),tot+=e[i].z;
for(i=0;i<n-1;i++)add(n-1,i,tot);
printf("%d",dmst(n-1)-tot);
}
}
int main(){
scanf("%d%d",&n,&m);
for(i=1;i<=n;i++){
scanf("%d%d",&x,&y);
a[i]=P(x,y);
}
for(i=0;i<m;i++){
scanf("%d%d%d%d",&x,&y,&z1,&z2);
e[i<<1]=E(x,y,z1);
e[i<<1|1]=E(y,x,z2);
}
for(i=0;i<m+m;i++)GetArea::g[e[i].x].insert(i);
GetArea::work();
DMST::work(cnt+2);
return 0;
}