結果

問題 No.2651 [Cherry 6th Tune B] C\mathbb{C}omplex комбинат
ユーザー 沙耶花
提出日時 2024-02-23 21:54:26
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,069 ms / 2,500 ms
コード長 1,424 bytes
コンパイル時間 4,286 ms
コンパイル使用メモリ 256,760 KB
最終ジャッジ日時 2025-02-19 19:46:17
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 39
権限があれば一括ダウンロードができます

ソースコード

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

#include <stdio.h>
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace atcoder;
using mint = modint998244353;
using namespace std;
#define rep(i,n) for (int i = 0; i < (n); ++i)
#define Inf32 1000000001
#define Inf64 1000000000000000001
int main(){
    vector<mint> kei = {-2,2,-8,2,-2,1,1,1,1,1,1,1,1};
vector<vector<int>> bunsi = {{2,0,2,0},{0,2,2,0},{1,1,1,1},{2,0,0,2},{0,2,0,2},{2,0,2,0},{0,2,2,0},{2,0,0,2},{0,2,0,2},{2,0,2,0},{2,0,0,2},{0,2,2,0}
    ,{0,2,0,2}};
    //cout<<kei.size()<<' '<<bunsi.size()<<endl;
    int _t;
    cin>>_t;
    
    rep(_,_t){
        int n;
        cin>>n;
        vector<mint> x(n),y(n);
        rep(i,n){
            int a,b;
            cin>>a>>b;
            x[i] = a,y[i] = b;
        }
        mint ans = 0;
        rep(k,kei.size()){
            mint sum = 0;
            rep(i,n){
                mint A = x[i] * x[i] + y[i] * y[i];
                {
                    mint t = 1;
                    t *= x[i].pow(bunsi[k][2]);
                    t *= y[i].pow(bunsi[k][3]);
                    if(k<5)t /= A;
                    else if(k>=9)t /= A*A;
                    ans += t * sum * kei[k];
                }
                {
                    mint t = 1;
                    t *= x[i].pow(bunsi[k][0]);
                    t *= y[i].pow(bunsi[k][1]);
                    if(k<5)t /= A;
                    else if(k<9)t /= A*A;
                    sum += t;
                }
                /*
                for(int j=i+1;j<n;j++){
                    mint t = 1;
                    
                    
                    
                    mint B = x[j] * x[j] + y[j] * y[j];
                    if(k<5){
                        t /= A*B;
                        
                    }
                    else if(k<9){
                        t /= A*A;
                        
                    }
                    else t /= B*B;
                    t *= kei[k];
                    ans += t;
                }
                */
            }
        }
        cout<<ans.val()<<endl;
    }
    
    return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0