CC和他的AE86

CC和他的AE86

题意：

概率dp

场上本来是用矩阵写的，无奈出了点偏差一直过不了样例！！！

#include <bits/stdc++.h>
using namespace std;

const int maxn=55;
const int maxe=3600;
struct Matrix{
    int n;
    double m[maxn][maxn];

    void init(int sz){
        n=sz;
        for(int i=0;i<n;i++){
            for(int j=0;j<n;j++){
                m[i][j]=0.0;
            }
        }
    }
    Matrix(int sz){init(sz);}

    void set_I(){
        for(int i=0;i<n;i++) m[i][i]=1.0;
    }

    Matrix operator* (const Matrix &a)const{
        Matrix ans(n);
        for(int k=0;k<n;k++){
            for(int i=0;i<n;i++){
                for(int j=0;j<n;j++){
                    ans.m[i][j]+=m[i][k]*a.m[k][j];
                }
            }
        }
        return ans;
    }
};

int deg[maxn];
int u[maxe],v[maxe];


int main(){
    int n,m;
    while(scanf("%d%d",&n,&m)!=EOF){
        memset(deg,0,sizeof(deg));
        for(int i=0;i<m;i++){
            scanf("%d%d",&u[i],&v[i]);
            u[i]--;v[i]--;
            deg[u[i]]++;
            deg[v[i]]++;
        }
        int q;
        scanf("%d",&q);
        while(q--){
            int ss,ee,d;
            scanf("%d%d%d",&ss,&ee,&d);
            ss--;ee--;
            Matrix base(n);
            for(int i=0;i<m;i++){
                int x=u[i],y=v[i];
                base.m[x][y]=1.0/deg[x];
                base.m[y][x]=1.0/deg[y];
                if(x==ee) base.m[x][y]=0;
                if(y==ee) base.m[y][x]=0;
            }
            Matrix ans(n);
            ans.set_I();
            while(d){
                if(d&1) ans=ans*base;
                d>>=1;
                base=base*base;
            }
            double res=0;
            for(int i=0;i<n;i++){
                if(i!=ee) res+=ans.m[ss][i];
            }
            printf("%.8lf
",1.0-res);
        }
    }
}

矩阵好像慢很多，下面dalao写的还没仔细看~

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int maxn=55;
const int maxm=1005;
const int maxe=7000005;
const int inf=0x3f3f3f3f;
int n,m,k;
double dp[maxn][maxn][105];
vector<int> e[maxn];
bool vis[maxn][maxn][105];
double DP(int u,int to,int times){
    if(times==-1)return 0;
    if(vis[u][to][times])return dp[u][to][times];
    if(u==to)return 1;
    int sz=e[u].size();
    double ans=0;
    for(int i=0;i<sz;i++){
        ans+=DP(e[u][i],to,times-1)/sz;
    }
    vis[u][to][times]=1;
    return dp[u][to][times]=ans;
}

int main(){
    int a,b,q,c;
    while(~scanf("%d%d",&n,&m)){
            for(int i=1;i<=n;i++)e[i].clear();
            memset(vis,0,sizeof(vis));
            for(int i=0;i<m;i++){
                scanf("%d%d",&a,&b);
                e[a].push_back(b);
                e[b].push_back(a);
            }
            scanf("%d",&q);
            for(int i=0;i<q;i++){
                scanf("%d%d%d",&a,&b,&c);
                printf("%.10lf
",DP(a,b,c));
            }
    }
    return 0;
}

哎，前三场个人赛全部爆炸！每场都会卡在一个小bug上无法自拔~~ 今天这场还差点爆零=_=

明天还有最后一场，但愿不爆零吧=_=

其实一个假期下来感觉已经提高不少了，一比才知道差距啊。。。。。。