ZOJ 3604 Help Me Escape（期望dp）

题目传送门

题意： 有n条道路，每条道路有一个危险值$c[i]$,你有一个初始战斗力$f$，如果你的战斗力大于一条道路的战斗力，那么你需要花$\frac{1+\sqrt{5}}{2}\ast c[i{]}^2$天逃离，否则，你每天战斗力增加$c[i]$。现在每天给你随机扔到一条道路上，直到你逃离这条路再继续下一条路。问你逃离所有道路的期望天数。

思路： $dp[f]$表示战斗力为$f$时，还需要$dp[f]$天才能逃离，直接记忆化搜索即可。

代码：

#include<bits/stdc++.h>
#define endl '\n'
#define null NULL
#define ls p<<1
#define rs p<<1|1
#define fi first
#define se second
#define mp make_pair
#define pb push_back
#define ll long long
#define int long long
#define lowbit(x) x&-x
#define pii pair<int,int>
#define ull unsigned long long
#define pdd pair<double,double>
#define sz(x) (int)(x).size()
#define all(x) (x).begin(),(x).end()
#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
char *fs,*ft,buf[1<<20];
#define gc() (fs==ft&&(ft=(fs=buf)+fread(buf,1,1<<20,stdin),fs==ft))?0:*fs++;
inline int read()
{
    int x=0,f=1;
    char ch=gc();
    while(ch<'0'||ch>'9')
    {
        if(ch=='-')
            f=-1;
        ch=gc();
    }
    while(ch>='0'&&ch<='9')
    {
        x=x*10+ch-'0';
        ch=gc();
    }
    return x*f;
}
using namespace std;
const int N=4e2+666;
const int inf=0x3f3f3f3f;
const int mod=1e9+7;
const double eps=1e-7;
const double PI=acos(-1);

double dp[20010];
int c[N],n;

double dfs(int f)
{
    if(dp[f]>0)
        return dp[f];
    dp[f]=0;
    for(int i=1;i<=n;i++)
    {
        if(f>c[i])
        {
            double temp=(1.0+sqrt(5))/2*c[i]*c[i];
            int t = (int)temp;
            dp[f]+=(double)t/n;
        }
        else
        {
            dp[f]+=(dfs(f+c[i])+1)/n;
        }
    }
    return dp[f];
}

int solve()
{
    int f;
    while(cin>>n>>f)
    {
        memset(dp,0,sizeof dp);
        for(int i=1; i<=n; i++)
        {
            cin>>c[i];
        }
        printf("%.3lf\n",dfs(f));
    }
    return 0;
}

signed main()
{
//    int _;
//    cin>>_;
//    while(_--)
        solve();
    return 0;
}