Make Them Equal
You are given an array $$$a$$$ consisting of $$$n$$$ positive integers, numbered from $$$1$$$ to $$$n$$$. You can perform the following operation no more than $$$3n$$$ times:  choose three integers $$$i$$$, $$$j$$$ and $$$x$$$ ($$$1 \le i, j \le n$$$; $$$0 \le x \le 10^9$$$);  assign $$$a_i := a_i - x \cdot i$$$, $$$a_j := a_j + x \cdot i$$$. After each operation, all elements of the array should be non-negative.Can you find a sequence of no more than $$$3n$$$ operations after which all elements of the array are equal?
The first line contains one integer $$$t$$$ ($$$1 \le t \le 10^4$$$) — the number of test cases. Then $$$t$$$ test cases follow.The first line of each test case contains one integer $$$n$$$ ($$$1 \le n \le 10^4$$$) — the number of elements in the array. The second line contains $$$n$$$ integers $$$a_1, a_2, \dots, a_n$$$ ($$$1 \le a_i \le 10^5$$$) — the elements of the array.It is guaranteed that the sum of $$$n$$$ over all test cases does not exceed $$$10^4$$$.
For each test case print the answer to it as follows:  if there is no suitable sequence of operations, print $$$-1$$$;  otherwise, print one integer $$$k$$$ ($$$0 \le k \le 3n$$$) — the number of operations in the sequence. Then print $$$k$$$ lines, the $$$m$$$-th of which should contain three integers $$$i$$$, $$$j$$$ and $$$x$$$ ($$$1 \le i, j \le n$$$; $$$0 \le x \le 10^9$$$) for the $$$m$$$-th operation. If there are multiple suitable sequences of operations, print any of them. Note that you don't have to minimize $$$k$$$.
import math
import sys
from typing import List
import heapq

def main():
    t = int(input())
    for _ in range(t):
        n = int(input())
        a = [int(x) for x in input().split()]
        s = sum(a)
        if s % n != 0:
            print(-1)
            continue
        s //= n
        print(3 * (n - 1))
        for i in range(n)[1:]:
            print(1, i + 1, (i + 1 - a[i] % (i + 1)) % (i + 1))
            a[i] += (i + 1 - a[i] % (i + 1)) % (i + 1)
            print(i + 1, 1, a[i] // (i + 1))
        for i in range(n)[1:]:
            print(1, i + 1, s)
if __name__ == '__main__':
    main()
for _ in range(int(input())):
    n = int(input())
    a = list(map(int, input().split()))
    s = sum(a)
    if s % n != 0:
        print(-1)
        continue
    t = s // n
    print(3 * (n - 1))
    for i in range(1, n):
        x = i + 1 - a[i] % (i + 1)
        if a[i] % (i + 1) == 0:
            x = 0
        a[i] += x
        print(1, i + 1, x)
        print(i + 1, 1, a[i] // (i + 1))
    for i in range(1, n):
        print(1, i + 1, t)
T = int(input())
for _ in range(T):
    n = int(input())
    y = list(map(int, input().split()))
    s = sum(y)
    if s % n == 0:
        res = list()
        for k in range(1, n):
            if y[k] % (k + 1) != 0:
                r = y[k] % (k + 1)
                x = k + 1 - r
                y[k] += x
                y[0] -= x
                res.append([1, k + 1, x])
            x = int(y[k] / (k + 1))
            y[k] -= x * (k + 1)
            y[0] += x * (k + 1)
            res.append([k + 1, 1, x])
        avg = int(s / n)
        for k in range(1, n):
            y[0] -= avg
            y[k] += avg
            res.append([1, k + 1, avg])
        print(len(res))
        for i in res:
            print(i[0], i[1], i[2])
    else:
        print(-1)