题目大意:
已知:f0=1−1e,fn=1−nfn−1f_0=1-\frac1e,f_n=1-nf_{n-1}f0=1−e1,fn=1−nfn−1
求fn,n≤104f_n,n\le10^4fn,n≤104精度要求10−410^{-4}10−4。
题解:
注意到fnf_nfn可以写成an−bne−1a_n-b_ne^{-1}an−bne−1的形式,a0=1,b0=1a_0=1,b_0=1a0=1,b0=1。
经过一些推导化简可以显然得知:
bn=(−n)bn−1=(−1)nn!an=∑i=0nn!(n−i)!(−1)i=n!(−1)n∑i=0n(−1)ii!b_n=(-n)b_{n-1}=(-1)^nn!\\
a_n=\sum_{i=0}^n\frac{n!}{(n-i)!}(-1)^i=n!(-1)^n\sum_{i=0}^n\frac{(-1)^i}{i!}
bn=(−n)bn−1=(−1)nn!an=i=0∑n(n−i)!n!(−1)i=n!(−1)ni=0∑ni!(−1)i
显然limitn→+inf∑i=0n(−1)ii!=e−1\mathrm{limit}_{n\rightarrow+\inf}\sum_{i=0}^n\frac{(-1)^i}{i!}=e^{-1}limitn→+infi=0∑ni!(−1)i=e−1
即limitn→+inffn=n!(−1)ne−1−n!(−1)ne−1=0\mathrm{limit}_{n\rightarrow+\inf}f_n=n!(-1)^ne^{-1}-n!(-1)^ne^{-1}=0limitn→+inffn=n!(−1)ne−1−n!(−1)ne−1=0
(顺带由刚刚的推导可以看出误差大约在O(1n)O\left(\frac1n\right)O(n1)左右的样子)
因此设fN=0f_{N}=0fN=0,并取N≥105N\ge10^5N≥105,然后反着推回来即可。
#include<bits/stdc++.h>
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define Rep(i,v) rep(i,0,(int)v.size()-1)
#define lint long long
#define ull unsigned lint
#define db long double
#define pb push_back
#define mp make_pair
#define fir first
#define sec second
#define gc getchar()
#define debug(x) cerr<<#x<<"="<<x
#define sp <<" "
#define ln <<endl
using namespace std;
typedef pair<int,int> pii;
typedef set<int>::iterator sit;
inline int inn()
{
int x,ch;while((ch=gc)<'0'||ch>'9');
x=ch^'0';while((ch=gc)>='0'&&ch<='9')
x=(x<<1)+(x<<3)+(ch^'0');return x;
}
const int N=500010;
db f[N];
int main()
{
int n=inn();f[N-1]=0;
for(int i=N-2;i>=0;i--)
f[i]=(1-f[i+1])/(i+1);
printf("%.4f\n",double(f[n]));
return 0;
}