QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#665150 | #5604. Triangle Containment | enze114514 | WA | 318ms | 41952kb | C++20 | 9.5kb | 2024-10-22 06:48:38 | 2024-10-22 06:48:38 |
Judging History
answer
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
using ld = long double;
#define pb push_back
const ld pi = 3.14159265358979323846;
const ll INF = 1e18;
const int mod = (int)1e9 + 7;
template<typename T>
T chmax(T a, T b) {
return a > b ? a : b;
}
template<typename T>
T chmin(T a, T b) {
return a > b ? b : a;
}
template<typename T>
struct Point{
T x, y;
ld eps;
Point() : x(0), y(0), eps(1e-9) {}
Point(T x, T y) : x(x), y(y), eps(1e-9) {}
void set_eps(T eps){
this->eps = eps;
}
Point operator+ (const Point& b){
return Point(x + b.x, y + b.y);
}
Point operator- (const Point& b){
return Point(x - b.x, y - b.y);
}
Point operator- (){
return Point(-x, -y);
}
Point operator* (T t) const{
return Point(x * t, y * t);
}
Point operator/ (T t) const{
return Point(x / t, y / t);
}
Point &operator+=(Point p) &{
x += p.x;
y += p.y;
return *this;
}
Point &operator-=(Point p) &{
x -= p.x;
y -= p.y;
return *this;
}
Point &operator*=(T v) &{
x *= v;
y *= v;
return *this;
}
Point &operator/=(T v) &{
x /= v;
y /= v;
return *this;
}
Point &operator=(const Point& b) &{
x = b.x;
y = b.y;
return *this;
}
friend Point operator+ (const Point& a, const Point& b){
return {a.x + b.x, a.y + b.y};
}
friend Point operator- (const Point& a, const Point& b){
return {a.x - b.x, a.y - b.y};
}
friend bool operator==(Point a, Point b){
return a.x == b.x && a.y == b.y;
}
int sign(T x){
if(fabs(x) < eps){
return 0;
}
if(x < 0){
return -1;
}
return 1;
}
bool cmp(T x, T y){
if(fabs(x - y) > eps){
return 0;
}
return 1;
}
bool cmp(const Point& a, const Point& b){
return cmp(a.x, b.x) && cmp(a.y, b.y);
}
T dot(const Point& a, const Point& b){
return a.x * b.x + a.y * b.y;
}
T square(Point a){
return dot(a, a);
}
T cross(const Point& a, const Point& b){
return a.x * b.y - a.y * b.x;
}
T cross(const Point& a, const Point& b, const Point& p){
return (b.x - a.x) * (p.y - a.y) - (b.y - a.y) * (p.x - a.x);
}
T get_len(const Point& a){
return sqrt(dot(a, a));
}
T get_angle(const Point& a, const Point& b){
return acos(dot(a, b) / get_len(a) / get_len(b));
}
T area(const Point& a, const Point& b, const Point& c){
return cross(b - a, c - a);
}
Point rotate(const Point& a, T angle){ //两个点就 b - a (b按a转)
T dx = a.x * cos(angle) + a.y * sin(angle);
T dy = -a.x * sin(angle) + a.y * cos(angle);
return Point(dx, dy);
}
Point intersect(const Point& p, const Point& v, const Point& q, const Point& w){
Point u = p - q;
T t = cross(w, u) / cross(v, w);
return p + v * t;
}
T point_dist(const Point& a, const Point& b){
return sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y));
}
T line_dist(const Point& a, const Point& b, const Point& p){
Point u = b - a, v = p - a;
return fabs(cross(u, v)) / get_len(u);
}
T get_slope(const Point& a, const Point& b){
if(b.y == a.y) return INF;
if(b.x == a.x) return 0;
return (b.y - a.y) / (b.x - a.x);
}
T circle_intersect(const Point& p1, const Point& p2, const T r1, const T r2){
ld d = sqrt((p2.x - p1.x) * (p2.x - p1.x) + (p2.y - p1.y) * (p2.y - p1.y));
if(d > r1 + r2 || d + chmin(r1, r2) <= chmax(r1, r2)){
return 0;
}
else if(d == r1 + r2){
return 1;
}
else{
return 2;
}
}
T seg_dist(const Point& a, const Point& b, const Point& p){
if(a == b){
return get_len(p - a);
}
Point u = b - a, v = p - a, w = p - b;
if(sign(dot(u, v)) < 0){
return get_len(v);
}
if(sign(dot(u, w)) > 0){
return get_len(w);
}
return line_dist(a, b, p);
}
Point projection(const Point& a, const Point& b, const Point& p){
Point v = b - a;
return a + v * dot(v, p - a) / get_len(v);
}
bool on_segment(const Point& a, const Point& b, const Point& p){
bool u = sign(cross(p - a, p - b)) == 0;
bool v = sign(dot(p - a, p - b)) <= 0;
return u && v;
}
bool seg_intersection(const Point& a1, const Point& a2, const Point& b1, const Point& b2){
T c1 = cross(a2 - a1, b1 - a1), c2 = cross(a2 - a1, b2 - a1);
T c3 = cross(b2 - b1, a2 - b1), c4 = cross(b2 - b1, a1 - b1);
return sign(c1) * sign(c2) <= 0 && sign(c3) * sign(c4) <= 0;
}
friend std::ostream &operator<<(ostream &os, Point p) {
return os << "(" << p.x << ", " << p.y << ")";
}
};
bool cmp(vector<ld> a, vector<ld> b){
return a[0] == b[0] ? a[1] < b[1] : a[0] > b[0];
}
template<class Operation, class Mark>
struct SegTree{
const int n;
vector<Operation> op;
vector<Mark> mrk;
SegTree(int n) : n(n), op(4 << __lg(n)), mrk(4 << __lg(n)){
function<void(int, int, int)> build = [&](int u, int l, int r){
op[u] = Operation();
if(l == r) return;
int mid = l + r >> 1;
build(u << 1, l, mid);
build(u << 1 | 1, mid + 1, r);
pushup(u);
};
build(1, 1, n);
}
void pushup(int u){
op[u] = op[u << 1] + op[u << 1 | 1];
}
void modify(int u, const Mark &mk){
op[u].modify(mk);
mrk[u].modify(mk);
}
void pushdown(int u) {
modify(u << 1, mrk[u]);
modify(u << 1 | 1, mrk[u]);
mrk[u] = Mark();
}
void update(int u, int l, int r, int x, const Operation &v){
if(l == r){
op[u] = v;
return;
}
int m = (l + r) >> 1;
pushdown(u);
if(x <= m){
update(u << 1, l, m, x, v);
}
else{
update(u << 1 | 1, m + 1, r, x, v);
}
pushup(u);
}
void update(int u, const Operation &v){
update(1, 1, n, u, v);
}
Operation query(int u, int l, int r, int x, int y){
if(x <= l && r <= y) {
return op[u];
}
int m = (l + r) >> 1;
Operation cur;
pushdown(u);
if(x <= m){
cur = query(u << 1, l, m, x, y);
}
if(y > m){
cur = cur + query(u << 1 | 1, m + 1, r, x, y);
}
return cur;
}
Operation query(int l, int r){
return query(1, 1, n, l, r);
}
void range_update(int u, int l, int r, int x, int y, const Mark &v){
if(l >= x && r <= y){
modify(u, v);
return;
}
int m = (l + r) >> 1;
pushdown(u);
if(x <= m){
range_update(u << 1, l, m, x, y, v);
}
if(y > m){
range_update(u << 1 | 1, m + 1, r, x, y, v);
}
pushup(u);
}
void range_update(int l, int r, const Mark &v){
return range_update(1, 1, n, l, r, v);
}
};
struct Mark{
void modify(const Mark &v){
}
};
struct Operation {
ll sum = 0;
Operation(ll p = 0){
sum = p;
}
void modify(const Mark &v){
}
};
Operation operator+(Operation a, Operation b){
return {a.sum + b.sum};
}
void solve() {
int n;
ld ox;
cin >> n >> ox;
vector<vector<int>> ct;
vector<vector<ld>> a;
vector<int> c;
vector<ld> b;
for(int i = 0; i < n; i++){
ld x, y, z;
cin >> x >> y >> z;
ct.pb({x, y, z});
c.pb(x);
c.pb(y);
}
c.pb(0);
c.pb(ox);
map<int, int> rec;
map<ld, int> mp;
sort(c.begin(), c.end());
c.erase(unique(c.begin(), c.end()), c.end());
for(int i = 0; i < c.size(); i++){
rec[c[i]] = i;
}
a.pb({0});
Point<ld> p1 = {rec[0], rec[0]}, p2 = {rec[ox], rec[0]};
for(int i = 0; i < n; i++){
ld x = rec[ct[i][0]], y = rec[ct[i][1]], z = ct[i][2];
Point<ld> p = {x, y};
ld dl = p.get_angle(p - p1, p2 - p1) * 180;
ld dr = p.get_angle(p - p2, p1 - p2) * 180;
a.pb({dl, dr, x, y, z, i + 1.0});
b.pb(dr);
}
sort(a.begin() + 1, a.end(), cmp);
sort(b.begin(), b.end());
b.erase(unique(b.begin(), b.end()), b.end());
//二位偏序
SegTree<Operation, Mark> tr(n + 1);
ll p[n + 1];
for(int i = 0; i < n; i++){
mp[b[i]] = i + 1;
}
for(int i = n; i >= 1; i--){
p[(int)a[i][5]] = tr.query(1, mp[a[i][1]]).sum;
tr.update(mp[a[i][1]], (int)a[i][4]);
}
for(int i = 1; i <= n; i++){
cout << p[i] << endl;
}
}
int main() {
ios::sync_with_stdio(false);
cin.tie(nullptr);
cout.tie(nullptr);
int t = 1;
// cin >> t;
while (t--) {
solve();
}
return 0;
}
详细
Test #1:
score: 100
Accepted
time: 1ms
memory: 3760kb
input:
5 8 -8 1 1 -1 10 2 0 3 4 7 1 8 8 2 16
output:
0 12 0 0 8
result:
ok 5 lines
Test #2:
score: 0
Accepted
time: 0ms
memory: 3876kb
input:
6 6 0 1 1 2 3 10 2 5 100 3 1 1000 3 5 10000 4 5 100000
output:
0 1000 1010 0 1010 1000
result:
ok 6 lines
Test #3:
score: 0
Accepted
time: 302ms
memory: 39728kb
input:
99999 1000000000 500002962 1 1 500025469 1 1 500044229 1 1 500026049 1 1 499983663 1 1 499965983 1 1 499988191 1 1 499987116 1 1 500029240 1 1 499975570 1 1 499973295 1 1 499986404 1 1 500023312 1 1 499964976 1 1 499952153 1 1 500046927 1 1 499951857 1 1 499984523 1 1 500038724 1 1 499991318 1 1 500...
output:
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ...
result:
ok 99999 lines
Test #4:
score: 0
Accepted
time: 218ms
memory: 39740kb
input:
100000 1000000000 -50000 1000000000 454290650 -49999 1000000000 208284433 -49998 1000000000 854275069 -49997 1000000000 627720731 -49996 1000000000 79147837 -49995 1000000000 614585061 -49994 1000000000 438660998 -49993 1000000000 300657551 -49992 1000000000 546865509 -49991 1000000000 353401129 -49...
output:
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ...
result:
ok 100000 lines
Test #5:
score: 0
Accepted
time: 228ms
memory: 39780kb
input:
100000 1000000000 -1000000000 1 423302241 -999999999 1 941931570 -999999998 1 801255926 -999999997 1 434775469 -999999996 1 784636342 -999999995 1 41794758 -999999994 1 768189347 -999999993 1 746924545 -999999992 1 259101843 -999999991 1 798620683 -999999990 1 447243634 -999999989 1 848852324 -99999...
output:
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ...
result:
ok 100000 lines
Test #6:
score: 0
Accepted
time: 278ms
memory: 40140kb
input:
99999 1000000000 499994038 499999999 998430190 500036272 499999997 789110725 499988970 499999999 119471973 500042096 499999996 855486238 499953314 499999995 464948333 499979222 499999999 573710457 500002347 499999999 385287206 500030797 499999998 589559003 500043266 499999996 394228619 500028049 499...
output:
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ...
result:
ok 99999 lines
Test #7:
score: -100
Wrong Answer
time: 318ms
memory: 41952kb
input:
99999 1000000000 500015023 90276 306948325 499983944 103118 727377493 500025915 268634 390844401 499969478 372636 395763733 500025195 253915 906667281 500002248 2021 592484584 500028251 319247 781019435 500002485 2470 479698120 500019573 153240 55967591 499996332 5381 572934141 500029257 342388 5702...
output:
13387209411328 13031712260900 21595673148600 22291044642584 21063241148691 2225793232815 23279890220189 2451076537722 16913947498857 3471009954296 23982851387523 16626157024751 30095021973687 19926537454464 24849458461846 31917289943481 4532554141357 5472863984472 6394119982793 29002604000179 328031...
result:
wrong answer 1st lines differ - expected: '15051696338423', found: '13387209411328'