//Krzysztof Boryczka
#pragma GCC optimize "O3"
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef long double ld;
typedef pair<int, int> ii;
typedef vector<int> vi;
typedef vector<ii> vii;
const int inf=0x3f3f3f3f;
const ll INF=0x3f3f3f3f3f3f3f3f;
#define FOR(i, b, e) for(int i=b; i<=e; i++)
#define FORD(i, b, e) for(int i=b; i>=e; i--)
#define SIZE(x) ((int)x.size())
#define pb push_back
#define st first
#define nd second
#define sp ' '
#define ent '\n'
const int N=1e5+5;
const ld EPS=1e-9;
int n;
pair<ld, ld> mam[N], dam[N];
ii m2[N], d2[N];
pair<ld, ld> akt;
void dolej(pair<ld, ld> x){
if(akt.nd < EPS) akt={0.0, 0.0};
if(x.nd < EPS) return;
else akt={(akt.st*akt.nd+x.st*x.nd)/(akt.nd+x.nd), akt.nd+x.nd};
}
void solve(){
int a, b, c;
cin>>n;
FOR(i, 1, n){
cin>>a>>b>>c;
m2[i]={b, a};
d2[i]={c, a};
}
sort(m2+1, m2+n+1);
sort(d2+1, d2+n+1);
FOR(i, 1, n) mam[i]=m2[i];
FOR(i, 1, n) dam[i]=d2[i];
akt={0.0, 0.0};
int pt=0;
FOR(i, 1, n){
while(pt<n && dam[i].nd > EPS){
pt++;
if(abs(mam[pt].st-dam[i].st) < EPS){
ld mink=min(mam[pt].nd, dam[i].nd);
mam[pt].nd-=mink;
dam[i].nd-=mink;
dolej(mam[pt]);
continue;
}
if(mam[pt].st < dam[i].st){
dolej(mam[pt]);
continue;
}
if(akt.nd < EPS){
break;
}
if(abs(dam[i].st-akt.st) < EPS){
ld mink=min(akt.nd, dam[i].nd);
akt.nd-=mink;
dam[i].nd-=mink;
if(akt.nd < EPS) akt={0.0, 0.0};
pt--;
continue;
}
ld prop=(dam[i].st-akt.st)/(mam[pt].st-dam[i].st);
if(prop > EPS){
ld leje=min({akt.nd, mam[pt].nd/prop, dam[i].nd/(prop+1.0)});
dam[i].nd-=leje*(prop+1.0);
akt.nd-=leje;
mam[pt].nd-=leje*prop;
}
if(mam[pt].nd > EPS) pt--;
}
if(dam[i].nd > EPS){
cout<<"NIE\n";
return;
}
}
cout<<"TAK\n";
}
int main(){
ios::sync_with_stdio(0);
cin.tie(0); cout.tie(0);
int tt; cin>>tt;
FOR(te, 1, tt)
solve();
return 0;
}
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 | //Krzysztof Boryczka #pragma GCC optimize "O3" #include <bits/stdc++.h> using namespace std; typedef long long ll; typedef long double ld; typedef pair<int, int> ii; typedef vector<int> vi; typedef vector<ii> vii; const int inf=0x3f3f3f3f; const ll INF=0x3f3f3f3f3f3f3f3f; #define FOR(i, b, e) for(int i=b; i<=e; i++) #define FORD(i, b, e) for(int i=b; i>=e; i--) #define SIZE(x) ((int)x.size()) #define pb push_back #define st first #define nd second #define sp ' ' #define ent '\n' const int N=1e5+5; const ld EPS=1e-9; int n; pair<ld, ld> mam[N], dam[N]; ii m2[N], d2[N]; pair<ld, ld> akt; void dolej(pair<ld, ld> x){ if(akt.nd < EPS) akt={0.0, 0.0}; if(x.nd < EPS) return; else akt={(akt.st*akt.nd+x.st*x.nd)/(akt.nd+x.nd), akt.nd+x.nd}; } void solve(){ int a, b, c; cin>>n; FOR(i, 1, n){ cin>>a>>b>>c; m2[i]={b, a}; d2[i]={c, a}; } sort(m2+1, m2+n+1); sort(d2+1, d2+n+1); FOR(i, 1, n) mam[i]=m2[i]; FOR(i, 1, n) dam[i]=d2[i]; akt={0.0, 0.0}; int pt=0; FOR(i, 1, n){ while(pt<n && dam[i].nd > EPS){ pt++; if(abs(mam[pt].st-dam[i].st) < EPS){ ld mink=min(mam[pt].nd, dam[i].nd); mam[pt].nd-=mink; dam[i].nd-=mink; dolej(mam[pt]); continue; } if(mam[pt].st < dam[i].st){ dolej(mam[pt]); continue; } if(akt.nd < EPS){ break; } if(abs(dam[i].st-akt.st) < EPS){ ld mink=min(akt.nd, dam[i].nd); akt.nd-=mink; dam[i].nd-=mink; if(akt.nd < EPS) akt={0.0, 0.0}; pt--; continue; } ld prop=(dam[i].st-akt.st)/(mam[pt].st-dam[i].st); if(prop > EPS){ ld leje=min({akt.nd, mam[pt].nd/prop, dam[i].nd/(prop+1.0)}); dam[i].nd-=leje*(prop+1.0); akt.nd-=leje; mam[pt].nd-=leje*prop; } if(mam[pt].nd > EPS) pt--; } if(dam[i].nd > EPS){ cout<<"NIE\n"; return; } } cout<<"TAK\n"; } int main(){ ios::sync_with_stdio(0); cin.tie(0); cout.tie(0); int tt; cin>>tt; FOR(te, 1, tt) solve(); return 0; } |
English