方程很简单，每一公里往上推就行

WA了2发，忘了单通道时的特判，还有n m傻傻分不清，忘了fixed什么的我好弱啊QAQ..

#include
    
     #define rep(i,j,k) for(int i=j;i<=k;i++)using namespace std;double dp[2][30010];vector
     
       vec[1003];int n,m,p,k,a,b;int main(){    ios::sync_with_stdio(0);    while(cin>>m>>k>>n>>p){        memset(dp,0,sizeof dp);dp[0][p]=1;        memset(vec,0,sizeof vec);        rep(i,1,k){            cin>>a>>b;            vec[b].push_back(a);        }        rep(i,1,n) sort(vec[i].begin(),vec[i].end());        vector
      
       ::iterator it;        if(m>=2) rep(i,0,n-1){            memset(dp[i+1&1],0,sizeof dp[i+1&1]);            rep(j,1,m) if( ( it = find(vec[i+1].begin(),vec[i+1].end(),j) ) !=vec[i+1].end() ) dp[i+1&1][j]=-6666666;            rep(j,1,m){                if(dp[i&1][j]<=0) continue;                if(j==1){                    if(dp[i+1&1][j+1]>=0) dp[i+1&1][j+1]+=1.0/2.0*dp[i&1][j];                    if(dp[i+1&1][j]>=0) dp[i+1&1][j]+=1.0/2.0*dp[i&1][j];                }                else if(j==m){                    if(dp[i+1&1][j-1]>=0) dp[i+1&1][j-1]+=1.0/2.0*dp[i&1][j];                    if(dp[i+1&1][j]>=0) dp[i+1&1][j]+=1.0/2.0*dp[i&1][j];                }                else{                    if(dp[i+1&1][j]>=0) dp[i+1&1][j]+=1.0/3.0*dp[i&1][j];                    if(dp[i+1&1][j-1]>=0) dp[i+1&1][j-1]+=1.0/3.0*dp[i&1][j];                    if(dp[i+1&1][j+1]>=0) dp[i+1&1][j+1]+=1.0/3.0*dp[i&1][j];                }            }        }        else{            bool flag=0;            rep(i,1,n)  if((it=find(vec[i].begin(),vec[i].end(),1))!=vec[i].end()) flag=1;            if(flag==1) cout<<"0.000000"<