大致題意：給你一個數列，裡面的數字要麼是0，要麼是1~m中的一個數字。你需要隨機的把這些0替換成1~m中的任意一個，然後再計算著整個數列的權值。權值定義為，現在問你最後權值的期望。

由於這個v陣列是輸入的，也就是說沒有什麼性質，所以這個權值的表示式也不能用數論的知識推出什麼結論。因此我們考慮暴力的dp。我們令dp[i][j][k][l]表示當前計算到第i位，且a[i]取j，a[i-1]取k，a[i-2]取l時的期望的分子。那麼，從這個狀態出發，列舉a[i+1]的數值，就看求出著四個數字的gcd，對應轉移過去即可。如此一來，時間複雜度就是。顯然是不可以接受的。

接下來，我們考慮進行優化。我們注意到，如果a[i]取2，a[i-1]取4和6都是一樣的，也就是說，每一維我們根本沒有必要從1..m的列舉，我們考慮列舉對應的gcd即可。由此，我們重新定義dp陣列，dp[i][j][k][l]表示當前計算到第i位，且a[i]取j，gcd(a[i],a[i-1])取k，gcd(k,a[i-2])取l時的期望的分子。但是這麼表示，列舉起來還是很複雜，於是我們考慮提前預處理所有合法的狀態，然後把每一個狀態對應一個數字，如此陣列就變成一個二維的dp[i][j]表示當前處理到第i個數字，然後最近三個數字的狀態是狀態j的期望的分子。

然後轉移目標就是下一位，然後對應修改三個gcd的值，轉移的量就是四個數字的gcd與dp[i][j]的乘積。這樣總的複雜度就是O(N*M*tot+M^3) 其中tot表示合法的狀態數目，M^3位初始化的複雜度。最後的結論是，合法狀態數不超過1471個。具體實現見程式碼：

#include<bits/stdc++.h>
#define INF 0x3f3f3f3f
#define IO ios::sync_with_stdio(0);cin.tie(0);cout.tie(0)
#define mod 1000000007
#define N 110
#define LL long long


using namespace std;

int dp[N][1500],gcd[N][N],id[N][N][N];
struct node{int x,y,z;} st[1500];
int n,m,tot,v[N],a[N],inv[N];

void init()
{
    tot=0; inv[1]=1;
    for(int i=2;i<N;i++)
        inv[i]=1LL*(mod-mod/i)*inv[mod%i]%mod;
    for(int i=1;i<=100;i++)
        for(int j=1;j<=i;j++)
            gcd[i][j]=gcd[j][i]=__gcd(i,j);
    for(int i=1;i<=100;i++)
        for(int j=1;j<=i;j++)
            if (i%j==0)
            for(int k=1;k<=j;k++)
            {
                if (j%k) continue;
                st[++tot]=node{i,j,k};
                id[i][j][k]=tot;
            }
}

int main()
{
    IO; init();
    int T; cin>>T;
    while(T--)
    {
        cin>>n>>m;
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=n;i++) cin>>a[i];
        for(int i=1;i<=m;i++) cin>>v[i];
        for(int i=1;i<=m;i++)
            if (!a[3]||a[3]==i)
            for(int j=1;j<=m;j++)
                if (!a[2]||a[2]==j)
                for(int k=1;k<=m;k++)
                {
                    if (a[1]&&a[1]!=k) continue;
                    int y=gcd[i][j],z=gcd[y][k];
                    int nxt=id[i][y][z];
                    dp[3][nxt]++;
                }
        for(int j=3;j<n;j++)
        {
            for(int i=1;i<=tot;i++)
            {
                if (!dp[j][i]) continue;
                int x=st[i].x,y=st[i].y,z=st[i].z;
                for(int k=1;k<=m;k++)
                {
                    if (a[j+1]&&k!=a[j+1]) continue;
                    int yy=gcd[k][x],zz=gcd[k][y];
                    int nxt=id[k][yy][zz];
                    dp[j+1][nxt]+=1LL*v[gcd[k][z]]*dp[j][i]%mod;
                    if (dp[j+1][nxt]>=mod) dp[j+1][nxt]-=mod;
                }
            }
        }
        LL ans=0;
        for(int i=1;i<=tot;i++)
        {
            ans+=dp[n][i];
            if (ans>=mod) ans-=mod;
        }
        for(int i=1;i<=n;i++)
            if (!a[i]) ans=ans*inv[m]%mod;
        cout<<ans<<endl;
    }
    return 0;
}