湖南集训day6

难度:☆☆☆☆☆☆☆☆

/*
对于第一问：f[i][j]表示前i个数，当前黑板上的数为j的概率
当前有三种情况
1．当前数不是j的倍数—>黑板上的数字改变。
2．当前数是j的倍数且当前数在前i个数中(已经选过)
3．当前数是j的倍数且没有选过
转移:f[i+1][j]=((j的倍数个数-i)*f[i][j]+f[i][gcd(j,k)])  的平均值  j的倍数个数-i是没选过的j的倍数。
对于第二问，考虑博弈论中sg函数。可知sg[i][1]二维含义同f数组)必定为0(最后黑板上剩下1必败)  sg[n][i]=0(选完了必败) 同样枚举上述三种情况，取后续状态mex值即可。
*/
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define eps 1e-8
#define N 1000
using namespace std;
inline int gcd(int a, int b)
{
    return b ==  ? a : gcd(b, a % b);
}
int sg[N + ][N + ], g[N + ][N + ], f[N + ], n, a[N + ];
double dp[N + ][N + ], ans = ;
bool getsg(int x, int y)
{
    if (x == ) return ;
    if (sg[x][y] != -) return sg[x][y];
    bool flag = ;
    if (f[x] > y) flag &= getsg(x, y + );
    for (int i = ; i <= n; i++)
        if (g[x][i] != x)
            flag &= getsg(g[x][i], y + );
    sg[x][y] = !flag;
    return sg[x][y];
}
int main()
{
    freopen("cards.in", "r", stdin);
    freopen("cards.out", "w", stdout);
    scanf("%d", &n);
    int mx = ;
    for (int i = ; i <= n; i++)
    {
        scanf("%d", a + i);
        mx = max(mx, a[i]);
        g[][i] = a[i];
    }
    for (int i = ; i <= mx; i++)
        for (int j = ; j <= n; j++)
            f[i] += (a[j] % i == ), g[i][j] = gcd(i, a[j]);
    dp[][] = ;
    for (int i = ; i <= n; i++)
        for (int j = ; j <= mx; j++)
            if (dp[i - ][j] > eps)
            {
                dp[i][j] += dp[i - ][j] * (f[j] - i + ) / (n - i + );
                for (int k = ; k <= n; k++)
                    if (g[j][k] != j)
                    {
                        if (g[j][k] != )
                            dp[i][g[j][k]] += dp[i - ][j] / (n - i + );
                        else
                            ans += (i +  & ) * dp[i - ][j] / (n - i + );
                    }
            }
    if (n & )
        for (int j = ; j <= mx; j++) ans += dp[n][j];
    printf("%.9lf ", ans);
    memset(sg, -, sizeof(sg));
    if (getsg(, )) puts("1.000000000");
    else puts("0.000000000");
    return ;
}

/*
对于60~80分，可以n^2暴力，断掉每条边时，O(N)求每个部分的直径，然后相乘。

正解：倒序加边，考虑两棵树合并的时候新直径一定是原来两个直径四个断点任意两个的路径。
所以可以首先求LCA倍增处理两点间路径，然后求最大。更新答案要用逆元。
*/
#include<iostream>
#include<cstdio>
#include<cstring>

#define N 100007
#define mod 1000000007
#define ll long long

using namespace std;
int a[N],dep[N],sum[N],fa[N][],head[N],del[N];
int D[N][],ans[N],f[N],len[N],end[];
int n,m,pre,cnt;
struct edge{
    int u,v,net;
}e[N<<];

inline int read()
{
    int x=,f=;char c=getchar();
    while(c>''||c<''){if(c=='-')f=-;c=getchar();}
    while(c>=''&&c<=''){x=x*+c-'';c=getchar();}
    return x*f;
}

inline void add(int u,int v)
{
    e[++cnt].u=u;e[cnt].v=v;e[cnt].net=head[u];head[u]=cnt;
}

void dfs(int u,int from)
{
    fa[u][]=from;dep[u]=dep[from]+;
    sum[u]=sum[from]+a[u];
    for(int i=;i<=;i++) fa[u][i]=fa[fa[u][i-]][i-];
    for(int i=head[u];i;i=e[i].net)
    {
        int v=e[i].v;
        if(v!=from)dfs(v,u);
    }return;
}

int find(int x)
{
    return x==f[x]?x:f[x]=find(f[x]);
}

int ksm(int a,int b)
{
    int res=;
    while(b)
    {
        if(b&) res=1ll*res*a%mod;
        b>>=;a=1ll*a*a%mod;
    }return res%mod;
}

int lca(int u,int v)
{
    if(dep[u]<dep[v]) swap(u,v);
    int t=dep[u]-dep[v];
    for(int i=;i<;i++)
      if(t&(<<i)) u=fa[u][i];
    if(u==v) return u;
    for(int i=;i>=;i--)
    {
        if(fa[u][i]!=fa[v][i])
        {
            u=fa[u][i];
            v=fa[v][i];
        }
    }return fa[u][];
}

