nbu 2416 奇怪的散步

题目链接：http://acm.nbu.edu.cn/v1.0/Problems/Problem.php?pid=2416

题目大意：

有一个色子，这个色子掷出1~K的概率都相等。每次掷出来是几点就往前走几步。请问当扔了N次色子后，正好走了L步的概率是多少？

1<=K<=20

1<=N<=100

1<=L<=2000

题目思路：

题目要求扔了N次骰子，走了L步的概率。

设这个概率为p[N][L]

那么我们知道p[N][L]=p[N-1][L-1]+p[N-1][L-2]+...+p[N-1][max(1,L-K)]。

显然p[0][0]=1,所以上面的递推公式，我们很容易求出p[N][L]。

复杂度为O(N*L*K)

代码：

#pragma comment(linker, "/STACK:102400000,102400000")
#include<stdio.h>
#include<string.h>
#include<math.h>
#include<stdlib.h>
#include<ctype.h>
#include<iostream>
#include<algorithm>
#include<stack>
#include<queue>
#include<map>
#include<set>
#include<vector>
#include<string>
using namespace std;
#define ll long long
#define clr(x,c,n) memset(x,c,sizeof(x[0])*(n))
#define clr_all(x,c) memset(x,c,sizeof(x))
#define IT iterator
#define ls rt<<1
#define rs ls|1
#define lson l,mid,ls
#define rson mid+1,r,rs
#define middle l+r>>1
#define MOD 1000000007
#define inf 0x3f3f3f3f
#define eps (1e-8)
#define PI 3.1415926535897932384626433832795
#define E 2.7182818284590452353602874713527
template <class T> T _min(T a,T b){return a<b? a:b;}
template <class T> T _max(T a,T b){return a>b? a:b;}
template <class T> T _abs(T a){return a>0? a:-a;}
template <class T> T _mod(T a,T m){return a<m? (a<0? (a%m+m)%m:a):a%m;}
template <class T> T _gcd(T a,T b){while(b){T t=b;b=a%b;a=t;}return a;}
template <class T> void _swap(T &a,T &b){T t=b;b=a;a=t;}
template <class T> void getmax(T &a,T b){a= a>b? a:b;}
template <class T> void getmin(T &a,T b){a= (a!=-1 && a<b)? a:b;}
int TS,cas=1;
const int M=2000+5;
int k,n,l;
double dp[2][M];

void run(){
	int i,j,o;
	int now=0;
	for(i=1;i<=l;i++) dp[0][i]= i<=k? 1.0/(1.0*k):0.0;
	for(i=2;i<=n;i++){
		now^=1;
		for(j=1;j<=l;j++) dp[now][j]=0.0;
		for(j=1;j<=l;j++)
			for(o=1;o<=k && o<j;o++)
				dp[now][j]+=dp[now^1][j-o]/(1.0*k);
	}
	printf("%.3lf\n",dp[now][l]);
}

void presof(){
}

int main(){
	//freopen("input.txt","r",stdin);
    //freopen("output.txt","w",stdout);
	presof();
	//run();
	while(~scanf("%d%d%d",&k,&n,&l)) run();
	//for(scanf("%d",&TS),cas=1;cas<=TS;cas++) run();
	return 0;
}