比赛 |
20181005 |
评测结果 |
AAAAAAAAAAAAAAAAAAAA |
题目名称 |
被破坏的电力系统 |
最终得分 |
100 |
用户昵称 |
观、一世沧桑如画 |
运行时间 |
0.071 s |
代码语言 |
C++ |
内存使用 |
11.94 MiB |
提交时间 |
2018-10-05 09:15:23 |
显示代码纯文本
#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <cmath>
using namespace std;
#define INF 999999999.0
#define MAX_N 1010
#define DIS(a,b) sqrt((pos[a].x-pos[b].x)*(long long)(pos[a].x-pos[b].x)+(pos[a].y-pos[b].y)*(long long)(pos[a].y-pos[b].y))
struct POS
{
int x,y;
} pos[MAX_N];
struct FLAG
{
int V[MAX_N];
int cnt;
} Flag[MAX_N];
int n,w;
double m;
bool road[MAX_N][MAX_N];
double dis[MAX_N][MAX_N];
double d[MAX_N][MAX_N];
double sd[MAX_N];
bool vis[MAX_N];
int sID,tID;
void dfs(int v,int flagID)
{
if(v==0) sID=flagID;
else if(v==n-1) tID=flagID;
vis[v]=true;
Flag[flagID].V[Flag[flagID].cnt++]=v;
for(int i=0;i<n;i++)
{
if(!vis[i] && road[v][i]) dfs(i,flagID);
}
}
#define MAX_Q 10000
#define ADD(x) x=(x+1)%MAX_Q
double SPFA(int s,int n)
{
int queue[MAX_Q],front,rear;
bool inq[MAX_N];
for(int i=0;i<n;i++) sd[i]=INF;
sd[s]=0;
front=rear=0;
queue[rear]=s;
inq[s]=true;
ADD(rear);
while(front!=rear)
{
int q=queue[front];
ADD(front);
inq[q]=false;
for(int u=0;u<n;u++)
{
if(sd[u]>sd[q]+d[q][u])
{
sd[u]=sd[q]+d[q][u];
if(inq[u]) continue;
queue[rear]=u;
ADD(rear);
inq[u]=true;
}
}
}
return sd[tID];
}
int Main()
{
freopen("pwrfail.in","r",stdin);
freopen("pwrfail.out","w",stdout);
scanf("%d%d",&n,&w);
scanf("%lf",&m);
for(int i=0;i<n;i++)
{
int a,b;
scanf("%d%d",&a,&b);
pos[i].x=a;pos[i].y=b;
}
for(int i=0;i<n;i++)
{
for(int j=0;j<n;j++)
{
d[i][j]=INF;
}
d[i][i]=0.0;
}
for(int i=0;i<w;i++)
{
int a,b;
scanf("%d%d",&a,&b);
a--;b--;
road[a][b]=true;
road[b][a]=true;
}
int flagcnt=0;
for(int i=0;i<n;i++)
{
if(!vis[i]) dfs(i,flagcnt++);
}
for(int i=0;i<n;i++)
{
for(int j=i+1;j<n;j++)
{
dis[j][i]=dis[i][j]=DIS(i,j);
}
}
for(int i=0;i<flagcnt;i++)
{
int icnt=Flag[i].cnt;
d[i][i]=0;
for(int j=i+1;j<flagcnt;j++)
{
int jcnt=Flag[j].cnt;
double minx=INF;
for(int k=0;k<icnt;k++)
{
for(int f=0;f<jcnt;f++)
{
minx=min(minx,dis[Flag[i].V[k]][Flag[j].V[f]]);
}
}
if(minx<=m)
d[j][i]=d[i][j]=minx;
else
d[j][i]=d[i][j]=INF;
}
}
double ans=SPFA(sID,flagcnt);
if(ans==INF)
printf("-1\n");
else
printf("%d",(int)(ans*1000.0));
return 0;
}
int haha=Main();
int main(){;}