比较实验
时间限制: 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;
}