結果

問題 No.2748 Strange Clock
ユーザー Yu_212
提出日時 2024-04-21 12:24:33
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 1,550 bytes
コンパイル時間 4,896 ms
コンパイル使用メモリ 278,780 KB
実行使用メモリ 152,480 KB
最終ジャッジ日時 2024-10-13 10:51:13
合計ジャッジ時間 11,872 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample -- * 2
other AC * 26 TLE * 1 -- * 10
権限があれば一括ダウンロードができます

ソースコード

diff # 

#pragma GCC target("avx2")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")

#include<bits/stdc++.h>
using namespace std;
using ll = long long;
using pll = pair<ll, ll>;
const ll mod = 998244353;
int n;

ll extgcd(ll a, ll b) {
  ll x0 = 1, x1 = 0, y0 = 0, y1 = 1;
  while (b != 0) {
    ll q = a / b;
    ll r = a % b;
    ll x2 = x0 - q * x1;
    ll y2 = y0 - q * y1;
    a = b;
    b = r;
    x0 = x1;
    x1 = x2;
    y0 = y1;
    y1 = y2;
  }
  return x0;
}

ll crt(ll r1, ll m1, ll r2, ll m2) {
  ll dpq1 = extgcd(m1, m2);
  ll m = m1 * m2;
  ll temp = (r2 - r1) * dpq1 % m2;
  return ((r1 + m1 * temp) % m + m) % m;
}


ll conv(ll v, int c) {
  ll r = 0;
  ll pc = 1;
  for (int i = 0; i < n; i++) {
    r += v % 3 * pc;
    pc *= c;
    v /= 3;
  }
  return r;
}

ll pw(ll a, ll b) {
  ll x = 1;
  while (b > 0) {
    if ((b & 1) == 1) {
      x = x * a;
    }
    a = a * a;
    b >>= 1;
  }
  return x;
}

int main() {
  cin.tie(0)->sync_with_stdio(false);
  cin >> n;
  ll m;
  cin >> m;
  ll p3 = pw(3, n);
  ll p4 = pw(4, n);
  ll p6 = pw(6, n);
  ll p12 = pw(12, n);
  map<ll, vector<ll>> mp;
  for (ll v3 = 0; v3 < p3; v3++) {
    ll v4 = conv(v3, 4);
    ll v6 = conv(v3, 6);
    ll v12 = crt(v3, p3, v4, p4);
    ll diff = (v12 % p6 + p6 - v6) % p6;
    mp[diff].push_back(v12);
  }
  ll ans = 0;
  for (auto &[k, vv]: mp) {
    sort(vv.begin(), vv.end());
    ll last = vv[vv.size() - 1] - p12;
    for (ll v : vv) {
      if (v - m > last) {
        ans++;
      }
      last = v;
    }
  }
  cout << ans << endl;
}
0