Problem p02679 (C++) — buggy submission, verdict: Runtime Error

#include <bits/stdc++.h>
using namespace std;

/*
イワシの入れ方は2^N- 1通り
ai * aj + bi * bj = 0
ai * aj = - bi * bj
biとajが0でないと仮定すると、
ai / bi = - bj / aj
割り切れるとも限らないから難しいかも

ai * aj + bi * bj = 0 かつ ai * ak + bi * bk = 0のとき、
ai * aj * ak + bi * bj * ak = 0
aj * (- bi * bk) + bi * bj * ak = 0
(bj * ak - aj * bk) * bi = 0
i) bi == 0
        ai * aj = 0, ai * ak = 0となるので、ai = 0 または aj == ak
ii) bj * ak - aj * bk == 0
        aj * bk = bj * ak
        aj : bj = ak : bkとなる

なんか比を保存しておけばよさそう
0は例外パターンで、両方0ならどれもダメ
片方0なら、もう一方が0でなければOK

全体から引いていく
その集合から2個以上選んで、他は任意を足していけばいい
ただ重複が生じてしまいそうなので、包除原理とかでできそうだけど計算量がだめ

普通に計算する
同じ集合からは1個しか出せない
選び方は集合の要素の個数

3
1 2
-1 1
2 -1
だったら...

なんか違う
A * x + B * y = 0
y = - (A / B) * x みたいになるはずで、xとyの比だけ見ればいいはず
A : B = (- Y) : X

一緒に入れちゃいけない集合が違うけど、例外除けば1:1対応になるはず
つまりパターンとしては、集合Aと集合Bが対応しているなら、集合Aのみ使うパターンか集合Bだけ使うパターンかどちらも使わないか

1:1対応にならない
{a, b}, {-b, a}, {b, -a}, {-a, -b}の4パターンが関係する
{a, b} -> {-b, a}と{b, -a}
{-b, a} -> {a, b}と{-a, -b}みたいになる
一緒にしちゃえば問題ないのでは？
aが必ず正と仮定すればいい

3
1 2
-1 1
2 -1
だったら集合は3つ
A = {{1, 2}}
B = {{2, -1}}
C = {{-1, 1}}
AとBが対応しているから、取りうるパターンは3通り
Cを使うか使わないかで2通り
0個を除くので-1
よって3 * 2 - 1 = 5?
なんかあってそう

つまり対応する集合のペアをさがしてあげればいい
ただ、{0, 0}はどれにおいてもダメなので無視して、最後に+sizeする
{0, Y}や{X, 0}については、この組み合わせがペア
*/

int mod = 1000000007;

long long gcd(long long a, long long b) {
  if (a < 0)
    a = -a;
  if (b < 0)
    b = -b;
  if (b == 0)
    return (a);
  return (gcd(b, a % b));
}

int main(void) {
  int n;
  cin >> n;
  multiset<pair<long long, long long>> iwashi;
  long long a, b;
  long long tmp_ans = 0;
  long long a_zero = 0;
  long long b_zero = 0;
  vector<long long> mod_two(n + 1, 1);
  for (int i = 1; i <= n; i++)
    mod_two[i] = mod_two[i - 1] * 2 % mod;
  for (int i = 0; i < n; i++) {
    cin >> a >> b;
    long long c = gcd(a, b);
    if (c != 0) {
      a /= c;
      b /= c;
    }
    if (a < 0) {
      a = -a;
      b = -b;
    }
    if (a == 0 && b == 0)
      tmp_ans++;
    else if (a == 0 || b == 0) {
      if (a == 0)
        a_zero++;
      if (b == 0)
        b_zero++;
    } else
      iwashi.insert({a, b});
  }
  long long ans = 1;
  ans *= (mod + mod_two[a_zero] + mod_two[b_zero] - 1) % mod;
  ans %= mod;
  map<pair<long long, long long>, bool> mp;
  auto ptr = iwashi.begin();
  int num, count;
  while (ptr != iwashi.end()) {
    a = (*ptr).first;
    b = (*ptr).second;
    if (!mp[{a, b}]) {
      num = iwashi.count({a, b});
      count = 0;
      count += iwashi.count({b, -a});
      count += iwashi.count({-b, a});
      ans *= (mod + mod_two[num] + mod_two[count] - 1) % mod;
      ans %= mod;
      mp[{a, b}] = true;
      mp[{b, -a}] = true;
      mp[{-b, a}] = true;
    }
    ptr++;
  }
  ans += tmp_ans;
  ans %= mod;
  cout << (mod + ans - 1) % mod << endl;
  return (0);
}


