程式師世界 >> 編程語言 >> C語言 >> C++ >> 關於C++ >> rockethon2015 C題 Second price auction 概率dp

rockethon2015 C題 Second price auction 概率dp

編輯：關於C++

題意：n個人去競拍一件商品，下面給出n個區間表示每個人出的價是區間中隨機的一個數（概率均等)則第一名需要付的錢是第二名的競拍價格（允許並列第一名）求支付的錢的期望。

/*********************************************************
  file name: C.cpp
  author : kereo
  create time:  2015年02月08日 星期日 22時26分08秒
*********************************************************/
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
using namespace std;
typedef long long ll;
const int sigma_size=26;
const int N=100+50;
const int MAXN=10000+50;
const int inf=0x3fffffff;
const double eps=1e-8;
const int mod=1000000000+7;
#define L(x) (x<<1)
#define R(x) (x<<1|1)
#define PII pair
#define mk(x,y) make_pair((x),(y))
int n;
int l[MAXN],r[MAXN];
double solve1(int x){ //第一 >x
    int flag,tag;
    double ans=0;
    for(int i=1;i<=n;i++){
        if(r[i]<=x)
            continue;
        flag=1,tag=0;
        for(int j=1;j<=n;j++){
            if(i == j)
                continue;
            if(l[j]>x){
                flag=0;
                break;
            }
            if(l[j]<=x && x<=r[j])
                tag=1;
        }
        if(flag == 0 || tag == 0)
            continue;
        double res=(r[i]-max(x,l[i]-1))*1.0/(r[i]-l[i]+1);
        double ans1=1.0,ans2=1.0; //ans1為全部<=x的概率,ans2為全部=l[j] && x<=r[j]){
                ans1*=(x-l[j]+1)*1.0/(r[j]-l[j]+1);
                ans2*=(x-l[j])*1.0/(r[j]-l[j]+1);
            }
        }
        ans+=(ans1-ans2)*res;
    }
    return ans;
}
double solve2(int x){ //第一 =x
    int flag=0;
    for(int i=1;i<=n;i++)
        if(xx
            ans+=solve2(i)*i; //第一人出價=x
        }
        printf(%.10f
,ans);
    }
	return 0;
}