結果
| 問題 |
No.2406 Difference of Coordinate Squared
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2023-08-04 23:23:56 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 1,030 bytes |
| コンパイル時間 | 3,887 ms |
| コンパイル使用メモリ | 254,628 KB |
| 最終ジャッジ日時 | 2025-02-15 23:17:27 |
|
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 47 WA * 8 |
ソースコード
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using ll = long long;
#define rep(i,n) for(int i=0;i<(int)(n);i++)
using mint = atcoder::modint998244353;
const int MOD=998244353;
const int MAX_N=2e6;
vector<ll> fac(MAX_N+1,1);
vector<ll> finv(MAX_N+1,1);
vector<ll> inv(MAX_N+1,1);
void comb_setup(){
for(int i=2;i<=MAX_N;i++){
fac.at(i)=(fac.at(i-1)*i)%MOD;
inv.at(i)=MOD-(inv.at(MOD%i)*(MOD/i))%MOD;
finv.at(i)=(finv.at(i-1)*inv.at(i))%MOD;
}
}
ll comb(int n,int k){
if(n<k) return 0;
return fac.at(n)*(finv.at(n-k)*finv.at(k)%MOD)%MOD;
}
vector<pair<ll,ll>> dv(ll n){
vector<pair<ll,ll>> rt;
for(ll i=1;i*i<=n;i++){
if(n%i==0){
rt.push_back({i,n/i});
}
}
return rt;
}
int main(){
comb_setup();
ll n,m;
cin>>n>>m;
m=abs(m);
vector<pair<ll,ll>> vp=dv(m);
mint ans=0;
mint dd=mint(4).inv().pow(n);
for(auto[v,u]:vp){
if((v+u)%2) continue;
if(u>n) continue;
if((n-u)%2) continue;
ans+=dd*comb(n,(n-u)/2)*comb(n,(n+v)/2)*(u==v?1:2);
}
ans*=2;
cout<<ans.val()<<endl;
}