結果

問題 No.3033 エルハートの数え上げ
ユーザー TKTYI
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#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;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0