結果

問題 No.2187 三立法和 mod 333
ユーザー Aging1986
提出日時 2023-01-13 22:05:52
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 31 ms / 100 ms
コード長 2,265 bytes
コンパイル時間 1,920 ms
コンパイル使用メモリ 196,064 KB
最終ジャッジ日時 2025-02-10 02:35:35
ジャッジサーバーID
(参考情報)
judge5 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 33
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<ll, ll> pll;
typedef pair<int, int> pii;
typedef pair<double, double> pdd;
typedef vector<int> vi;
typedef vector<ll> vll;
typedef vector<double> vd;
typedef vector<vi> vvi;
typedef vector<vll> vvll;
typedef vector<pii> vpii;
typedef vector<pll> vpll;
typedef vector<pdd> vpdd;
typedef vector<vd> vvd;
#define yn(ans) printf("%s\n", (ans)?"Yes":"No");
#define YN(ans) printf("%s\n", (ans)?"YES":"NO");
template<class T> bool chmax(T &a, T b) {
	if (a >= b) return false;
	a = b; return true;
}
template<class T> bool chmin(T &a, T b) {
	if (a <= b) return false;
	a = b; return true;
}
#define FOR(i, s, e, t) for ((i) = (s); (i) < (e); (i) += (t)) 
#define REP(i, e) for (int i = 0; i < (e); ++i) 
#define REP1(i, s, e) for (int i = (s); i < (e); ++i)
#define RREP(i, e) for (int i = (e); i >= 0; --i)
#define RREP1(i, e, s) for (int i = (e); i >= (s); --i)
#define all(v) v.begin(), v.end()
#define pb push_back
#define qb pop_back
#define pf push_front
#define qf pop_front
#define maxe max_element
#define mine min_element
ll inf = 1e18;
#define DEBUG printf("%d\n", __LINE__); fflush(stdout);
template<class T> void print(vector<T> &v, bool withSize) {
	if (withSize) cout << v.size() << endl;
	REP(i, v.size()) cout << v[i] << " "; 
	cout << "\n";
}
mt19937_64 rng((unsigned int) chrono::steady_clock::now().time_since_epoch().count());

int __FAST_IO__ = []() {
	std::ios::sync_with_stdio(0);
	std::cin.tie(0);
	std::cout.tie(0);
	return 0;
}();

int main() {
    int A;
    cin >> A;
    vvi p(333);
    REP1(i, 1, 4445) {
    	ll q = 1ll * i * i % 333 * i % 333;
    	p[q].pb(i);
    }
    
    ll tot = 1ll * 4444 * 4444 * 4444 * 4444;
    ll ans = 0;
    REP(i, 333) REP(j, 333) {
    	int t = (A + 666 - i - j) % 333;
    	for (auto a: p[i]) {
    		for (int k = 0, l = (int)p[t].size() - 1; k < p[j].size(); ++k) {
    			int b = p[j][k];
    			ll sum = 1ll * a * a * a * a + 1ll * b * b * b * b;
    			if (sum > tot) break;
    			while (l >= 0 && sum + 1ll * p[t][l] * p[t][l] * p[t][l] * p[t][l] > tot) --l;
    			if (l == -1) break;
    			ans += l + 1;
    		}
    	}
   // 	printf("%d %d %lld\n", i, j, ans);
    }
    
    printf("%lld\n", ans);
    return 0;
}
0