답안 #260154

# 제출 시각 아이디 문제 언어 결과 실행 시간 메모리
260154 2020-08-09T11:52:30 Z tincamatei Printed Circuit Board (CEOI12_circuit) C++14
40 / 100
100 ms 12388 KB
#include <bits/stdc++.h>

const int MAX_N = 100000;

struct Fraction {
  int a, b;

  Fraction() {
    a = 0;
    b = 1;
  }

  Fraction(int _a, int _b) {
    a = _a;
    b = _b;
    if(b < 0) {
      a = -a;
      b = -b;
    }
  }

  bool operator== (const Fraction x) const {
    return (long long)a * x.b == (long long)b * x.a;
  }

  bool operator< (const Fraction x) const {
    return (long long)a * x.b < b * x.a;
  }
};

struct Point {
  int x, y;
  Fraction slope;
} points[MAX_N];
bool marked[MAX_N];

static inline long long detarea(int a, int b) {
  return (long long)points[a].x * points[b].y - (long long)points[a].y * points[b].x;
}

static inline long long quadarea(int a, int b, int c, int d) {
  return detarea(a, b) + detarea(b, c) + detarea(c, d) + detarea(d, a);
}

static inline long long triarea(int a, int b, int c) {
  return detarea(a, b) + detarea(b, c) + detarea(c, a);
}

struct Segment {
  int a, b;

  bool operator< (const Segment x) const {
    if(b == x.a)
      return true;
    if(a == x.b)
      return false;
    long long area = quadarea(a, b, x.b, x.a);
    if(area == 0)
      return points[a].x < points[x.a].x;
    return area < 0; 
  }
} segments[MAX_N];

struct Event {
  Fraction t;
  int type, segmentId;
};

bool cmp(Event a, Event b) {
  return a.t < b.t || (a.t == b.t && a.type < b.type);
}

int main() {
  int N;
  scanf("%d", &N);
  for(int i = 0; i < N; ++i) {
    scanf("%d%d", &points[i].x, &points[i].y);
    points[i].slope = Fraction(points[i].y, points[i].x);
    segments[i] = {i, (i + 1) % N};
  }

  std::vector<Event> ev;

  for(int i = 0; i < N; ++i) {
    if(points[segments[i].b].slope < points[segments[i].a].slope)
      std::swap(segments[i].a, segments[i].b);
    else if(points[segments[i].b].slope == points[segments[i].a].slope &&
            points[segments[i].b].x > points[segments[i].a].x)
      std::swap(segments[i].a, segments[i].b);
      
    ev.push_back({points[segments[i].a].slope, 0, i});
    ev.push_back({points[segments[i].b].slope, 2, i});
    ev.push_back({points[i].slope, 1, i});
  }

  std::sort(ev.begin(), ev.end(), cmp);
  std::set<Segment> segs;

  for(int i = 0; i < 3 * N; ++i) {
    //printf("Event: %d %d-%d\n", ev[i].type,
    //                            segments[ev[i].segmentId].a + 1,
    //                            segments[ev[i].segmentId].b + 1);
    if(ev[i].type == 0) {
      segs.insert(segments[ev[i].segmentId]);
    } else if(ev[i].type == 2)
      segs.erase(segments[ev[i].segmentId]);
    else {
      Segment firstseg = *segs.begin();
      if(points[firstseg.a].slope == points[firstseg.b].slope)
        marked[ev[i].segmentId] = (ev[i].segmentId == firstseg.a);
      else
        marked[ev[i].segmentId] = (ev[i].segmentId == firstseg.a ||
                                   ev[i].segmentId == firstseg.b);
    }

    //printf("Debug segments:\n");
    //for(auto it: segs) {
    //  printf("[%d %d]\n", it.a + 1, it.b + 1);
    //}
  }

  //printf("??%d %d\n", segments[2] < segments[3], segments[3] < segments[2]);
  //printf("%lld\n", quadarea(segments[2].a, segments[2].b, segments[3].b, segments[3].a));
  //printf("%lld\n", quadarea(segments[2].a, segments[2].b, segments[3].a, segments[3].b));

  int res = 0;
  for(int i = 0; i < N; ++i)
    if(marked[i])
      ++res;

  printf("%d\n", res);
  for(int i = 0; i < N; ++i)
    if(marked[i])
      printf("%d ", i + 1);
	return 0;
}

Compilation message

circuit.cpp: In function 'int main()':
circuit.cpp:75:8: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
   scanf("%d", &N);
   ~~~~~^~~~~~~~~~
circuit.cpp:77:10: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
     scanf("%d%d", &points[i].x, &points[i].y);
     ~~~~~^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# 결과 실행 시간 메모리 Grader output
1 Correct 2 ms 1920 KB Output is correct
2 Correct 2 ms 2048 KB Output is correct
3 Correct 4 ms 2304 KB Output is correct
4 Correct 6 ms 2304 KB Output is correct
5 Incorrect 12 ms 2936 KB Output isn't correct
6 Incorrect 14 ms 2936 KB Output isn't correct
7 Incorrect 22 ms 3448 KB Output isn't correct
8 Incorrect 9 ms 2808 KB Output isn't correct
9 Correct 11 ms 2684 KB Output is correct
10 Incorrect 11 ms 2936 KB Output isn't correct
11 Incorrect 12 ms 3320 KB Output isn't correct
12 Correct 20 ms 3448 KB Output is correct
13 Incorrect 36 ms 4840 KB Output isn't correct
14 Correct 45 ms 4852 KB Output is correct
15 Correct 50 ms 7152 KB Output is correct
16 Execution timed out 111 ms 12244 KB Time limit exceeded
17 Execution timed out 125 ms 12388 KB Time limit exceeded
18 Execution timed out 33 ms 4088 KB Time limit exceeded (wall clock)
19 Execution timed out 33 ms 4088 KB Time limit exceeded (wall clock)
20 Execution timed out 31 ms 4224 KB Time limit exceeded (wall clock)