Morass

Shopping With Grandma II

Jul 24th, 2016
138
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 1.30 KB | None | 0 0
  1. #include <bits/stdc++.h>
  2. using namespace std;
  3. #define PB push_back
  4. #define ZERO (1e-10)
  5. #define INF (1<<29)
  6. #define CL(A,I) (memset(A,I,sizeof(A)))
  7. #define DEB printf("DEB!\n");
  8. #define D(X) cout<<"  "<<#X": "<<X<<endl;
  9. #define EQ(A,B) (A+ZERO>B&&A-ZERO<B)
  10. typedef long long ll;
  11. typedef double ld;
  12. typedef pair<ll,ll> pll;
  13. typedef vector<int> vi;
  14. typedef pair<int,int> ii;
  15. typedef vector<ii> vii;
  16. #define IN(n) int n;scanf("%d",&n);
  17. #define FOR(i, m, n) for (int i(m); i < n; i++)
  18. #define REP(i, n) FOR(i, 0, n)
  19. #define F(n) REP(i, n)
  20. #define FF(n) REP(j, n)
  21. #define FT(m, n) FOR(k, m, n)
  22. #define aa first
  23. #define bb second
  24. void ga(int N,int *A){F(N)scanf("%d",A+i);}
  25. #define MX (101)
  26. ld dp[2*MX][MX],dd[2*MX][MX];
  27. ld dyn(int s,int b){
  28.     if(!b)return s;
  29.     ld &v(dp[s][b]);
  30.     if(v>-ZERO)return v;
  31.     if(!s)return v=dyn(1,b-1);
  32.     return v=(s*dyn(s-1,b)+b*dyn(s+1,b-1))/(s+b);
  33. }
  34. ld dyp(int s,int b){
  35.     if(!b)return 0;
  36.     ld &v(dd[s][b]);
  37.     if(v>-ZERO)return v;
  38.     if(!s)return v=dyp(1,b-1)+1;
  39.     return v=(s*dyp(s-1,b)+b*dyp(s+1,b-1))/(s+b)+1;
  40. }
  41. int A,B;
  42. int main(void){
  43.     F(2*MX)FF(MX)dp[i][j]=-666;
  44.     F(2*MX)FF(MX)dd[i][j]=-666;
  45.     dyn(MX-1,MX-1);
  46.     while(scanf("%d%d",&A,&B),A|B){
  47.         printf("%.12lf %.12lf\n",A?dp[B][A]:B,dyp(B,A));
  48.     }
  49.     return 0;
  50. }
Add Comment
Please, Sign In to add comment