银河历59451年，在银河系有许许多多已被人类殖民的星系。如果想要在行

第一行是N，M，分别表示这个银河系内的行星系的个数

总共Q行，每一行一个实数，表示花花第i次需要计算的答

6 6
3 4
5 6
1 3
4 4
3 3
2 4
1 2
1 3
2 3
2 4
3 5
5 6
3
3 6
2 4
4 6

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define RG register
using namespace std;
const int N=300050;
int n,m,cnt,tot,root[N],head[N],to[N],nxt[N],fa[N],dep[N],size[N],son[N],top[N],cir[N],rt[N];
int tx[N],ty[N],vis[N];
int gi()
{int x=0,flag=1;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-') flag=-1;ch=getchar();}while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();return x*flag;
}
struct data{  int pre[6];  void update(RG int u,RG int v){  pre[0]++;pre[1]+=u;pre[2]+=u*u;  pre[3]+=v;pre[4]+=v*v;pre[5]+=u*v;  }  double query(){RG double n=pre[0],totx=pre[1],toty=pre[3],totx2=pre[2],toty2=pre[4],totxy=pre[5];RG double xar=totx/n,yar=toty/n;RG double a=totx2-2*xar*totx+n*xar*xar;RG double b=2*yar*totx+2*xar*toty-2*totxy-2*n*xar*yar;RG double c=toty2-2*yar*toty+n*yar*yar;RG double A=4.0,B=-4*(a+c),C=4*a*c-b*b;return (-B-sqrt(B*B-4*A*C))/2/A;}
}a[N],b[N];
data operator +(data x,data y){  for(RG int i=0;i<6;i++) x.pre[i]+=y.pre[i]; return x;
}
data operator -(data x,data y){  for(RG int i=0;i<6;i++) x.pre[i]-=y.pre[i]; return x;
}
void lnk(RG int x,RG int y){  to[++tot]=y;nxt[tot]=head[x];head[x]=tot;
}
inline void dfs1(int x,int f){  vis[x]=1;size[x]=1;son[x]=0;rt[x]=f;  a[x]=a[fa[x]]; a[x].update(tx[x],ty[x]);  for(RG int i=head[x];i;i=nxt[i]){  RG int y=to[i];  if (!vis[y]&&y!=fa[x]){  dep[y]=dep[x]+1;fa[y]=x;  dfs1(y,f);size[x]+=size[y];  if(size[y]>size[son[x]]) son[x]=y;  }  }
}
inline void dfs2(RG int x,RG int f){  top[x]=f;  if(son[x]) dfs2(son[x],f);  for (RG int i=head[x]; i; i=nxt[i]){  RG int y=to[i];  if(x==fa[y] && y!=son[x]) dfs2(y,y);  }
}
inline int lca(RG int x,RG int y){while(top[x]!=top[y]){if(dep[top[x]]<dep[top[y]]) swap(x,y);x=fa[top[x]];}if(dep[x]<dep[y]) swap(x,y);return y;
}
inline void work(){  memset(vis,0,sizeof(vis));  for(RG int x=1;x<=n;x++)  for(RG int i=head[x];i;i=nxt[i]){if(fa[x]!=to[i]&&fa[to[i]]!=x){  RG int y=to[i];  if(dep[x]>dep[y]) swap(x,y);  for(;y!=x;y=fa[y]){  root[++cnt]=y;cir[y]=cnt;vis[y]=1;  b[cnt]=b[cnt-1];b[cnt].update(tx[y],ty[y]);  }  root[++cnt]=x;cir[x]=cnt;vis[x]=1;  b[cnt]=b[cnt-1];b[cnt].update(tx[x],ty[x]);  return;  }}
}
int main(){n=gi(),m=gi();RG int x,y;  for(RG int i=1;i<=n;i++) tx[i]=gi(),ty[i]=gi();  for(RG int i=1;i<=m;i++){x=gi(),y=gi();lnk(x,y); lnk(y,x);  }  dfs1(1,1);  memset(vis,0,sizeof(vis));  if(n==m) work();else root[cnt=1]=1;  memset(fa,0,sizeof(fa));   for(RG int i=1;i<=cnt;i++){  dfs1(root[i],root[i]); dfs2(root[i],root[i]);  }  RG int Q=gi();data u,v;  for(RG int i=1;i<=Q;i++){  x=gi(),y=gi();  if (rt[x]==rt[y]){  u=a[x]+a[y]-a[lca(x,y)]-a[fa[lca(x,y)]];  printf("%.5f\n",u.query());  } else{  if (cir[rt[x]]>cir[rt[y]]) swap(x,y);  u=a[x]-a[rt[x]]+a[y]-a[rt[y]]+b[cir[rt[y]]]-b[cir[rt[x]]-1];  v=a[x]-a[rt[x]]+a[y]-a[rt[y]]+b[cir[rt[x]]]+b[cnt]-b[cir[rt[y]]-1];  printf("%.5f\n",min(u.query(),v.query()));  }  }  return 0;
}