【POJ 3744】 Scout YYF I

【题目链接】

http://poj.org/problem?id=3744

【算法】

概率DP + 矩阵乘法

【代码】

#include <algorithm>
#include <bitset>
#include <cctype>
#include <cerrno>
#include <clocale>
#include <cmath>
#include <complex>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <ctime>
#include <deque>
#include <exception>
#include <fstream>
#include <functional>
#include <limits>
#include <list>
#include <map>
#include <iomanip>
#include <ios>
#include <iosfwd>
#include <iostream>
#include <istream>
#include <ostream>
#include <queue>
#include <set>
#include <sstream>
#include <stdexcept>
#include <streambuf>
#include <string>
#include <utility>
#include <vector>
#include <cwchar>
#include <cwctype>
#include <stack>
#include <limits.h>
using namespace std;

int i,n;
double p,g,ans;
double x[];

struct Matrix
{
        double mat[][];
} m;

inline Matrix multipy(Matrix &a,Matrix b)
{
        int i,j,k;
        Matrix res;
        memset(res.mat,,sizeof(res.mat));
        for (i = ; i <= ; i++)
        {
                for (j = ; j <= ; j++)
                {
                        for (k = ; k <= ; k++)
                        {
                                res.mat[i][j] += a.mat[i][k] * b.mat[k][j];
                        }
                }
        }
        return res;
}
inline Matrix power(Matrix a,int n)
{
        int i,j;
        Matrix res,b = a;
        for (i = ; i <= ; i++)
        {
                for (j = ; j <= ; j++)
                {
                        res.mat[i][j] = (i == j);
                }
        }
        while (n)
        {
                if (n & ) res = multipy(res,b);
                b = multipy(b,b);
                n >>= ;
        }
        return res;
}

int main()
{

        while (scanf("%d%lf",&n,&p) != EOF)
        {
                for (i = ; i <= n; i++) scanf("%lf",&x[i]);
                sort(x+,x+n+);
                m.mat[][] = p; m.mat[][] =  - p;
                m.mat[][] = ; m.mat[][] = ;
                ans = 1.0;
                x[] = ;
                for (i = ; i <= n; i++)
                {
                        m = power(m,x[i]-x[i-]-);
                        g = m.mat[][];
                        ans *= ( - g);
                        m.mat[][] = p; m.mat[][] =  - p;
                        m.mat[][] = ; m.mat[][] = ;
                }
                printf("%.7lf\n",ans);
        }

        return ;

}