#include<bits/stdc++.h>
#define ll long long
using namespace std;

const int maxn = 2100000;

int n,m,S;
double ans;
double sum[maxn],se[maxn];
int d[maxn];
vector< pair<int,int> >E[maxn];

int dfn[maxn],dfi;
void dfs(const int x,const int fa)
{
	dfn[x]=++dfi;
	for(auto [y,i]:E[x]) if(y!=fa)
		dfs(y,x);
}

int main()
{
	ios_base::sync_with_stdio(false);
	cin.tie(0);
	
	int Tcase; cin>>Tcase;
	while(Tcase--)
	{
		cin>>n>>S;
		for(int i=1;i<=n;i++)
		{
			sum[i]=0;
			d[i]=0;
			vector< pair<int,int> >_;
			E[i].swap(_);
		}
		for(int i=1;i<=n*2;i++) se[i]=0;
				
		//map< pair<int,int>,int >mp;
		for(int i=1;i<n;i++)
		{
			int x,y; cin>>x>>y;
			d[x]++; d[y]++;
			E[x].emplace_back(y,i);
			E[y].emplace_back(x,i);
			//mp[make_pair(x,y)]=i;
			//mp[make_pair(y,x)]=n+i;
		}
		dfi=0; dfs(1,0);
		ans=0;
		
		cin>>m;
		vector<int>xi(m+5);
		for(int i=1;i<=m;i++) cin>>xi[i];
		
		sum[S]=1;
		for(int i=1,las=0;i<=m;i++)
		{
			int x=xi[i];
			if(x==S && ans==0)
			{
				ans=1;
			}
			if(las)
			{
				int ei;
				if(dfn[las]<dfn[x]) ei=dfn[x];
				else ei=n+dfn[x];
				double temp= (sum[las]-se[ei])/d[las];
				se[ei]=sum[las];
				ans+=temp;
				sum[x]+=temp;
			}
			las=x;
		}
		cout<<fixed<<setprecision(12)<<ans<<'\n';
	}
	
	return 0;
}