#include<bits/stdc++.h>
#define VAR(i,n) __typeof(n) i = (n)
#define loop(i,j,s) for(int i=j;i<s;i++)
#define loopback(i,j,s) for(int i=j;i>=s;i--)
#define foreach(i,c) for(VAR(i,(c).begin());i!=(c).end();i++)
#define pln( x ) cout << x << "\n"
#define ps( x ) cout << x << " "
#define entr cout << "\n"
#define pcnt(i) __builtin_popcount(i)
#define ll long long
#define pb push_back
#define mp make_pair
#define ff first
#define ss second
#define SIZE(c) (c).size()
#define ALL(c) (c).begin(), (c).end()
using namespace std;
typedef vector<int> VI;
typedef pair<int, int> PII;
typedef vector<vector<int> > VVI;
const ll INFTY=2000000000000ll;
const int MAX=500100;
const int MOD=10000000;
void coutTab(int* tab,int n){
loop(i,0,n){
cout<<tab[i]<<" ";
}
cout<<"\n";
}
template<class T> void coutVec(vector<T> tab){
for(T t : tab){
cout<<t<<" ";
}
cout<<"\n";
}
//------------------------------------------
ll l[MAX];
ll a[MAX];
ll b[MAX];
int n;
const double EPS = 1e-6;
int comp(double a, double b) {
if(abs(a-b) < EPS) return 0;
return (a<b) ? -1 : 1;
}
set<int>::iterator last_lower_eq_than(set<int> & A, int a) {
auto it = A.upper_bound(a);
///ps("tryue ");pln((it==A.begin()));
if(it == A.begin()) return A.end();
it--;
return it;
}
bool check() {
map<int, double> cA;
map<int, double> cB;
set<int> A;
set<int> B;
loop(i,0,n) {
cA[a[i]]+=l[i];
A.insert(a[i]);
cB[b[i]]+=l[i];
B.insert(b[i]);
}
while(!A.empty() || !B.empty()) {
int et = *B.rbegin();
int t2 = *A.rbegin();
set<int>::iterator it1;
if((it1 = last_lower_eq_than(A, et)) == A.end()) return false;
int t1 = *it1;
if(t2 < et) return false;
//ps(t1);ps(t2);pln(et);
//ps("test ");
//for(auto i : A) ps(i);
//entr;
//for(auto i : B) ps(i);
//entr;
//entr;
if(t1 == et) {
double r = min(cA[t1], cB[et]);
cA[t1] -= r;
cB[et] -= r;
} else if(t2 == et) {
double r = min(cA[t2], cB[et]);
cA[t2] -= r;
cB[et] -= r;
} else {
double ax = cB[et]*(et - t2)/(t1-t2);
ax = min(ax, cA[t1]);
double bx = ax*(t1-et)/(et-t2);
bx = min(bx, cA[t2]);
ax = bx*(et-t2)/(t1-et);
cA[t1] -= ax;
cA[t2] -= bx;
cB[et] -= (ax+bx);
}
if(comp(cA[t1],0) == 0) {
A.erase(t1);
}
if(comp(cA[t2],0) == 0) {
A.erase(t2);
}
if(comp(cB[et],0) == 0) {
B.erase(et);
}
}
if(!A.empty() || !B.empty()) return false;
return true;
}
int main(){
ios_base::sync_with_stdio(0);
int t;
cin>>t;
while(t--) {
cin>>n;
loop(i,0,n) cin>>l[i]>>a[i]>>b[i];
if(check()) pln("TAK");
else pln("NIE");
}
}
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 | #include<bits/stdc++.h> #define VAR(i,n) __typeof(n) i = (n) #define loop(i,j,s) for(int i=j;i<s;i++) #define loopback(i,j,s) for(int i=j;i>=s;i--) #define foreach(i,c) for(VAR(i,(c).begin());i!=(c).end();i++) #define pln( x ) cout << x << "\n" #define ps( x ) cout << x << " " #define entr cout << "\n" #define pcnt(i) __builtin_popcount(i) #define ll long long #define pb push_back #define mp make_pair #define ff first #define ss second #define SIZE(c) (c).size() #define ALL(c) (c).begin(), (c).end() using namespace std; typedef vector<int> VI; typedef pair<int, int> PII; typedef vector<vector<int> > VVI; const ll INFTY=2000000000000ll; const int MAX=500100; const int MOD=10000000; void coutTab(int* tab,int n){ loop(i,0,n){ cout<<tab[i]<<" "; } cout<<"\n"; } template<class T> void coutVec(vector<T> tab){ for(T t : tab){ cout<<t<<" "; } cout<<"\n"; } //------------------------------------------ ll l[MAX]; ll a[MAX]; ll b[MAX]; int n; const double EPS = 1e-6; int comp(double a, double b) { if(abs(a-b) < EPS) return 0; return (a<b) ? -1 : 1; } set<int>::iterator last_lower_eq_than(set<int> & A, int a) { auto it = A.upper_bound(a); ///ps("tryue ");pln((it==A.begin())); if(it == A.begin()) return A.end(); it--; return it; } bool check() { map<int, double> cA; map<int, double> cB; set<int> A; set<int> B; loop(i,0,n) { cA[a[i]]+=l[i]; A.insert(a[i]); cB[b[i]]+=l[i]; B.insert(b[i]); } while(!A.empty() || !B.empty()) { int et = *B.rbegin(); int t2 = *A.rbegin(); set<int>::iterator it1; if((it1 = last_lower_eq_than(A, et)) == A.end()) return false; int t1 = *it1; if(t2 < et) return false; //ps(t1);ps(t2);pln(et); //ps("test "); //for(auto i : A) ps(i); //entr; //for(auto i : B) ps(i); //entr; //entr; if(t1 == et) { double r = min(cA[t1], cB[et]); cA[t1] -= r; cB[et] -= r; } else if(t2 == et) { double r = min(cA[t2], cB[et]); cA[t2] -= r; cB[et] -= r; } else { double ax = cB[et]*(et - t2)/(t1-t2); ax = min(ax, cA[t1]); double bx = ax*(t1-et)/(et-t2); bx = min(bx, cA[t2]); ax = bx*(et-t2)/(t1-et); cA[t1] -= ax; cA[t2] -= bx; cB[et] -= (ax+bx); } if(comp(cA[t1],0) == 0) { A.erase(t1); } if(comp(cA[t2],0) == 0) { A.erase(t2); } if(comp(cB[et],0) == 0) { B.erase(et); } } if(!A.empty() || !B.empty()) return false; return true; } int main(){ ios_base::sync_with_stdio(0); int t; cin>>t; while(t--) { cin>>n; loop(i,0,n) cin>>l[i]>>a[i]>>b[i]; if(check()) pln("TAK"); else pln("NIE"); } } |
English