UVa 10288 - Coupons(数学期望 + 递推)
链接:
https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=1229
题意:
大街上到处在卖彩票,一元钱一张。购买撕开它上面的锡箔,你会看到一个漂亮的图案。
图案有n种,如果你收集到所有n(n≤33)种彩票,就可以得大奖。
请问,在平均情况下,需要买多少张彩票才能得到大奖呢?如n=5时答案为137/12。
分析:
已有k个图案,令s=k/n,拿一个新的需要t次的概率:(s^(t-1))(1-s);
因此平均需要的次数为(1-s)(1 + 2s + 3s^2 + 4s^3 + …) = (1-s)E,
而sE = s + 2s^2 + 3s^3 + … = E-(1+s+s^2+…),移项得(1-s)E = 1+s+s^2+… = 1/(1-s) = n/(n-k)
换句话说,已有k个图案:平均拿n/(n-k)次就可多搜集一个,所以总次数为:n(1/n+1/(n-1)+1/(n-2)+…+1/2+1/1)
代码:
1 import java.io.*; 2 import java.util.*; 3 4 public class Main { 5 Scanner cin = new Scanner(new BufferedInputStream(System.in)); 6 7 long gcd(long a, long b) { 8 return b == 0 ? a : gcd(b, a%b); 9 } 10 11 long lcm(long a, long b) { 12 return a / gcd(a,b) * b; 13 } 14 15 int len(long v) { 16 int res = 0; 17 do { v /= 10; res++; } while(v > 0); 18 return res; 19 } 20 21 void printCh(char c, int x) { 22 while(x --> 0) System.out.print(c); 23 } 24 25 void output(long a, long b, long c) { 26 if(b == 0) { System.out.println(a); return; } 27 printCh(' ', len(a)+1); System.out.println(b); 28 System.out.print(a + " "); printCh('-', len(c)); System.out.println(); 29 printCh(' ', len(a)+1); System.out.println(c); 30 } 31 32 void MAIN() { 33 while(cin.hasNext()) { 34 int n = cin.nextInt(); 35 long b = 0, c = 1; 36 for(int i = 2; i <= n; i++) c = lcm(c, i); 37 for(int i = 0; i < n; i++) b += c / (n-i) * n; 38 output(b/c, b%c/gcd(b%c,c), c/gcd(b%c,c)); 39 } 40 } 41 42 public static void main(String args[]) { new Main().MAIN(); } 43 }