#include <bits/stdc++.h>

using namespace std;
using pii = pair <int, int>;

inline int untitled(const pii &a) {
    if (a.first < 0) {
        // x < 0
        return a.second < 0 ? 3 : 2;
    }else {
        return a.second < 0 ? 4 : 1;
    }
}

inline bool cmp(const pii &a, const pii &b) {
    int a_ = untitled(a), b_ = untitled(b);
    if (a_ != b_) return a_ < b_;
    // dya / dxa < dyb / dxb
    return 1ll * a.second * b.first < 1ll * a.first * b.second;
}

const int MAX_N = 1000 + 5;

int N;
struct Segment{
    int x1, y1, x2, y2;
    Segment (int x1_ = 0, int y1_ = 0, int x2_ = 0, int y2_ = 0) : x1(x1_), y1(y1_), x2(x2_), y2(y2_) {}
    inline bool operator < (const Segment &comp) const {
        return cmp(pair(x2 - x1, y2 - y1), pair(comp.x2 - comp.x1, comp.y2 - comp.y1));
        // return atan2(y2 - y1, x2 - x1) < atan2(comp.y2 - comp.y1, comp.x2 - comp.x1);
    }
}seg[MAX_N];
pii node[MAX_N << 1];

inline bool check(const Segment &seg, int x0, int y0) {
    int dx1 = x0 - seg.x1, dy1 = x0 - seg.y1;
    int dx2 = seg.x2 - x0, dy2 = seg.y2 - y0;
    // dy1 / dx1 = dy2 / dx2
    return 1ll * dy1 * dx2 == 1ll * dx1 * dy2 && min(seg.x1, seg.x2) <= x0 && x0 <= max(seg.x1, seg.x2);
}

void solve() {
    cin >> N;
    for (int i = 1; i <= N; i ++) {
        int x1, y1, x2, y2;
        cin >> x1 >> y1 >> x2 >> y2;
        node[i] = pair(x1, y1);
        node[i + N] = pair(x2, y2);
        seg[i] = Segment(x1, y1, x2, y2);
    }
    int ans = 0;
    sort(node + 1, node + (N << 1) + 1);
    int M = unique(node + 1, node + (N << 1) + 1) - node - 1;
    for (int i = 1; i <= M; i ++) {
        auto [x0, y0] = node[i];
        vector <Segment> seg1;
        vector <pii> seg2;
        for (int j = 1; j <= N; j ++) {
            if (seg[j].x1 == x0 && seg[j].y1 == y0) {
                seg2.emplace_back(seg[j].x2 - seg[j].x1, seg[j].y2 - seg[j].y1);
            }else if (seg[j].x2 == x0 && seg[j].y2 == y0) {
                seg2.emplace_back(seg[j].x1 - seg[j].x2, seg[j].y1 - seg[j].y2);
            }else if (check(seg[j], x0, y0)) {
                seg1.push_back(seg[j]);
            }
        }
        if (!(seg1.empty() || seg2.empty())) {
            sort(seg1.begin(), seg1.end());
            sort(seg2.begin(), seg2.end(), cmp);
            auto ed = unique(seg2.begin(), seg2.end());
            int n = ed - seg2.begin();
            int i = 0, j = n - 1;
            // fprintf(stderr, "n = %d\n", n);
            // for (int i = 0; i < n; i ++) fprintf(stderr, "i = %d, dx = %d, dy = %d\n", i, seg2[i].first, seg2[i].second);
            for (int k = 0; k < (int)seg1.size(); k ++) {
                pii dir0 = pair(seg1[k].x2 - seg1[k].x1, seg1[k].y2 - seg1[k].y1);
                pii dir1 = pair(seg1[k].x1 - seg1[k].x2, seg1[k].y1 - seg1[k].y2);
                // fprintf(stderr, "k = %d, dir0 = (%d, %d), dir1 = (%d, %d)\n", k, dir0.first, dir0.second, dir1.first, dir1.second);
                if (cmp(dir0, seg2[i])) i = 0;
                if (cmp(seg2[j], dir1)) j = n - 1;
                // i < dir0, j > dir1
                while (i < n && cmp(seg2[i], dir0)) i ++;
                i --; if (i < 0) i = n - 1;
                while (j >= 0 && cmp(dir1, seg2[j])) j --;
                j ++; if (j == n) j = 0;
                int cnt0 = j <= i ? i - j + 1 : n - (j - i - 1), cnt1 = n - cnt0;
                ans += (cnt0 * (cnt0 - 1) + cnt1 * (cnt1 - 1)) / 2;
            }
        }
    }
    cout << ans << '\n';
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(0); cout.tie(0);
    
    int T; cin >> T;
    while (T --) solve();
    
    return 0;
}