#include <bits/stdc++.h>
using namespace std;
#define int long long
#pragma GCC optimize("O3,unroll-loops")

#define rep(i, a, b) for(int i = a; i < (b); ++i)
#define all(x) begin(x), end(x)
#define sz(x) (int)(x).size()
typedef long long ll;
typedef pair<int, int> pii;
typedef vector<int> vi;
template <class T> int sgn(T x) { return (x > 0) - (x < 0); }
template<class T>
struct Point {
	typedef Point P;
	T x, y;
	explicit Point(T x=0, T y=0) : x(x), y(y) {}
	bool operator<(P p) const { return tie(x,y) < tie(p.x,p.y); }
	bool operator==(P p) const { return tie(x,y)==tie(p.x,p.y); }
	P operator+(P p) const { return P(x+p.x, y+p.y); }
	P operator-(P p) const { return P(x-p.x, y-p.y); }
	P operator*(T d) const { return P(x*d, y*d); }
	P operator/(T d) const { return P(x/d, y/d); }
	T dot(P p) const { return x*p.x + y*p.y; }
	T cross(P p) const { return x*p.y - y*p.x; }
	T cross(P a, P b) const { return (a-*this).cross(b-*this); }
	T dist2() const { return x*x + y*y; }
	double dist() const { return sqrt((double)dist2()); }
	// angle to x-axis in interval [-pi, pi]
	double angle() const { return atan2(y, x); }
	P unit() const { return *this/dist(); } // makes dist()=1
	P perp() const { return P(-y, x); } // rotates +90 degrees
	P normal() const { return perp().unit(); }
	// returns point rotated 'a' radians ccw around the origin
	P rotate(double a) const {
		return P(x*cos(a)-y*sin(a),x*sin(a)+y*cos(a)); }
	friend ostream& operator<<(ostream& os, P p) {
		return os << "(" << p.x << "," << p.y << ")"; }
};

typedef Point<int> P;

P fix(P a){
	if(a.y < 0 || (a.y == 0 && a.x < 0)){
		a = a*-1;
	}
	return a;
}

template<class T>
T polygonArea2(vector<Point<T>>& v) {
	T a = v.back().cross(v[0]);
	rep(i,0,sz(v)-1) a += v[i].cross(v[i+1]);
	return a;
}

template<class P> bool onSegment(P s, P e, P p) {
	return p.cross(s, e) == 0 && (s - p).dot(e - p) <= 0;
}

string f(int a){
	string s = to_string(a/2);
	s.push_back('.');
	if(a%2==1)s.push_back('5');
	else s.push_back('0');
	return s;
}
int pos(int x, int y){
	if (y < 0) return -1;
        if (y == 0 && 0 <= x) return 0;
        return 1;
}

