結果

問題 No.2651 [Cherry 6th Tune B] $\mathbb{C}$omplex комбинат
ユーザー chro_96chro_96
提出日時 2024-02-23 23:02:34
言語 C
(gcc 13.3.0)
結果
AC  
実行時間 358 ms / 2,500 ms
コード長 1,503 bytes
コンパイル時間 1,372 ms
コンパイル使用メモリ 31,872 KB
実行使用メモリ 6,656 KB
最終ジャッジ日時 2024-09-29 08:19:34
合計ジャッジ時間 13,241 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 39
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <stdio.h>

long long power_mod (long long a, long long b, long long mod_num) {
  long long ans = 1LL;
  
  if (b > 0LL) {
    ans = power_mod(a, b/2LL, mod_num);
    ans = (ans * ans) % mod_num;
    if (b%2LL == 1LL) {
      ans = (ans * (a % mod_num)) % mod_num;
    }
  }
  
  return ans;
}

int main () {
  int t = 0;
  int n = 0;
  long long x[300000] = {};
  long long y[300000] = {};
  
  int res = 0;
  
  long long mod_num = 998244353LL;
  
  res = scanf("%d", &t);
  while (t > 0) {
    long long ans = 0LL;
    long long acc[6] = {};
    res = scanf("%d", &n);
    for (int i = 0; i < n; i++) {
      res = scanf("%lld", x+i);
      res = scanf("%lld", y+i);
    }
    for (int i = n-1; i >= 0; i--) {
      long long tmp[6] = {};
      long long abs_sq = (x[i]*x[i]+y[i]*y[i])%mod_num;
      long long inv = power_mod(abs_sq, mod_num-2LL, mod_num);
      tmp[0] = inv;
      tmp[1] = abs_sq;
      tmp[2] = ((x[i]*x[i]+mod_num-(y[i]*y[i])%mod_num)*inv)%mod_num;
      tmp[3] = ((2LL*x[i]*y[i])*inv)%mod_num;
      tmp[4] = tmp[2];
      tmp[5] = (mod_num-tmp[3])%mod_num;
      ans += (tmp[0]*acc[1])%mod_num;
      ans += (tmp[1]*acc[0])%mod_num;
      ans += mod_num-(tmp[2]*acc[4])%mod_num;
      ans += (tmp[3]*acc[5])%mod_num;
      ans += mod_num-(tmp[4]*acc[2])%mod_num;
      ans += (tmp[5]*acc[3])%mod_num;
      for (int j = 0; j < 6; j++) {
        acc[j] += tmp[j];
        acc[j] %= mod_num;
      }
    }
    printf("%lld\n", ans%mod_num);
    t--;
  }
  
  return 0;
}
0