POJ[3744]Scout YYF I 矩阵乘法

放传送门的同时告诉大家我是SB

可以轻松推出$f_i=f_{i-1}*p+f_{i-1}*(1-p)$
没雷就矩乘转移
有雷就矩乘到雷的前面，再跳一步
问题迎刃而解
！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！

输出用%f！！！！！！起点有雷要特判！！！！！！调试语句要删！！！！！！读入优化就是傻逼！！！！！！

代码：

#include<algorithm>
#include<iostream>
#include<ctype.h>
#include<cstdio>
#define N 25
using namespace std;
inline int read(){
    int x=0,f=1;char c;
    do c=getchar(),f=c=='-'?-1:f; while(!isdigit(c));
    do x=(x<<3)+(x<<1)+c-'0',c=getchar(); while(isdigit(c));
    return x*f;
}///读入优化sb东西！！！！！！
double f[3];
double p;
int n,pre,t,a[N];
bool flag;
struct Matrix{
    double a[3][3];
    inline void mul(Matrix);
}q,b,sum;
inline void Matrix::mul(Matrix c){
    Matrix b;
    b.a[1][1]=a[1][1]*c.a[1][1]+a[1][2]*c.a[2][1];
    b.a[1][2]=a[1][1]*c.a[1][2]+a[1][2]*c.a[2][2];
    b.a[2][1]=a[2][1]*c.a[1][1]+a[2][2]*c.a[2][1];
    b.a[2][2]=a[2][1]*c.a[1][2]+a[2][2]*c.a[2][2];
    a[1][1]=b.a[1][1];a[1][2]=b.a[1][2];
    a[2][1]=b.a[2][1];a[2][2]=b.a[2][2];
    return;
}
inline void init(){
    f[1]=0;f[2]=1;
    q.a[1][1]=0;q.a[1][2]=1-p;
    q.a[2][1]=1;q.a[2][2]=p;
    pre=1;
    flag=false;
}
inline Matrix Pow(Matrix a,int k){
    sum=a;
    k--;
    while(k){
        if(k&1) sum.mul(a);
        a.mul(a);
        k>>=1;
    }
    return sum;
}
inline void mul(Matrix a){
    double t,t1;
    t=f[1]*a.a[1][1]+f[2]*a.a[2][1];
    t1=f[1]*a.a[1][2]+f[2]*a.a[2][2];///蒟蒻这里没赋中间变量直接转移跪爆
    f[1]=t;
    f[2]=t1;
}
int main(){
    while(scanf("%d%lf",&n,&p)==2){
        init();///pre代表当前f[2]代表的位置
        for(int i=1;i<=n;i++){
            scanf("%d",&a[i]);
            if(a[i]==1) flag=true;///特判无解
        }
        sort(a+1,a+n+1);
        for(int i=1;i<n;i++)
            if(a[i]+1==a[i+1]) flag=true;///特判无解*2
        if(flag || a[1]==1){
            printf("0.0000000\n");
            continue;
        }
        for(int i=1;i<=n;i++){
            t=a[i]-1-pre;
            if(t) mul(Pow(q,t));
            pre=a[i]+1;
            f[1]=0;f[2]=f[2]*(1-p);
        }
        printf("%.7f\n",f[2]);///不要用%lf输出！
    }
return 0;
}