signed main(){
	cin.tie(nullptr)->sync_with_stdio(false);
	int t;
	cin >> t;
	while(t--){
		int n;
		cin >> n;
		if(n==0)return 0;
		map<P,int>hm;
		int m = 0;
		vector<pair<int,int>>a(n);
		vector<P>arr;
		
		for(int i = 0; i<n; i++){
			int x,y;
			cin >> x >> y;
			P p{x,y};
			cin >> x >> y;
			P p2{x,y};
			int u = 0; int v = 0;
			if(!hm.count(p)){
				hm[p] = m;
				m++;
				arr.push_back(p);
			}
			if(!hm.count(p2)){
				hm[p2] = m;
				m++;
				arr.push_back(p2);
			}
			u = hm[p];
			v = hm[p2];
			a[i] = {u,v};
		}
		vector<int>have(m*m);
		assert(arr.size()==m);
		vector<vector<int>>adj(m);
		vector<int>hid(m);
		for(int i = 0; i<n; i++){
			int u = a[i].first;
			int v = a[i].second;
			adj[u].push_back(v);
			adj[v].push_back(u);
			//cout << u << ' ' << v << '\n';
		}
		for(int i = 0; i<m; i++){
			sort(adj[i].begin(),adj[i].end(),[&](int a, int b){
				int ax = arr[a].x;
				int ay = arr[a].y;
				int bx = arr[b].x;
				int by = arr[b].y;
				if(pos(ax,ay) != pos(bx,by))return pos(ax,ay) < pos(bx,by);
				return 0 < (ax * by - ay * bx);
			});
		}
		sort(arr.begin(),arr.end(), [](P a, P b){
			return P{a.x,-a.y} < P{b.x,-b.y};
		});
		vector<pair<int,int>>vals;
		vector<int>idx(m);
		for(int i = 0; i<m; i++){
			hid[hm[arr[i]]] = i;
		}
		for(int i = 0; i<m; i++){
			idx[i] = i;
		}
		for(int i = 0; i<n; i++){
			int u = hid[a[i].first];
			int v = hid[a[i].second];
			have[u*m+v]++;
			have[v*m+u]++;
		}
		for(int i = 0; i<m; i++){
			for(int j = i+1; j<m; j++){
				vals.push_back({i,j});
			}
		}
		sort(vals.begin(),vals.end(), [&](pair<int,int>a, pair<int,int> b){
			P vec1 = (arr[a.first]-arr[a.second]).perp();
			P vec2 = (arr[b.first]-arr[b.second]).perp();
			vec1 = fix(vec1);
			vec2 = fix(vec2);
			if(vec1.cross(vec2)!=0)
				return vec1.cross(vec2) > 0;
			return make_pair(a.second,a.first) < make_pair(b.second,b.first);
		});
		int ans = 0;
		for(auto [c,d] : vals){
			int x = idx[c];
			int y = idx[d];
			
			int l = min(x,y);
			int r = max(x,y);
			//cout << l << ' ' << r << '\n';
			/*
			cout << c << ' ' << d << '\n';
			cout << arr[l] << ' ' << arr[r] << '\n';
			for(int i = 0; i<m; i++){
				cout << arr[i] << ' ';
			}
			cout << '\n';
			*/
			if(have[c*m+d]){
				int mul = have[c*m+d];
				vector<bool>bad(m);
				for(int i = 0; i<m; i++){
					if(i==l || i==r)bad[i] = true;
					else if(arr[i].cross(arr[l],arr[r]) == 0){
						bad[i] = true;
					}
				}
				for(int i = 0; i<m; i++){
					if(i==l || i==r)continue;
					//
					if(onSegment(arr[l],arr[r],arr[i])){
						//cout << arr[l] << ' ' << arr[r] << ' ' << arr[i] << '\n';
						int u = hm[arr[i]];
						int cl = 0;
						int cr = 0;
						for(int nxt: adj[u]){
							int id = idx[hid[nxt]];
							if(bad[id])continue;
							if(id < l){
								cl++;
								//lt.push_back(arr[id]-arr[i]);
							}
							else{
								if(id < r){
									//cout << id << ' ' << l << ' ' << r << '\n';
									//cout << arr[id] << ' ' << arr[l] << ' ' << arr[r] << '\n';
								}
								assert(id > r);
								cr++;
								//rt.push_back(arr[id]-arr[i]);
							}
						}
						int cur = 0;
						for(int j = 1; j<adj[u].size(); j++){
							int id1 = idx[hid[adj[u][j-1]]];
							int id2 = idx[hid[adj[u][j]]];
							if(bad[id2] || bad[id1]){
								cur = 0;
								continue;
							}
							if((id1 < l && id2 > r) || (id1 > r && id2 < l)){
								cur = 0;
								continue;
							}
							if(arr[i].cross(arr[id1],arr[id2]) == 0){
								cur++;
								ans-=cur;
								//cout << arr[i] << ' ' << arr[id1] << ' ' << arr[id2] << '\n';
							}
							else{
								cur = 0;
							}
						}
						//cout << cl << ' ' << cr << '\n';
						ans+=cl*(cl-1)/2 * mul;
						ans+=cr*(cr-1)/2 * mul;
					}
				}
				//cout << arr[l] << ' ' << arr[r] << '\n';
			}
			
			swap(idx[c],idx[d]);
			swap(arr[x],arr[y]);
		}
		cout << ans << '\n';
	}
	return 0;
}