比较实验

比较实验

时间限制: 1 Sec 内存限制: 512 MB

题目描述
小D被邀请到实验室，做一个跟图片质量评价相关的主观实验。实验用到的图片集一共有N张图片，编号为1到N。实验分若干轮进行，在每轮实验中，小D会被要求观看某两张随机选取的图片，然后小D需要根据他自己主观上的判断确定这两张图片谁好谁坏，或者这两张图片质量差不多。
用符号“<”、“>”和“=”表示图片x和y（x、y为图片编号）之间的比较：如果上下文中x和y是图片编号，则x

题解

首先，题目说的这堆东西，其实就是一棵树。。。
可以考虑将相等的节点缩成一个点，这些点不会影响答案。
定义状态f[i][j]，表示以i为根的子树内有j个<号。
考虑两棵子树的合并。
两棵子树合并时可能会产生一些等号，故可以产生转移式：
f[i][c]=g[i][a]*g[k][b]*C(c,i)*C(i,i+j-c)。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<cmath>
#include<algorithm>
#define ll long long
#define mod 1000000007
#define N 105
using namespace std;
int n,m,pos[N],top[N],fa[N],f[N][N],size[N];
int k,cnt,ans,la[N],ff[N],C[N][N],d[N],q[N];
struct node{int a,b;}map[N];
int add(int a,int b){return (a+=b)>=mod?a-mod:a;}

void adde(int a,int b)
{
  map[++k]=(node){a,b};
  ff[k]=la[a];la[a]=k;d[b]++;
}

void get_C()
{
  for(int i=0;i<=n;i++)C[i][0]=1;
  for(int i=1;i<=n;i++)
    for(int j=1;j<=n;j++)
      C[i][j]=add(C[i-1][j],C[i-1][j-1]);
}

void dfs(int x)
{
  int g[N];f[x][0]=1;size[x]=1;
  for(int a=la[x];a;a=ff[a])
    if(fa[x]!=map[a].b)
    {
      fa[map[a].b]=x;dfs(map[a].b);
      memset(g,0,sizeof(g));
      for(int i=0;i<size[x];i++)
        for(int j=1;j<=size[map[a].b];j++)
          for(int k=max(i,j);k<=i+j;k++)
            g[k]=add(g[k],(ll)f[x][i]*f[map[a].b][j]%mod*C[k][j]%mod*C[j][i+j-k]%mod);
      size[x]+=size[map[a].b];
      for(int i=0;i<size[x];i++)f[x][i]=g[i];
    }
  for(int i=size[x];i;i--)f[x][i]=f[x][i-1];
}

bool pd()
{
  int l=1,r=2;q[1]=0;
  while(l<r)
  {
    int x=q[l];l++;
    for(int a=la[x];a;a=ff[a])
      if(--d[map[a].b]==0)q[r]=map[a].b,r++;
  }
  return r>cnt;
}

int find(int x)
{
  if(!top[x])return x;
  return top[x]=find(top[x]);
}

int main()
{
  int a,b;char ch;
  scanf("%d%d",&n,&m);
  for(int i=1;i<=m;i++)
  {
    scanf("%d %c%d",&a,&ch,&b);
    if(ch=='<')swap(a,b);
    if(ch=='=')top[find(b)]=find(a);
    else fa[a]=b;
  }
  for(int i=1;i<=n;i++)
  {
    int x=find(i);
    if(x==i)pos[i]=++cnt;
    else pos[i]=pos[x];
  }
  for(int i=1;i<=n;i++)
    if(find(i)==i)adde(pos[fa[i]],pos[i]);
  if(!pd()){printf("0\n");return 0;}
  get_C();dfs(0);
  for(int i=0;i<=size[0];i++)ans=add(ans,f[0][i]);
  printf("%d\n",ans);
  return 0;
}