結果
問題 | No.2187 三立法和 mod 333 |
ユーザー |
![]() |
提出日時 | 2023-04-01 23:51:10 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 58 ms / 100 ms |
コード長 | 1,618 bytes |
コンパイル時間 | 2,895 ms |
コンパイル使用メモリ | 250,972 KB |
実行使用メモリ | 9,216 KB |
最終ジャッジ日時 | 2024-09-24 13:44:37 |
合計ジャッジ時間 | 5,580 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 33 |
ソースコード
#include <bits/stdc++.h>//#include <atcoder/all>using namespace std;//using namespace atcoder;//using mint = modint1000000007;//const int mod = 1000000007;//using mint = modint998244353;//const int mod = 998244353;//const int INF = 1e9;//const long long LINF = 1e18;#define rep(i, n) for (long long i = 0; i < (n); ++i)#define rep2(i,l,r)for(long long i=(l);i<(r);++i)#define rrep(i, n) for (long long i = (n-1); i >= 0; --i)#define rrep2(i,l,r)for(long long i=(r-1);i>=(l);--i)#define all(x) (x).begin(),(x).end()#define allR(x) (x).rbegin(),(x).rend()#define endl "\n"#define P pair<int,int>template<typename A, typename B> inline bool chmax(A & a, const B & b) { if (a < b) { a = b; return true; } return false; }template<typename A, typename B> inline bool chmin(A & a, const B & b) { if (a > b) { a = b; return true; } return false; }int main() {ios::sync_with_stdio(false);cin.tie(nullptr);int a; cin >> a;long long lim = 4444LL * 4444 * 4444 * 4444;vector<long long>v4(4445, 1);rep(i, v4.size()) {rep(j, 4)v4[i] *= i;}vector<vector<int>>v3(333, vector<int>(4446));rep2(i, 1, 4445) {int x = (i * i * i) % 333;v3[x][i]++;}rep(i, 333)rep(j, 4445) {v3[i][j + 1] += v3[i][j];}long long ans = 0;rep2(i, 1, 4445) {int k = 4444;//ok maxrep2(j, 1, i + 1) {if (v4[i] + v4[j] > lim)break;long long y = 3330000000000000 - i * i*i - j * j*j + a;y %= 333;while (v4[i] + v4[j] + v4[k] > lim) {k--;if (0 == k)break;}if (0 == k)break;ans += v3[y][k];if (i != j)ans += v3[y][k];}}cout << ans << endl;return 0;}