poj 1217 FOUR QUARTERS

# include <cstdio># include <cstring># include <iostream># include <map>using namespace std;# define enpre(a,b) (((a)<<6)|(b))# define getp1(a) ((a)>>6)# define getp2(a) ((a)&63)int main(){     double res[21][2];     //p[i][j]，A得到i-1分，B得到j-1分的概率，以后为了处理方便，可以将得-1,0,1,2分看作得0,1,2,3分     double p[4][4]={{0,0,0,1.0/16},{1.0/16,3.0/8,1.0/8,0},{1.0/8,3.0/16,0,0},{1.0/16,0,0,0}};     memset(res,0,sizeof(res));     map<int,double> table[21];     table[0][enpre(0,0)]=1.0;//边界情况：A得0分，B得0分的概率为1     for(int step=1;step<=21;step++)     {         for(map<int,double>::iterator ita=table[step-1].begin();ita!=table[step-1].end();ita++)         {  int p1=getp1(ita->first),p2=getp2(ita->first);  if(step!=21)  {      for(int i=0;i<4;i++)          for(int j=0;j<4;j++)   table[step][enpre(p1+i,p2+j)]+=(ita->second)*p[i][j];//状态转移  }  if(p1<p2)//累加A输给B的概率      res[step-1][0]+=ita->second;  else if(p1==p2)//累加打平的概率      res[step-1][1]+=ita->second;         }     }     printf("Round   A wins    B wins    Tien");     for(int i=1;i<=20;i++)//输出格式注意,如果需要在printf内输出%的话要用%%         printf("%5d%10.4f%%%9.4f%%%9.4f%%n",i,(1-res[i][0]-res[i][1])*100,res[i][0]*100,res[i][1]*100);     return 0; }