結果
問題 |
No.3033 エルハートの数え上げ
|
ユーザー |
|
提出日時 | 2025-02-21 22:42:07 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 11 ms / 2,000 ms |
コード長 | 1,195 bytes |
コンパイル時間 | 3,835 ms |
コンパイル使用メモリ | 278,628 KB |
実行使用メモリ | 6,824 KB |
最終ジャッジ日時 | 2025-02-21 22:42:15 |
合計ジャッジ時間 | 5,198 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 20 |
ソースコード
#include<bits/stdc++.h> using namespace std; #define rep(i,a,b) for(int i=(a);i<(b);i++) #define all(a) begin(a),end(a) #define sz(a) (int)(a).size() typedef long long ll; typedef vector<int> vi; typedef pair<int,int> pii; const ll mod=998244353; ll pow_mod(ll a,ll n){ a%=mod; ll res=1; while(n){ if(n%2)res=res*a%mod; a=a*a%mod;n/=2; } if(res<0)res+=mod; return res; } int main(){ cin.tie(0)->sync_with_stdio(0); cin.exceptions(cin.failbit); int N,M;cin>>N>>M; vector<ll>A(M),B(M),C(M),D(M); rep(i,0,M)cin>>A[i]>>B[i]>>C[i]>>D[i]; ll a=0,b=0,c=0,d=0; //f(1),f(2),f(-1),f(-2) rep(x,-15,16)rep(y,-15,16)rep(z,-15,16){ bool Int=1,Cl=1; rep(i,0,M){ if(A[i]*x+B[i]*y+C[i]*z+D[i]<0)Cl=0; if(A[i]*x+B[i]*y+C[i]*z+D[i]<=0)Int=0; } if(Int)c--; if(Cl)a++; } rep(x,-30,31)rep(y,-30,31)rep(z,-30,31){ bool Int=1,Cl=1; rep(i,0,M){ if(A[i]*x+B[i]*y+C[i]*z+2*D[i]<0)Cl=0; if(A[i]*x+B[i]*y+C[i]*z+2*D[i]<=0)Int=0; } if(Int)d--; if(Cl)b++; } ll p=8*a-2*b+8*c-2*d; ll q=8*a-b-8*c+d; ll r=-2*a+2*b-2*c+2*d; ll s=-2*a+b+2*c-d; ll f=p-N%mod*q%mod+N%mod*N%mod*r%mod-N%mod*N%mod*N%mod*s%mod; f*=-pow_mod(12,mod-2); f%=mod; if(f<0)f+=mod; cout<<f<<endl; }