Random Ranking

Imagine a real contest or exam of n participants. Every participant will get a particular score. We can predict the standings board more or less, if we do some statistics on their previous performance.

Let’s say the score of the participants will be uniformly distributed in interval [l i, r i] (the score can be a real number). Can you predict the standings board according to these data? In other words you should say for each participant the probability that he gets some fixed place in the scoreboard. The participants are sorted by increasing of their scores in the scoreboard. So, the participant with the largest score gets the last place.

Input

The first line contains integer n (1  ≤ n  ≤ 80), showing how many participants we have. Each of the next n lines contains our predictions, the i-th line contains a pair of integers l i, r i (0 ≤ l i < r i ≤ 109) as the distributed interval for participant i.

Consider the participants numbered from 1 to n in some way.

Output

Output a distributed matrix a of order n. The element a ij of the matrix is the probability that participant i has rank j.

Your answer will considered correct if it has at most 10 - 6 absolute or relative error.

Examples

input

2
1 6
4 9

output

0.9200000000 0.080 
0.080 0.9200000000

input

8
0 2
1 3
2 4
3 5
4 6
5 7
6 8
7 9

output

0.875 0.125 0 0 0 0 0 0 
0.125 0.750 0.125 0 0 0 0 0
0 0.125 0.750 0.125 0 0 0 0
0 0 0.125 0.750 0.125 0 0 0
0 0 0 0.125 0.750 0.125 0 0
0 0 0 0 0.125 0.750 0.125 0
0 0 0 0 0 0.125 0.750 0.125
0 0 0 0 0 0 0.125 0.875

Note

The score probability distribution is continuous, which means, there is no possibility for a draw.

Solution:

#include <stdio.h>
#include <memory.h>
#include <vector>
#include <algorithm>

using namespace std;

double ans[88][88], res[88][88][88], fp[88][88][88], fs[88][88][88];

int main() {
int n, l[88], r[88];
scanf("%d", &n);
for (int i=0;i<n;i++) scanf("%d %d", l+i, r+i);
  vector <int> bounds;
for (int i=0;i<n;i++) {
    bounds.push_back(l[i]);
    bounds.push_back(r[i]);
  }
  sort(bounds.begin(), bounds.end());
  for (int i=0;i<n;i++)
    for (int j=0;j<n;j++)
      for (int k=0;k<n;k++) res[i][j][k] = 0.0;
  int sz = bounds.size();
  for (int q=0;q<sz-1;q++) {
    int st = bounds[q], fin = bounds[q+1];
    if (st == fin) continue;
    int m = 0, c[88];
    double pin[88], pbig[88], psmall[88];
    int add = 0;
    for (int i=0;i<n;i++)
      if (l[i] <= st && fin <= r[i]) {
        m++;
        c[m] = i;
        pin[m] = 1.0 * (fin - st) / (r[i] - l[i]);
        pbig[m] = 1.0 * (st - l[i]) / (r[i] - l[i]);
        psmall[m] = 1.0 - pin[m] - pbig[m];
      } else
      if (r[i] <= st) add++;
    memset(fp, 0, sizeof(fp));
    fp[0][0][0] = 1.0;
    for (int i=0;i<m;i++)
      for (int x=0;x<=i;x++)
        for (int y=0;y<=i-x;y++) {
          double ft = fp[i][x][y];
          fp[i+1][x][y] += ft * psmall[i+1];
          fp[i+1][x][y+1] += ft * pin[i+1];
          fp[i+1][x+1][y] += ft * pbig[i+1];
        }
    memset(fs, 0, sizeof(fs));
    fs[m+1][0][0] = 1.0;
    for (int i=m+1;i>1;i--)
for (int x=0;x<=m+1-i;x++)
        for (int y=0;y<=m+1-i-x;y++) {
          double ft = fs[i][x][y];
          fs[i-1][x][y] += ft * psmall[i-1];
          fs[i-1][x][y+1] += ft * pin[i-1];
          fs[i-1][x+1][y] += ft * pbig[i-1];
        }
    for (int t=m/2;t>=1;t--) {
for (int i=t-1;i<m-1;i++) {
        for (int x=0;x<=i+1;x++)
          for (int y=0;y<=i+1-x;y++) fp[i+1][x][y] = 0.0;
        for (int x=0;x<=i;x++)
          for (int y=0;y<=i-x;y++) {
            double ft = fp[i][x][y];
            fp[i+1][x][y] += ft * psmall[i+2];
            fp[i+1][x][y+1] += ft * pin[i+2];
            fp[i+1][x+1][y] += ft * pbig[i+2];
          }
      }
      for (int x=0;x<=m-1;x++)
        for (int y=0;y<=m-1-x;y++) res][x+add][y] += pin[t] * fp[m-1][x][y];
    }
    for (int t=m/2+1;t<=m;t++) {
      for (int i=t+1;i>2;i--) {
for (int x=0;x<=m+2-i;x++)
          for (int y=0;y<=m+2-i-x;y++) fs[i-1][x][y] = 0.0;
        for (int x=0;x<=m+1-i;x++)
          for (int y=0;y<=m+1-i-x;y++) {
            double ft = fs[i][x][y];
            fs[i-1][x][y] += ft * psmall[i-2];
            fs[i-1][x][y+1] += ft * pin[i-2];
            fs[i-1][x+1][y] += ft * pbig[i-2];
          }
      }
      for (int x=0;x<=m-1;x++)
        for (int y=0;y<=m-1-x;y++) res][x+add][y] += pin[t] * fs[2][x][y];
    }
  }
  for (int i=0;i<n;i++)
    for (int j=0;j<n;j++) ans[i][j] = 0.0;
  for (int i=0;i<n;i++)
    for (int x=0;x<n;x++)
      for (int y=0;y<n-x;y++)
        for (int z=0;z<=y;z++) ans[i][x+z] += res[i][x][y] / (y+1);
  for (int i=0;i<n;i++) {
    for (int j=0;j<n;j++) printf("%.6lf ", ans[i][j]);
    printf("\n");
  }
  return 0;
}