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 ll = long long;
using ld = long double;
using namespace std;
const int MAXN = 4001;
const int C = 250;
const int P = 1000;
const ld EE = 2.7182818284590452353602874713527;
const ld INF = 1e+20;
struct Array {
  ld a[P + MAXN];
  Array() {
    for (int i = 0; i < P + MAXN; i++) a[i] = -INF;
  }
  ld& operator[](int x) { return a[x + P + 3]; }
};
Array p[2][P + 2];
int x[C];
ld fact[MAXN];
ld Poisson(int p, int k) { return logl(p) * k - logl(EE) * p - fact[k]; }
void PreCalc() {
  for (int i = 1; i <= P; i++) {
    for (int j = 0; j < MAXN; j++) {
      ld pp = Poisson(i, j);
      p[1][i][j - i] = Poisson(i, j);
    }
  }
  for (int i = 1; i <= P; i++) {
    for (int j = 0; j <= 2 * i; j++) {
      ld lg = -logl(2 * i + 1);
      p[0][i][j - i] = lg;
    }
  }
}
void Solve() {
  ld mx = -INF;
  int tp = -1;
  int cn = -1;
  for (int i = 0; i < C; i++) {
    cin >> x[i];
    x[i] = min(x[i], MAXN - 1);
  }
  for (int t = 0; t < 2; t++) {
    for (int i = 1; i <= P; i++) {
      ld c = 0;
      for (int j = 0; j < C; j++) {
        c += p[t][i][x[j]];
      }
      if (c > mx) mx = c, tp = t, cn = i;
    }
  }
  if (tp == 0)
    cout << "uniform\n";
  else
    cout << "poisson\n";
}
int main() {
  ios_base::sync_with_stdio(false);
  cout.setf(ios::fixed);
  cout.precision(20);
  fact[0] = 0;
  for (int i = 1; i < MAXN; i++) fact[i] = fact[i - 1] + logl(i);
  PreCalc();
  int v;
  cin >> v;
  for (int i = 0; i < v; i++) Solve();
  return 0;
}