#include <bits/stdc++.h>
using namespace std;
struct p {
  int l, r, d;
} q[1000010];
;
long long a[1000010];
long long ans;
int sz, n, k;
int cnt[1000100];
long long otv[1000100];
inline void del(long long k) {
  ans = ans - (((cnt[k] * 1ll * cnt[k] * 1ll) -
                ((cnt[k] - 1) * 1ll * (cnt[k] - 1) * 1ll)) *
               1ll * k * 1ll);
  cnt[k]--;
}
inline void add(long long k) {
  ans = ans + ((((cnt[k] + 1) * 1ll * (cnt[k] + 1) * 1ll) -
                ((cnt[k]) * 1ll * (cnt[k]) * 1ll)) *
               1ll * k * 1ll);
  cnt[k]++;
}
bool cmp(p p1, p p2) {
  if ((p1.l / sz) != (p2.l / sz)) return (p1.l / sz) < (p2.l / sz);
  return p1.r < p2.r;
}
int main() {
  scanf("%d%d", &n, &k);
  for (int i = 1; i <= n; i++) scanf("%I64d", &a[i]);
  sz = int(sqrt(n)) + 1;
  for (int i = 1; i <= k; i++) scanf("%d%d", &q[i].l, &q[i].r), q[i].d = i;
  sort(q + 1, q + k + 1, cmp);
  for (int i = q[1].l; i <= q[1].r; i++)
    ans = ans + ((((cnt[a[i]] + 1) * 1ll * (cnt[a[i]] + 1) * 1ll) -
                  ((cnt[a[i]]) * 1ll * (cnt[a[i]]) * 1ll)) *
                 1ll * a[i] * 1ll),
    cnt[a[i]]++;
  int x = q[1].r, y = q[1].l;
  for (int i = 1; i <= k; i++) {
    while (y < q[i].l) del(a[y++]);
    while (x < q[i].r) add(a[++x]);
    while (y > q[i].l) add(a[--y]);
    while (x > q[i].r) del(a[x--]);
    otv[q[i].d] = ans;
  }
  for (int i = 1; i <= k; i++) printf("%I64d\n", otv[i]);
  return 0;
}
