結果

問題 No.3403 Count 1210 Sequence
コンテスト
ユーザー Mike scarf
提出日時 2025-12-12 16:04:40
言語 C++23
(gcc 13.3.0 + boost 1.89.0)
結果
AC  
実行時間 1,297 ms / 2,000 ms
コード長 1,355 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,921 ms
コンパイル使用メモリ 276,484 KB
実行使用メモリ 8,476 KB
最終ジャッジ日時 2025-12-12 16:05:13
合計ジャッジ時間 31,487 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 31
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:27:10: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   27 |     scanf("%lld",&t);
      |     ~~~~~^~~~~~~~~~~
main.cpp:29:14: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   29 |         scanf("%lld%lld",&n,&a);
      |         ~~~~~^~~~~~~~~~~~~~~~~~

ソースコード

diff #
raw source code 

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const ll mod=998244353;
const int maxn=300005;
ll fac[maxn+5],inv[maxn+5];
ll qmksm(ll a,ll b){
    ll res;
    res=1;
    while(b){
        if(b&1) res=res*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return res;
}
ll C(ll n,ll k){
    if(k<0 or k>n) return 0;
    return fac[n]*inv[k]%mod*inv[n-k]%mod;
}
int main(){
    ll t,n,a,ans,i,k,cur,zc,zc1,zc2;
    fac[0]=1;
    for(i=1;i<=maxn;i++) fac[i]=fac[i-1]*i%mod;
    inv[maxn]=qmksm(fac[maxn],mod-2);
    for(i=maxn-1;i>=0;i--) inv[i]=inv[i+1]*(i+1)%mod;
    scanf("%lld",&t);
    while(t--){
        scanf("%lld%lld",&n,&a);
        ans=0;
        cur=n-1;
        for(i=1;i*i<=a;i++){
            if(a%i==0){
                k=i;
                zc=a/k;
                if((cur-zc)%2==0){
                    zc1=(cur-zc)/2;
                    zc2=(cur+zc+2)/2;
                    ans=(ans+C(cur,zc1)-C(cur,zc2)+mod)%mod;
                }
                if(i*i!=a){
                    k=a/i;
                    zc=a/k;
                    if((cur-zc)%2==0){
                        zc1=(cur-zc)/2;
                        zc2=(cur+zc+2)/2;
                        ans=(ans+C(cur,zc1)-C(cur,zc2)+mod)%mod;
                    }
                }
            }
        }
        printf("%lld\n",ans);
    }
    return 0;
}
0