結果

問題 No.1105 Many Triplets
ユーザー hiro71687khiro71687k
提出日時 2023-03-10 20:01:00
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 3 ms / 2,000 ms
コード長 1,241 bytes
コンパイル時間 4,676 ms
コンパイル使用メモリ 253,300 KB
最終ジャッジ日時 2025-02-11 07:30:12
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 25
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;
using ll=long long;
using ld=long double;
ld pie=3.141592653589793;
ll inf=155550555555555;
ll mod=1000000007;
ll modpow(ll x, ll n) {
  x = x%mod;
  if(n==0) return 1;  //再帰の終了条件

  else if(n%2==1) {
    return (x*modpow(x, n-1))%mod;  //nが奇数ならnを1ずらす
  }
  else return modpow((x*x)%mod, n/2)%mod;  //nが偶数ならnが半分になる
}
int main(){
    ll n;
    cin >> n;
    ll a,b,c;
    cin >> a >> b >> c;
    n-=2;
    vector<tuple<ll,ll,ll>>t(6);
    t[0]={a-b,b-c,c-a};
    t[1]={a-2*b+c,a+b-2*c,-2*a+b+c};
    t[2]={-3*b+3*c,3*a-3*c,3*b-3*a};
    t[3]={-3*a-3*b+6*c,6*a-3*b-3*c,-3*a+6*b-3*c};
    t[4]={-9*a+9*c,9*a-9*b,9*b-9*c};
    t[5]={-18*a+9*b+9*c,9*a-18*b+9*c,9*a+9*b-18*c};
    for (ll i = 0; i < t.size(); i++)
    {
        get<0>(t[i])%=mod;
        get<1>(t[i])%=mod;
        get<2>(t[i])%=mod;
    }
    ll z=n/6;
    ll x=modpow(-27,z);
    x+=mod;
    x%=mod;
    n%=6;
    z=n;
    ll aa=(get<0>(t[z])*x)%mod,bb=(get<1>(t[z])*x)%mod ,cc=(get<2>(t[z])*x)%mod;
    aa+=mod;
    aa%=mod;
    bb+=mod;
    bb%=mod;
    cc+=mod;
    cc%=mod;
    cout << aa << ' ' << bb << ' ' << cc << endl;
}
0