电路稳定性(cir.pas/c/cpp)（栈的应用）

电路稳定性(cir.pas/c/cpp)

【问题描述】
有一个电路，上面有 n 个元件。已知 i 损坏耳断开的概率是 pi(i=1…n,0≤pi≤1)。求电路
断路的概率。
元件的连接方式十分简单，对电路表示如下：
1、 一个元件是最小的电路，用字母表里的第 i 个字母表示元件 i；
2、 k 个电路组成串联电路表示为：电路 1,电路 2,……,电路 k。
3、 k 个电路组成并联电路表示为： (电路 1)(电路 2)……(电路 k)。
【输入文件】
输入文件第一行是一个整数 n（ 1≤n≤26），表示共有 n 个元件；第二行是表示电路
的字符串；最后 n 行，每行有一个实数 pi(i=1…n,0≤pi≤1)，表示该元件的断路的概率。
【输出文件】
输出文件只有一个实数，表示整个电路的断路概率，精确到小数点后四位。
【文件样例】
cir.in cir.out
5 0.2992
(A,B)((C)(D),E)
0.2
0.3
0.4
0.5
0.6
【思路】
这个题相当于一个模拟题，其实就是把符号存入符号栈，符号的位置存入位置栈1，数字存入数字栈，数字的位置存入位置栈2。然后每当遇到一个“）”的时候，就往前进行清除处理，对上一个“(”到这一个“）”之间的元素进行处理，这中间肯定是只有数字和符号“，”组成，所以进行串联的话就是1-∏（1-Pi）的概率。
然后如果遇见连续的数字的话（那些括号已经被处理过但还留下数字的数字），就表明它是并联的，对应概率就应该是∏Pi。
当然了，还有一些细节需要注意
while(npos>pos[top2]&&npos&&top2>1)
表示的是当前处理到的数字在字符串上的位置如果在尚未处理的最后一个符号在字符串上的位置的后面，就说明我们需要处理这个数。
然后就是输出精度的问题了。具体代码感觉写的挺清楚了。
这个题我用了四个手打栈，据说可以不用记录位置的两个栈的，还请众多dalao赐教。

【代码】

#include<iostream>
#include<cstring>
#include<cmath>
using namespace std;
int n,m,i,j;
double a[10001];
string s;
double num[10001];
int top1;
int sym[10001],top2;
int po[10001],pos[10001];


int r()
{
    int ans=0,f=1;
    char ch=getchar();
    while(ch<'0'||ch>'9')
    {
//      if(ch=='-')
//      f=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
    {
        ans*=10;
        ans+=ch-'0';
        ch=getchar();
    }
    return ans*f;
}

void clearr(int x)
{
    int npos=x;
    int kk=0;
    while(npos>pos[top2]&&npos&&top2>1)
    {
        double ans=num[top1];
        while(sym[top2]!=1&&top1>1&&pos[top2]<po[top1-1])
        {
            ans=(num[top1])*(num[top1-1]);
            num[--top1]=ans;
        }
        while(sym[top2]!=1)
        {
            ans=(1-num[top1-1])*(1-num[top1]);
            num[--top1]=1-ans;
            top2--;
//          kk=1;
        }
//      if(kk)
//      kk=0;
//      else
        top2--;

        if(po[top1-1]>pos[top2]+1&&sym[top2]!=2)
        {
            top1--;
            num[top1]*=num[top1+1];
        }
        npos=po[top1-1];
    }
}

int main()
{
//  freopen("cir.in","r",stdin);
//  freopen("cir.out","w",stdout);
    n=r();
    cin>>s;
    for(i=1;i<=n;i++)
    scanf("%lf",&a[i]);
    int len=s.size();
    sym[++top2]=1;
    pos[top2]=-1;
    s[len]=')';
    for(i=0;i<=len;i++)
    {
        if(s[i]=='(')
        sym[++top2]=1,pos[top2]=i;
        else if(s[i]==')')
        {
            clearr(i-1);
        }
        else if(s[i]==',')
        sym[++top2]=2,pos[top2]=i;
        else
        {num[++top1]=a[s[i]-'A'+1];
        po[top1]=i;}
    }
    double aans=1;

    printf("%.4f",num[1]);
    return 0;
}
/*
num 0.44 0.4 
pos 4 8
sym ( ( 
pos -1 6 7

6
(A,B)(E,(C)(D),F)
0.2
0.3
0.4
0.5
0.6
0.8
*/