結果

問題 No.74 貯金箱の退屈
ユーザー fumiphysfumiphys
提出日時 2019-04-15 20:17:42
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 3 ms / 5,000 ms
コード長 2,762 bytes
コンパイル時間 1,820 ms
コンパイル使用メモリ 178,756 KB
実行使用メモリ 4,348 KB
最終ジャッジ日時 2023-10-22 06:46:10
合計ジャッジ時間 3,360 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,348 KB
testcase_01 AC 2 ms
4,348 KB
testcase_02 AC 2 ms
4,348 KB
testcase_03 AC 2 ms
4,348 KB
testcase_04 AC 2 ms
4,348 KB
testcase_05 AC 2 ms
4,348 KB
testcase_06 AC 2 ms
4,348 KB
testcase_07 AC 2 ms
4,348 KB
testcase_08 AC 2 ms
4,348 KB
testcase_09 AC 2 ms
4,348 KB
testcase_10 AC 2 ms
4,348 KB
testcase_11 AC 2 ms
4,348 KB
testcase_12 AC 2 ms
4,348 KB
testcase_13 AC 2 ms
4,348 KB
testcase_14 AC 3 ms
4,348 KB
testcase_15 AC 2 ms
4,348 KB
testcase_16 AC 2 ms
4,348 KB
testcase_17 AC 2 ms
4,348 KB
testcase_18 AC 2 ms
4,348 KB
testcase_19 AC 2 ms
4,348 KB
testcase_20 AC 2 ms
4,348 KB
testcase_21 AC 2 ms
4,348 KB
testcase_22 AC 2 ms
4,348 KB
testcase_23 AC 2 ms
4,348 KB
testcase_24 AC 2 ms
4,348 KB
testcase_25 AC 2 ms
4,348 KB
testcase_26 AC 2 ms
4,348 KB
testcase_27 AC 2 ms
4,348 KB
testcase_28 AC 2 ms
4,348 KB
testcase_29 AC 2 ms
4,348 KB
testcase_30 AC 2 ms
4,348 KB
testcase_31 AC 2 ms
4,348 KB
testcase_32 AC 2 ms
4,348 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

// includes
#include <bits/stdc++.h>

// macros
#define ll long long int
#define pb emplace_back
#define mk make_pair
#define pq priority_queue
#define FOR(i, a, b) for(int i=(a);i<(b);++i)
#define rep(i, n) FOR(i, 0, n)
#define rrep(i, n) for(int i=((int)(n)-1);i>=0;i--)
#define all(x) (x).begin(),(x).end()
#define sz(x) ((int)(x).size())
#define UNIQUE(v) v.erase(unique(v.begin(), v.end()), v.end())
#define FI first
#define SE second
using namespace std;

//  types
typedef pair<int, int> P;
typedef pair<ll, int> Pl;
typedef pair<ll, ll> Pll;
typedef pair<double, double> Pd;
 
// constants
const int inf = 1e9;
const ll linf = 1LL << 50;
const double EPS = 1e-10;
const int mod = 1e9 + 7;

// solve
template <class T>bool chmax(T &a, const T &b){if(a < b){a = b; return 1;} return 0;}
template <class T>bool chmin(T &a, const T &b){if(a > b){a = b; return 1;} return 0;}
template <typename T> istream &operator>>(istream &is, vector<T> &vec){for(auto &v: vec)is >> v; return is;}


const int MAX_H = 510, MAX_W = 510;
struct MatrixB {
  int H, W;
  bitset<MAX_W> val[MAX_H];
  MatrixB(int m = 1, int n = 1): H(m), W(n){}
  inline bitset<MAX_W>& operator[](size_t i){return val[i];}
  int gauss_jordan(bool is_extended = false){
    int rank = 0;
    for(int col = 0; col < W; col++){
      if(col == W - 1 && is_extended)break;
      int pivot = -1;
      for(int row = rank; row < H; row++){
        if(val[row][col]){
          pivot = row;
          break;
        }
      }
      if(pivot == -1)continue;
      swap(val[pivot], val[rank]);
      for(int row = 0; row < H; row++){
        if(row != rank && val[row][col])val[row] ^= val[rank];
      }
      rank++;
    }
    return rank;
  }
};

int linear_eq(MatrixB &A, vector<int> b, vector<int> & res){
  int m = A.H, n = A.W;
  MatrixB M(m, n + 1);
  for(int i = 0; i < m; i++){
    for(int j = 0; j < n; j++){
      M[i][j] = A[i][j];
    }
    M[i][n] = b[i];
  }
  int rank = M.gauss_jordan(true);

  for(int row = rank; row < m; row++){
    if(M[row][n])return -1;
  }

  res.assign(n, 0);
  for(int i = 0; i < rank; i++)res[i] = M[i][n];
  return rank;
}

set<int> v[101];

int main(int argc, char const* argv[])
{
  ios_base::sync_with_stdio(false);
  cin.tie(0);
  int n;
  cin >> n;
  vector<int> d(n), w(n);
  cin >> d >> w;
  rep(i, n){
    int in = (i + d[i]) % n;
    v[in].insert(i);
    in = (i - d[i]) % n;
    if(in < 0)in += n;
    v[in].insert(i);
  }
  MatrixB M(n, n);
  vector<int> b(n);
  rep(i, n){
    for(auto itr = v[i].begin(); itr != v[i].end(); ++itr){
      M[i][*itr] = 1;
    }
    if(w[i] == 0)b[i] = 1;
    else b[i] = 0;
  }
  vector<int> res;
  int rank = linear_eq(M, b, res);
  if(rank == -1)cout << "No" << endl;
  else cout << "Yes" << endl;
	return 0;
}
0