結果

問題 No.2348 Power!! (Easy)
ユーザー maksimmaksim
提出日時 2023-06-09 22:57:14
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 449 ms / 5,000 ms
コード長 1,360 bytes
コンパイル時間 1,660 ms
コンパイル使用メモリ 168,616 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-06-10 14:07:08
合計ジャッジ時間 6,305 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 11 ms
5,248 KB
testcase_01 AC 182 ms
5,376 KB
testcase_02 AC 171 ms
5,376 KB
testcase_03 AC 182 ms
5,376 KB
testcase_04 AC 27 ms
5,376 KB
testcase_05 AC 448 ms
5,376 KB
testcase_06 AC 449 ms
5,376 KB
testcase_07 AC 388 ms
5,376 KB
testcase_08 AC 394 ms
5,376 KB
testcase_09 AC 397 ms
5,376 KB
testcase_10 AC 388 ms
5,376 KB
testcase_11 AC 393 ms
5,376 KB
testcase_12 AC 387 ms
5,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;
#define int long long
const int p=998244353;
int po(int a,int b){
  int x = a, res = 1;
  int h = b;
  while (h > 0) {
    if (h & 1) {res *= x;res%=p;}
    x *= x;x%=p;
    h >>= 1;
  }
  return res;
}
int inverse(int a,int m) {
  int u = 0, v = 1;
  while (a != 0) {
    int t = m / a;
    m -= t * a; swap(a, m);
    u -= t * v; swap(u, v);
  }
  return u;
}
int32_t main()
{
    ios_base::sync_with_stdio(false);cin.tie(0);cout.tie(0);
    int h=(1<<12)*7*17;int sz=((p-1)/h);
    int t;cin>>t;while(t--)
    {
    int a,n;cin>>a>>n;
    int pos[sz];for(int i=0;i<sz;++i) pos[i]=po(a,i*h);
    int invs[sz];for(int i=0;i<sz;++i) invs[i]=inverse(pos[i]-1,p);
    int res=0;
    int cur1=1;int cur2=1;int t=1;int u1=1;int u=1;
    for(int i=0;i<h;++i)
    {
        if(i>=n) continue;
        int o=(n-1-i)/h+1;
        int l=i*i;int r=i*i+o*i*2*h;
        int step=2*i*h;int u=pos[(2*i)%sz];
        if(i==0 || (n-i-1)/h!=(n-i)/h)
        {
            t=po(a,o*i*2*h);u1=po(a,o*2*h);
        }
        else
        {
            t*=u1;t%=p;
        }
        if(u!=1) {int inva=invs[(2*i)%sz];int v=(cur1*t)%p;res+=(v-cur1)*inva;res%=p;}
        else {res+=cur1*o;res%=p;}
        cur1*=cur2;cur1%=p;cur1*=cur2;cur1%=p;cur1*=a;cur1%=p;cur2*=a;cur2%=p;
    }
    cout<<(res%p+p)%p<<'\n';
    }
    return 0;
}
0