Problem: Your task is the exact same as for the easy version. But this time, the marmots subtract the village's population P from their random number before responding to Heidi's request.

Also, there are now villages with as few as a single inhabitant, meaning that <image>.

Can you help Heidi find out whether a village follows a Poisson or a uniform distribution?

Input

Same as for the easy and medium versions. But remember that now 1 ≤ P ≤ 1000 and that the marmots may provide positive as well as negative integers.

Output

Output one line per village, in the same order as provided in the input. The village's line shall state poisson if the village's distribution is of the Poisson type, and uniform if the answers came from a uniform distribution.

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
string determine_type(vector<int> a) {
  double P = accumulate(a.begin(), a.end(), 0) / 250.;
  sort(a.begin(), a.end());
  vector<double> Fn(2000), Fp(2000), Fu(2000);
  int cur = 0;
  double supu = 0.0, supp = 0.0;
  for (int i = 0; i < 2000; ++i) {
    Fn[i] = i == 0 ? 0 : Fn[i - 1];
    while (cur < 250 && a[cur] < i) {
      Fn[i] += 1 / 250.;
      ++cur;
    }
    Fp[i] = (i == 0 ? 0 : Fp[i - 1]) + exp(i * log(P) - P - lgamma(i + 1));
    Fu[i] = (i == 0 ? 0 : Fu[i - 1]) + (i > 2 * P ? 0 : 1. / (2 * P + 1));
    supu = max(supu, abs(Fu[i] - Fn[i]));
    supp = max(supp, abs(Fp[i] - Fn[i]));
  }
  return supu < supp ? "uniform" : "poisson";
}
int determine_P(vector<int> a) {
  double P;
  if (determine_type(a) == "poisson") {
    P = accumulate(a.begin(), a.end(), 0) / 250.;
  } else {
    sort(a.begin(), a.end());
    P = (a.front() + a.back()) / 2.0;
  }
  return int(round(P));
}
int call = 0;
string determine_type_minus_P(vector<int> a) {
  sort(a.begin(), a.end());
  int mi = a[0], ma = a[249], ra = ma - mi + 1;
  double d1 = mi + ra * 0.25, d2 = mi + ra * 0.5, d3 = mi + ra * 0.75;
  int q[] = {0, 0, 0, 0};
  for (int i = 0; i < 250; ++i) {
    if (a[i] < d1) {
      q[0] += 1;
    } else if (a[i] >= d1 && a[i] < d2) {
      q[1] += 1;
    } else if (a[i] >= d2 && a[i] < d3) {
      q[2] += 1;
    } else {
      q[3] += 1;
    }
  }
  double H = 0.0;
  for (int i = 0; i < 4; ++i) {
    H += q[i] / 250. * log(q[i] / 250.);
  }
  return (H < -1.25 && q[1] + q[2] - (q[0] + q[3]) < 70 || ra <= 3) ? "uniform"
                                                                    : "poisson";
}
void check() {
  uniform_int_distribution<> Pdistr(10, 1000);
  default_random_engine engine;
  vector<int> arg(250);
  int perror = 0, uerror = 0;
  int pPerror = 0, uPerror = 0;
  int pmPerror = 0, umPerror = 0;
  for (int test = 0; test < 10000; ++test) {
    int P = Pdistr(engine);
    poisson_distribution<> pdistr(P);
    for (int i = 0; i < 250; ++i) {
      arg[i] = pdistr(engine);
    }
    if (determine_type(arg) != "poisson") {
      ++perror;
    }
    if (abs(determine_P(arg) / double(P) - 1) > 0.05) {
      ++pPerror;
    }
    for (int i = 0; i < 250; ++i) {
      arg[i] -= P;
    }
    if (determine_type_minus_P(arg) != "poisson") {
      ++pmPerror;
    }
    uniform_int_distribution<> udistr(0, 2 * P);
    for (int i = 0; i < 250; ++i) {
      arg[i] = udistr(engine);
    }
    if (determine_type(arg) != "uniform") {
      ++uerror;
    }
    if (abs(determine_P(arg) / double(P) - 1) > 0.05) {
      cout << "real=" << P << " "
           << " predicted=" << determine_P(arg) << endl;
      ++uPerror;
    }
    for (int i = 0; i < 250; ++i) {
      arg[i] -= P;
    }
    if (determine_type_minus_P(arg) != "uniform") {
      ++pmPerror;
    }
  }
  cout << call << endl;
  cout << perror << " " << uerror << endl;
  cout << pPerror << " " << uPerror << endl;
  cout << pmPerror << " " << umPerror << endl;
}
int main() {
  int v;
  cin >> v;
  for (int test = 0; test < v; ++test) {
    vector<int> a(250);
    for (int i = 0; i < 250; ++i) {
      cin >> a[i];
    }
    cout << determine_type_minus_P(a) << endl;
  }
}