#include <cstdio>
#include <vector>
#include <algorithm>
using namespace std;
long double e = (long double)1e-9;
int main(){
int t;
scanf("%d", &t);
while (t--){
int n;
scanf("%d", &n);
std::vector<std::pair<long long, long double>> a, b;
long long e1 = 0, e2 = 0;
for (int i = 0; i < n; i++){
long long _l, _a, _b;
scanf("%lld%lld%lld", &_l, &_a, &_b);
a.push_back(make_pair(_a, (long double)_l));
b.push_back(make_pair(_b, (long double)_l));
e1 += _l * _a;
e2 += _l * _b;
}
if (e1 != e2){
printf("NIE\n");
continue;
}
sort(a.begin(), a.end());
sort(b.begin(), b.end());
/*for (int i = 0; i < n; i++){
printf("(%lld %lf) ", b[i].first, (double)b[i].second);
}
printf("\n");
for (int i = 0; i < n; i++){
printf("(%lld %lf) ", a[i].first, (double)a[i].second);
}
printf("\n");*/
int p = 0, q = 0;
bool ok = true;
for (int i = 0; i < (int)b.size(); i++){
//printf("cc %d %d %lld %lf %d %lld %lf\n", i, p, a[p].first, (double)a[p].second, i, b[i].first, (double)b[i].second);
while ((q < (int)a.size()) && (a[q].first < b[i].first)){
q++;
}
//printf("cd %d\n", q);
while (q < (int)a.size()){
if (a[q].first == b[i].first){
if (a[q].second >= b[i].second){
a[q].second -= b[i].second;
b[i].second = 0;
break;
}
else {
b[i].second -= a[q].second;
a[q].second = 0;
q++;
continue;
}
}
else if ((a[p].first >= b[i].first) && (b[i].second > e)){
//printf("h1 %d %lld %lf %d %lld %lf\n", p, a[p].first, (double)a[p].second, i, b[i].first, (double)b[i].second);
ok = false;
break;
}
else {
long double ratio = (long double)(b[i].first - a[p].first) / (a[q].first - a[p].first);
if ((a[p].second <= b[i].second * (1 - ratio)) && (a[p].second * ratio <= a[q].second * (1 - ratio))){
b[i].second -= a[p].second / (1 - ratio);
a[q].second -= a[p].second * ratio / (1 - ratio);
a[p].second = 0;
}
else if ((a[q].second <= b[i].second * ratio) && (a[q].second * (1 - ratio) <= a[p].second * ratio)){
b[i].second -= a[q].second / ratio;
a[p].second -= a[q].second * (1 - ratio) / ratio;
a[q].second = 0;
}
else {
a[p].second -= b[i].second * (1 - ratio);
a[q].second -= b[i].second * ratio;
b[i].second = 0;
}
if (a[p].second < e){
p++;
}
if (a[q].second < e){
//printf("wro %d %d %d %lf %lld %lld %lld %lf %lf %lf\n", q, p, i, (double)ratio, a[p].first, b[i].first, a[q].first, (double)a[p].second, (double)b[i].second, (double)a[q].second);
q++;
}
if (b[i].second < e){
break;
}
}
}
if ((!ok) || ((q >= (int)a.size()) && (!(b[i].second < e)))){
//if (t == 100 - 57)
//printf("h2 %d %d %d %d %d %.20lf\n", i, n, q, (int)a.size(), (int)ok, (double)b[i].second);
ok = false;
break;
}
}
if (!ok){
printf("NIE\n");
}
else {
printf("TAK\n");
}
}
}
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 | #include <cstdio> #include <vector> #include <algorithm> using namespace std; long double e = (long double)1e-9; int main(){ int t; scanf("%d", &t); while (t--){ int n; scanf("%d", &n); std::vector<std::pair<long long, long double>> a, b; long long e1 = 0, e2 = 0; for (int i = 0; i < n; i++){ long long _l, _a, _b; scanf("%lld%lld%lld", &_l, &_a, &_b); a.push_back(make_pair(_a, (long double)_l)); b.push_back(make_pair(_b, (long double)_l)); e1 += _l * _a; e2 += _l * _b; } if (e1 != e2){ printf("NIE\n"); continue; } sort(a.begin(), a.end()); sort(b.begin(), b.end()); /*for (int i = 0; i < n; i++){ printf("(%lld %lf) ", b[i].first, (double)b[i].second); } printf("\n"); for (int i = 0; i < n; i++){ printf("(%lld %lf) ", a[i].first, (double)a[i].second); } printf("\n");*/ int p = 0, q = 0; bool ok = true; for (int i = 0; i < (int)b.size(); i++){ //printf("cc %d %d %lld %lf %d %lld %lf\n", i, p, a[p].first, (double)a[p].second, i, b[i].first, (double)b[i].second); while ((q < (int)a.size()) && (a[q].first < b[i].first)){ q++; } //printf("cd %d\n", q); while (q < (int)a.size()){ if (a[q].first == b[i].first){ if (a[q].second >= b[i].second){ a[q].second -= b[i].second; b[i].second = 0; break; } else { b[i].second -= a[q].second; a[q].second = 0; q++; continue; } } else if ((a[p].first >= b[i].first) && (b[i].second > e)){ //printf("h1 %d %lld %lf %d %lld %lf\n", p, a[p].first, (double)a[p].second, i, b[i].first, (double)b[i].second); ok = false; break; } else { long double ratio = (long double)(b[i].first - a[p].first) / (a[q].first - a[p].first); if ((a[p].second <= b[i].second * (1 - ratio)) && (a[p].second * ratio <= a[q].second * (1 - ratio))){ b[i].second -= a[p].second / (1 - ratio); a[q].second -= a[p].second * ratio / (1 - ratio); a[p].second = 0; } else if ((a[q].second <= b[i].second * ratio) && (a[q].second * (1 - ratio) <= a[p].second * ratio)){ b[i].second -= a[q].second / ratio; a[p].second -= a[q].second * (1 - ratio) / ratio; a[q].second = 0; } else { a[p].second -= b[i].second * (1 - ratio); a[q].second -= b[i].second * ratio; b[i].second = 0; } if (a[p].second < e){ p++; } if (a[q].second < e){ //printf("wro %d %d %d %lf %lld %lld %lld %lf %lf %lf\n", q, p, i, (double)ratio, a[p].first, b[i].first, a[q].first, (double)a[p].second, (double)b[i].second, (double)a[q].second); q++; } if (b[i].second < e){ break; } } } if ((!ok) || ((q >= (int)a.size()) && (!(b[i].second < e)))){ //if (t == 100 - 57) //printf("h2 %d %d %d %d %d %.20lf\n", i, n, q, (int)a.size(), (int)ok, (double)b[i].second); ok = false; break; } } if (!ok){ printf("NIE\n"); } else { printf("TAK\n"); } } } |
English