标题效果
有着n巫妖。m精灵。k木。他们都有自己的位置坐标表示。冷却时间，树有覆盖范围。
假设某个巫妖攻击精灵的路线（他俩之间的连线）经过树的覆盖范围，表示精灵被树挡住巫妖攻击不到。求巫妖杀死所有精灵的时间。若无法所有杀死输出-1；
解题思路：
推断巫妖能否打到精灵用线段与点的最短距离来推断，若最短距离小于树的覆盖范围，就攻击不到。
最小时间能够跑费用流来解决，也能够二分图的最优匹配。
以下是代码：
#include <set> #include <map> #include <queue> #include <math.h> #include <vector> #include <string> #include <stdio.h> #include <string.h> #include <stdlib.h> #include <iostream> #include <cctype> #include <algorithm> #define eps 1e-6 #define pi acos(-1.0) #define inf 107374182 #define inf64 1152921504606846976 #define lc l,m,tr<<1 #define rc m + 1,r,tr<<1|1 #define zero(a) fabs(a)<eps #define iabs(x) ((x) > 0 ? (x) : -(x)) #define clear1(A, X, SIZE) memset(A, X, sizeof(A[0]) * (min(SIZE,sizeof(A)))) #define clearall(A, X) memset(A, X, sizeof(A)) #define memcopy1(A , X, SIZE) memcpy(A , X ,sizeof(X[0])*(SIZE)) #define memcopyall(A, X) memcpy(A , X ,sizeof(X)) #define max( x, y ) ( ((x) > (y)) ?(x) : (y) ) #define min( x, y ) ( ((x) < (y)) ? (x) : (y) ) using namespace std; struct Point { double x,y; } hitted[205]; struct node1 { double r; int t; Point p; } att[205]; struct node3 { double r; Point p; } tree[205]; int n,m,K; double Distance(Point a, Point b) { return sqrt((a.x-b.x)*(a.x-b.x)*1.0 + 1.0*(a.y-b.y)*(a.y-b.y)); } double xmult(Point p1, Point p2, Point p) { return (p1.x-p.x)*(p2.y-p.y) - (p2.x-p.x)*(p1.y-p.y); } double disptoseg(Point p, Point a ,Point b) { Point t = p; t.x += a.y-b.y, t.y += b.x-a.x; if(xmult(a,t,p)*xmult(b,t,p)>eps) return Distance(p,a) < Distance(p,b)? Distance(p,a) : Distance(p,b); return fabs(xmult(p, a, b))/Distance(a, b); } int head[450],cnt; struct node { int u,v,w,f,next; } edge[230*230*20]; void addedge(int u,int v,int w,int f) { edge[cnt].u=u; edge[cnt].v=v; edge[cnt].w=w; edge[cnt].f=f; edge[cnt].next=head[u]; head[u]=cnt++; edge[cnt].u=v; edge[cnt].v=u; edge[cnt].w=-w; edge[cnt].f=0; edge[cnt].next=head[v]; head[v]=cnt++; } int pre[450]; bool vis[450]; bool spfa() { int i; //clearall(pre,-1); pre[n+m+1]=-1; pre[n+m]=-1; int dis[450]; clearall(dis,0x3f3f); queue <int >q; dis[n+m]=0; vis[n+m]=true; q.push(n+m); while(!q.empty()) { int t=q.front(); q.pop(); i=head[t]; vis[t]=false; while(i!=-1) { if(edge[i].f>0&&dis[edge[i].v]>dis[t]+edge[i].w) { dis[edge[i].v]=dis[t]+edge[i].w; pre[edge[i].v]=i; if(!vis[edge[i].v]) { vis[edge[i].v]=true; q.push(edge[i].v); } } i=edge[i].next; } } if(pre[n+m+1]==-1)return false; return true; } bool cnthit[205]; double input() { char s=getchar(); double ans=0; while((s<'0'||s>'9')&&s!='-') { s=getchar(); } bool flat=false; if(s=='-') { flat=true; s=getchar(); } while(s>='0'&&s<='9') { ans*=10; ans+=s-'0'; s=getchar(); } if(flat)ans=-ans; return ans; } int main() { int T; scanf("%d",&T); while(T--) { scanf("%d%d%d",&n,&m,&K); for(int i=0; i<n; i++) { att[i].p.x=input(); att[i].p.y=input(); att[i].r=input(); att[i].t=input(); } for(int i=0; i<m; i++) { hitted[i].x=input(); hitted[i].y=input(); } for(int i=0; i<K; i++) { tree[i].p.x=input(); tree[i].p.y=input(); tree[i].r=input(); } clearall(head,-1); clearall(cnthit,0); cnt=0; for(int i=0; i<n; i++) { int cnt1=0; for(int j=0; j<m; j++) { bool flat=true; if(Distance(att[i].p, hitted[j]) <= att[i].r) { for(int k=0; k<K; k++) { if(disptoseg(tree[k].p,att[i].p,hitted[j])<=tree[k].r) { flat=false; break; } } if(flat) { cnthit[j]=true; addedge(i,n+j,0,1); addedge(n+m,i,att[i].t*cnt1,1); cnt1++; } } } } for(K=0; K<m; K++) { if(cnthit[K]==false) { break; } } if(K!=m) { puts("-1"); continue; } for(int i=0; i<m; i++) { addedge(i+n,n+m+1,0,1); } int ans=0; while(spfa()) { int max1=inf; int p=pre[n+m+1]; while(p!=-1) { max1=min(max1,edge[p].f); p=pre[edge[p].u]; } p=pre[n+m+1]; while(p!=-1) { edge[p].f-=max1; edge[p^1].f+=max1; p=pre[edge[p].u]; } } for(int i=head[n+m];i!=-1;i=edge[i].next) { if(edge[i].f==0)ans=max(ans,edge[i].w); } printf("%d ",ans); } return 0; }
版权声明：本文博主原创文章，博客，未经同意不得转载。