2004 TCO QR3 div1 1000（简单条件概率）

题意：
首先给了个基本遗传模型。
string p1a, p1b, p2a, p2b 分别是父母的各自两条DNA链。
然后依赖关系，i depends on j 表示 i 只能在 j 呈显性的时候表现出显性。依赖关系没有环。
求出，孩子每对基因表现显隐的概率。
思路：
处理依赖，就用类似拓扑的方法就行了。。

int n;
double P[55][2];
vector<int> g[55];
class GeneticCrossover
{
        public:
        double cross(string p1a, string p1b, string p2a, string p2b, vector <int> dom, vector <int> rec, vector <int> dependencies)
        {
        n = p1a.length();
        queue<int> q;
        for (int i=0;i<n;++i) g[i].clear();
        for (int i=0;i<n;++i) {
            if (dependencies[i] != -1) {
                g[dependencies[i]].push_back(i);
            } else {
                q.push(i);
            }
        }

        while (!q.empty()) {
            int fr = q.front();q.pop();
            double px, py;
            py = ( (islower(p1a[fr]) != 0) + (islower(p1b[fr]) != 0) ) / 2.0 * ( (islower(p2a[fr]) != 0) + (islower(p2b[fr]) != 0) ) / 2.0;
            px = 1.0 - py;
            if (dependencies[fr] == -1) {
                P[fr][1] = px, P[fr][0] = py;
            } else {
                int u = dependencies[fr];
                P[fr][1] = px * P[u][1];
                P[fr][0] = py + px * P[u][0];
            }
            for (int i=0;i<g[fr].size();++i) q.push(g[fr][i]);
        }
        double ret = 0;
        for (int i=0;i<n;++i) ret += P[i][1] * dom[i] + P[i][0] * rec[i];
        return ret;
        }
}