Error: 0

Fixed submission:

#include <bits/stdc++.h>
using namespace std;

/*
イワシの入れ方は2^N- 1通り
ai * aj + bi * bj = 0
ai * aj = - bi * bj
biとajが0でないと仮定すると、
ai / bi = - bj / aj
割り切れるとも限らないから難しいかも

ai * aj + bi * bj = 0 かつ ai * ak + bi * bk = 0のとき、
ai * aj * ak + bi * bj * ak = 0
aj * (- bi * bk) + bi * bj * ak = 0
(bj * ak - aj * bk) * bi = 0
i) bi == 0
        ai * aj = 0, ai * ak = 0となるので、ai = 0 または aj == ak
ii) bj * ak - aj * bk == 0
        aj * bk = bj * ak
        aj : bj = ak : bkとなる

なんか比を保存しておけばよさそう
0は例外パターンで、両方0ならどれもダメ
片方0なら、もう一方が0でなければOK

全体から引いていく
その集合から2個以上選んで、他は任意を足していけばいい
ただ重複が生じてしまいそうなので、包除原理とかでできそうだけど計算量がだめ

普通に計算する
同じ集合からは1個しか出せない
選び方は集合の要素の個数

3
1 2
-1 1
2 -1
だったら...

なんか違う
A * x + B * y = 0
y = - (A / B) * x みたいになるはずで、xとyの比だけ見ればいいはず
A : B = (- Y) : X

一緒に入れちゃいけない集合が違うけど、例外除けば1:1対応になるはず
つまりパターンとしては、集合Aと集合Bが対応しているなら、集合Aのみ使うパターンか集合Bだけ使うパターンかどちらも使わないか

1:1対応にならない
{a, b}, {-b, a}, {b, -a}, {-a, -b}の4パターンが関係する
{a, b} -> {-b, a}と{b, -a}
{-b, a} -> {a, b}と{-a, -b}みたいになる
一緒にしちゃえば問題ないのでは？
aが必ず正と仮定すればいい

3
1 2
-1 1
2 -1
だったら集合は3つ
A = {{1, 2}}
B = {{2, -1}}
C = {{-1, 1}}
AとBが対応しているから、取りうるパターンは3通り
Cを使うか使わないかで2通り
0個を除くので-1
よって3 * 2 - 1 = 5?
なんかあってそう

つまり対応する集合のペアをさがしてあげればいい
ただ、{0, 0}はどれにおいてもダメなので無視して、最後に+sizeする
{0, Y}や{X, 0}については、この組み合わせがペア
*/

int mod = 1000000007;

long long gcd(long long a, long long b) {
  if (a < 0)
    a = -a;
  if (b < 0)
    b = -b;
  if (b == 0)
    return (a);
  return (gcd(b, a % b));
}

int main(void) {
  int n;
  cin >> n;
  multiset<pair<long long, long long>> iwashi;
  long long a, b;
  long long tmp_ans = 0;
  long long a_zero = 0;
  long long b_zero = 0;
  vector<long long> mod_two(n + 1, 1);
  for (int i = 1; i <= n; i++)
    mod_two[i] = mod_two[i - 1] * 2 % mod;
  for (int i = 0; i < n; i++) {
    cin >> a >> b;
    long long c = gcd(a, b);
    a /= c;
    b /= c;
    if (a < 0) {
      a = -a;
      b = -b;
    }
    if (a == 0 && b == 0)
      tmp_ans++;
    if (a == 0 || b == 0) {
      if (a == 0)
        a_zero++;
      if (b == 0)
        b_zero++;
    } else
      iwashi.insert({a, b});
  }
  long long ans = 1;
  ans *= (mod + mod_two[a_zero] + mod_two[b_zero] - 1) % mod;
  ans %= mod;
  auto ptr = iwashi.begin();
  while (ptr != iwashi.end()) {
    a = (*ptr).first;
    b = (*ptr).second;
    int num = iwashi.count({a, b});
    int count = 0;
    count += iwashi.count({b, -a});
    iwashi.erase({b, -a});
    count += iwashi.count({-b, a});
    iwashi.erase({-b, a});
    ans *= (mod + mod_two[num] + mod_two[count] - 1) % mod;
    ans %= mod;
    ptr++;
  }
  ans += tmp_ans;
  ans %= mod;
  cout << (mod + ans - 1) % mod << endl;
  return (0);
}