3
1 1 5
1 3 42
1 2 42
 2
42 23
1 2 10
 10
45792697 27183921 4791060 77070594 27286110 66957065 19283638 2944837 26775511 65379701
5 8 25
8 10 94
1 4 14
9 7 96
4 2 16
2 6 58
2 3 99
6 8 21
8 9 87
 5
3 1 2 6 5
1 2 3
2 3 1
4 3 9
5 3 1
 10
26 31 21 57 56 44 54 58 79 97
2 1 277
1 8 91
1 7 380
4 1 649
10 5 799
10 1 940
8 6 627
6 3 177
9 6 685
 10
8 100 62 38 26 51 26 23 77 32
6 8 29
4 1 82
7 8 66
7 10 67
8 9 98
6 1 14
1 5 48
3 1 58
2 9 90
 2
23 42
1 2 10
 1
42
 3
1 5 1
1 3 42
1 2 42
 10
1 1 1 1 1 1 1 1 1 1
1 2 5
1 3 5
1 4 5
1 5 5
1 6 5
1 7 5
1 8 5
1 9 5
1 10 5
 3
1 1 1
1 2 42
2 3 42
 3
5 1 1
1 3 42
1 2 24

3 1042.0638260857 1 727.3238618387 6 283959805416.40496826
 3 192.0000000000 1 13952642.4776933547 8 544327.1927847408 2 727.3238618387 1 0.0000000000 2 1042.0638260857 1 100.6230589875 2 544.3822186663 1 389.7666169867
#include <bits/stdc++.h>
using namespace std;
template <class T>
inline void read(T &x) {
  x = 0;
  char c = getchar();
  int f = 1;
  while (!isdigit(c)) {
    if (c == '-') f = -1;
    c = getchar();
  }
  while (isdigit(c)) {
    x = x * 10 + c - '0';
    c = getchar();
  }
  x *= f;
}
template <class T>
inline void umin(T &x, T y) {
  x = x < y ? x : y;
}
template <class T>
inline void umax(T &x, T y) {
  x = x > y ? x : y;
}
inline unsigned int R() {
  static unsigned int seed = 416;
  return seed ^= seed >> 5, seed ^= seed << 17, seed ^= seed >> 13;
}
const int N = 266666;
int n, w[N];
int ans1;
double ans2 = 1e100, cur, ans[N];
bool gg[N];
int sz[N], T;
vector<pair<int, int> > e[N];
void getsz(int u, int fa) {
  sz[u] = 1;
  for (register int i = (0); i <= (((int)e[u].size()) - 1); i++)
    if (!gg[e[u][i].first] && e[u][i].first != fa)
      getsz(e[u][i].first, u), sz[u] += sz[e[u][i].first];
}
void getcen(int u, int fa, pair<int, int> &res) {
  int mx = T - sz[u];
  for (register int i = (0); i <= (((int)e[u].size()) - 1); i++)
    if (!gg[e[u][i].first] && e[u][i].first != fa)
      getcen(e[u][i].first, u, res), umax(mx, sz[e[u][i].first]);
  res = min(res, pair<int, int>(mx, u));
}
void getans(int u, int fa, long long dis, int bel) {
  cur += w[u] * dis * sqrt(dis);
  ans[bel] += w[u] * sqrt(dis);
  for (register int i = (0); i <= (((int)e[u].size()) - 1); i++)
    if (e[u][i].first != fa)
      getans(e[u][i].first, u, dis + e[u][i].second, bel);
}
void solve(int rt) {
  getsz(rt, 0);
  T = sz[rt];
  pair<int, int> res(1e9, rt);
  getcen(rt, 0, res);
  assert(res.first <= T / 2);
  rt = res.second;
  gg[rt] = 1;
  cur = 0;
  double sum = 0;
  cerr << rt << ' ' << T << ' ' << res.first << endl;
  for (register int i = (0); i <= (((int)e[rt].size()) - 1); i++) {
    ans[e[rt][i].first] = 0,
    getans(e[rt][i].first, rt, e[rt][i].second, e[rt][i].first),
    sum += ans[e[rt][i].first];
  }
  if (cur < ans2) ans2 = cur, ans1 = rt;
  for (register int i = (0); i <= (((int)e[rt].size()) - 1); i++)
    if (ans[e[rt][i].first] * 2 > sum && !gg[e[rt][i].first]) {
      solve(e[rt][i].first);
      break;
    }
}
int main() {
  read(n);
  for (register int i = (1); i <= (n); i++) read(w[i]);
  for (register int i = (1); i <= (n - 1); i++) {
    int u, v, w;
    read(u);
    read(v);
    read(w);
    e[u].push_back(pair<int, int>(v, w));
    e[v].push_back(pair<int, int>(u, w));
  }
  solve(1);
  printf("%d %.11lf\n", ans1, ans2);
  return 0;
}
