Cod sursa(job #2629897)

Utilizator segtreapMihnea Andreescu segtreap Data 23 iunie 2020 10:53:17
Problema Regiuni Scor 40
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 2.83 kb
// while (clock()<=69*CLOCKS_PER_SEC)
// #pragma comment(linker, "/stack:200000000")
// #pragma GCC optimize("O3")
// #pragma GCC target ("avx2")
// #pragma GCC optimize("Ofast")
// #pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")
// #pragma GCC optimize("unroll-loops")
// #include <bits/stdc++.h>
// #include <ext/pb_ds/assoc_container.hpp>
// #include <ext/pb_ds/tree_policy.hpp>
// #include <cmath>
// #include <functional>
// #include <fstream>
// #include <algorithm>
// #include <string>
// #include <set>
// #include <map>
// #include <list>
// #include <time.h>
// #include <math.h>
// #include <random>
// #include <deque>
// #include <queue>
// #include <cassert>
// #include <unordered_map>
// #include <unordered_set>
// #include <iomanip>
// #include <bitset>
// #include <sstream>
// #include <chrono>
// #include <cstring>
///#include <vector>
#include <cstdio>
#include <iostream>
#include <vector>

using namespace std;

typedef long long ll;
typedef long double ld;
///mt19937 rng(chrono::high_resolution_clock::now().time_since_epoch().count());
bool HOME = 1;
const string TASK_NAME = "regiuni";
#define files(TASK_NAME) read(TASK_NAME); write(TASK_NAME)
#define read(TASK_NAME) freopen((TASK_NAME + ".in").c_str(), "r", stdin)
#define write(TASK_NAME) freopen((TASK_NAME + ".out").c_str(), "w", stdout)
#define all(c) (c).begin(), (c).end()
#define sz(c) (int) (c).size()
#define left left228
#define right right228
#define y1 y1228
#define mp make_pair
#define pb push_back
#define y2 y2228
#define rank rank228

struct Point {
  ld x, y;
};

ll f(Point a, Point b) {
  return (a.x - b.x) * (a.y + b.y);
}

ll f(Point a, Point b, Point c) {
  return f(a, b) + f(b, c) + f(c, a);
}

const int N = 1000 + 7;
int n, m;
ld a[N], b[N], c[N];
Point p[N];
vector<int> co[N];

int main() {
#ifdef INFOARENA
    HOME = 0;
#endif
  if (HOME) {
    freopen ("input", "r", stdin);
  } else {
    files(TASK_NAME);
  }
  cin >> n >> m;
  for (int i = 1; i <= n; i++) {
    cin >> a[i] >> b[i] >> c[i];
  }
  for (int i = 1; i <= m; i++) {
    cin >> p[i].x >> p[i].y;
    co[1].push_back(i);
  }
  int cnt = 1;
  for (int i = 1; i <= n; i++) {
    Point p1, p2;
    if (b[i] == 0) {
      p1 = {-c[i] / a[i], 0};
      p2 = {-c[i] / a[i], 1};
    } else {
      p1 = {0, -(c[i]) / b[i]};
      p2 = {1, -(c[i] + a[i]) / b[i]};
    }
    int cnt2 = cnt;
    for (int j = 1; j <= cnt; j++) {
      vector<int> a, b;
      for (auto &it : co[j]) {
        if (f(p1, p2, p[it]) < 0) {
          a.push_back(it);
        } else {
          b.push_back(it);
        }
      }
      if (a.empty() || b.empty()) {
        continue;
      }
      co[j] = a;
      co[++cnt2] = b;
    }
    cnt = cnt2;
  }
  cout << cnt << "\n";
  return 0;
}