题意就不说了,我们也错了那么多次,这个题有两个问,其实比较难的是第一个问,但是刚好第一个问我在上一周曾经问过龙哥,那么第一个问就秒出答案,对于第一个问的精度问题,当暴力算时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;
}