int getlen(int u,int v)
{
    int L=lca(u,v);
    return sum[u]+sum[v]-*sum[L]+a[L];
}

int main()
{
    freopen("forest.in","r",stdin);
    freopen("forest.out","w",stdout);
    int x,y;
    n=read();pre=;
    for(int i=;i<=n;i++)
    {
        a[i]=read();f[i]=i;
        pre=(ll)pre*a[i]%mod;
        D[i][]=D[i][]=i;len[i]=a[i];
    }
    for(int i=;i<n;i++)
    {
        x=read();y=read();
        add(x,y);add(y,x);
    }
    dfs(,);int t=n;ans[n]=pre;
    for(int i=;i<n;i++) del[i]=read();

    for(int i=n-;i;i--)
    {
        int id=del[i],u=e[id*-].u,v=e[id*-].v;
        u=find(u);v=find(v);
        if(len[u]<len[v]) swap(u,v);
        int tmax=len[u];
        for(int j=;j<;j++) end[j]=D[u][j];
        for(int j=;j<;j++)
          for(int k=;k<;k++)
            {
                int l=getlen(D[u][j],D[v][k]);
                if(l>tmax)
                {
                    tmax=l;
                    end[]=D[u][j];end[]=D[v][k];
                }
            }

        pre=(ll) pre*ksm(len[u],mod-)%mod;
        pre=(ll) pre*ksm(len[v],mod-)%mod;
        f[v]=u;len[u]=tmax;
        for(int j=;j<;j++) D[u][j]=end[j];
        pre=(ll) pre*len[u]%mod;
        ans[--t]=pre;
    }
    for(int i=;i<=n;i++) printf("%d\n",ans[i]);
    return ;
}

/*
这题确实是恶心，并且貌似没有部分分......
让我们来瞻仰一下std做法
考虑两列的情况。若两列颜色分别为A,B,则A独有的颜色就是A—A∩B ，B同理。
若是多列那还是设两边两列为A,B，中间多列为C，那根据题目结论可以知道C一定是A∩B的子集。枚举A中独有颜色个数，B中独有颜色个数与A中相同。若有j独有，i共有C(K，i)*C(k-i+1，j)*C(k-i-j，j)
因为每次选择都必须是恰好那些颜色，不能少，所以用总方案数减去不是恰好的就可以了。
*/
# include<iostream>
# include<cstdio>
# include<cstring>
# include<cstdlib>

using namespace std;
const int pp=;
int c[][],f[],p[],ni[];
int n,m,k,nn;

inline int power(int x,int n)
{
    int ans=,tmp=x;
    while (n)
    {
          if (n&) ans=(long long)ans*tmp%pp;
          tmp=(long long)tmp*tmp%pp;n>>=;
    }
    return ans;
}

void Count_c()
{
     for (int i=;i<=nn;i++) c[i][]=;
     for (int i=;i<=nn;i++)
      for (int j=;j<=i;j++)
      {
          c[i][j]=c[i-][j-]+c[i-][j];
          if (c[i][j]>=pp) c[i][j]-=pp;
      }
}

void Count_p()
{
     int mm=(m-)*n;
     for (int i=;i<=nn;i++)
      p[i]=power(i,mm);
}

void Count_f()
{
     f[]=;f[]=;
     for (int i=;i<=nn;i++)
     {
         f[i]=power(i,n);
         for (int j=;j<i;j++)
         {
             f[i]-=(long long)f[j]*c[i][j]%pp;
             if (f[i]<=-pp) f[i]+=pp;
         }
         if (f[i]<) f[i]+=pp;
     }
}

void Count_ni()
{
     ni[]=;
     for (int i=;i<=nn;i++)
     ni[i]=power(i,pp-);
}

int main()
{
    freopen("photo.in","r",stdin);
    freopen("photo.out","w",stdout);
    scanf("%d%d%d",&n,&m,&k);
    nn=min(n,k);
    if (m==)
       printf("%d\n",power(k,n));
    else
    {
        Count_c();Count_p();
        Count_f();Count_ni();
        long long tmp=,tmp1=,sum=,sum1;
        for (int s=;s<=nn;s++)
        {
            tmp=tmp*ni[s]%pp;
            tmp=tmp*(k-s+)%pp;
            tmp1=;sum1=;
            for (int j=;j<=s;j++)
            {
                sum1+=tmp1*c[s][s-j]%pp*p[s-j]%pp;
                if (sum1>=pp) sum1-=pp;
                tmp1=tmp1*ni[j+]%pp;
                if (k-s<j+) break;
                tmp1=tmp1*(k-s-j)%pp;
            }
            sum+=tmp*f[s]%pp*f[s]%pp*sum1%pp;
            if (sum>=pp) sum-=pp;
        }
        printf("%d\n",sum);
    }
    fclose(stdin);
    fclose(stdout);
    return ;
}