点击这里查看原题
首先暴力判断哪些巫妖可以打到哪些小精灵,先判距离是否小于半径,然后判断树是否阻挡。判断树的方法是先计算出巫妖和小精灵连线的一般式ax+by+c=0,然后据此计算树的圆心到两点连线的距离
接下来就是二分时间time,源点与所有巫妖之间连容量为time/t+1的边,巫妖和能打到的精灵之间连容量为1的边,精灵和汇点连容量为1的边。然后跑最大流看最大流是不是m
/*
User:Small
Language:C++
Problem No.:1822
*/
#include<bits/stdc++.h>
#define ll long long
#define inf 999999999
using namespace std;
int n,m,t,tot,fir[405],dep[405];
bool can[205],ok[205][205];
struct edge{
int v,w,nex;
}e[100005];
struct witch{
int x,y,r,t;
}wy[205];
struct elf{
int x,y;
}jl[205];
struct tree{
int x,y,r;
}sh[205];
int sqr(int x){
return x*x;
}
void add(int u,int v,int w){
e[++tot]=(edge){v,w,fir[u]};
fir[u]=tot;
e[++tot]=(edge){u,0,fir[v]};
fir[v]=tot;
}
bool bfs(int s,int t){
memset(dep,-1,sizeof(dep));
dep[s]=0;
queue<int> q;
q.push(s);
while(!q.empty()){
int u=q.front();
q.pop();
for(int i=fir[u];i;i=e[i].nex){
int v=e[i].v;
if(e[i].w&&dep[v]==-1){
dep[v]=dep[u]+1;
q.push(v);
}
}
}
return dep[t]!=-1;
}
int dfs(int u,int f,int t){
if(u==t) return f;
int sum=0;
for(int i=fir[u];i;i=e[i].nex){
int v=e[i].v;
if(e[i].w&&dep[v]==dep[u]+1){
int d=dfs(v,min(f,e[i].w),t);
sum+=d;
f-=d;
e[i].w-=d;
e[i^1].w+=d;
if(!f) return sum;
}
}
if(sum==0) dep[u]=-1;
return sum;
}
int dinic(int s,int t){
int res=0;
while(bfs(s,t)){
res+=dfs(s,inf,t);
}
return res;
}
bool check(int t){
tot=1;
memset(fir,0,sizeof(fir));
for(int i=1;i<=n;i++){
add(1,i+2,t/wy[i].t+1);
for(int j=1;j<=m;j++) if(ok[i][j]) add(2+i,2+n+j,1);
}
for(int i=1;i<=m;i++) add(2+n+i,2,1);
return dinic(1,2)==m;
}
int main(){
freopen("data.in","r",stdin);//
scanf("%d%d%d",&n,&m,&t);
for(int i=1;i<=n;i++) scanf("%d%d%d%d",&wy[i].x,&wy[i].y,&wy[i].r,&wy[i].t);
for(int i=1;i<=m;i++) scanf("%d%d",&jl[i].x,&jl[i].y);
for(int i=1;i<=t;i++) scanf("%d%d%d",&sh[i].x,&sh[i].y,&sh[i].r);
for(int i=1;i<=n;i++){
int ax=wy[i].x,ay=wy[i].y;
for(int j=1;j<=m;j++){
int bx=jl[j].x,by=jl[j].y;
if(sqr(ax-bx)+sqr(ay-by)>sqr(wy[i].r)) continue;
int a=by-ay,b=ax-bx,c=ay*bx-by*ax;//计算巫妖和小精灵的连线的一般式ax+by+c=0
double d=sqrt(sqr(a)+sqr(b));
bool flag=0;
for(int k=1;k<=t;k++){
double dis=fabs(a*sh[k].x+b*sh[k].y+c)/d;
if(dis<sh[k].r){
flag=1;
break;
}
}
if(!flag){
can[j]=1;
ok[i][j]=1;
}
}
}
for(int i=1;i<=m;i++)
if(!can[i]){
printf("-1\n");
return 0;
}
int l=0,r=4e6,mid;
while(l<=r){
mid=l+r>>1;
if(check(mid)) r=mid-1;
else l=mid+1;
}
printf("%d\n",r+1);
return 0;
}