2017川大校赛A题

题意就不说了，我们也错了那么多次，这个题有两个问，其实比较难的是第一个问，但是刚好第一个问我在上一周曾经问过龙哥，那么第一个问就秒出答案，对于第一个问的精度问题，当暴力算时n是小于10的，答案要求保留前6位小数，没问题。然后第二个问就是一个比较简单的概率DP。

题目求抽m次刚好抽到n张卡的概率，那么对于第m-1次，其实只有n-1张卡，这是一个关键的问题。

令一个dp[i][j] 表示抽 i 次抽到 j 张卡的概率，那么最终的答案就是  dp[m-1][n-1]*(1/n)；这个就比较好理解了。

为什么用这个状态来表示，因为对于这个题，如果令dp为刚好抽到的概率其实很难，不好推出转移方程，那么用这个比较容易求出的子问题，其实更简单。

具体转移我就不写了，都比较简单。

注意最后一些特殊数据的处理，毕竟n,m都可能等于0.

#include<stdio.h>
#include<algorithm>
#include<string>
#include<string.h>
#include<queue>
#include<vector>
#include<stack>
#include<math.h>
#include<map>
#include<iostream>
using namespace std;
#define maxn 1005
int n,m;
double dp[50][50];
int main()
{
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        double ans1=0.0;
        for(int i=1;i<=n;i++)
            ans1+=(double)n/(double)i;
        dp[1][1]=1.0;
        dp[0][0]=1.0;
        for(int i=1;i<=m-1;i++)
        {
            for(int j=1;j<=n-1;j++)
            {
                if(i==1&&j==1)
                    continue;
                dp[i][j]=dp[i-1][j-1]*(double)(n-j+1)/(double)n;
                dp[i][j]+=dp[i-1][j]*(double)j/(double)n;
            }
        }
        if(n==0&&m==0)
            printf("%.6f %.6f\n",ans1,1.0);
        else if(n==0&&m>0)
             printf("%.6f %.6f\n",ans1,0.0);
        else
            printf("%.6f %.6f\n",ans1,dp[m-1][n-1]/(double)n);
    }
return 0;
}