結果

問題 No.1136 Four Points Tour
ユーザー leafirby
提出日時 2020-07-19 00:14:49
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 12 ms / 2,000 ms
コード長 1,957 bytes
コンパイル時間 2,109 ms
コンパイル使用メモリ 205,580 KB
最終ジャッジ日時 2025-01-12 00:38:11
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 41
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
/*#include <boost/multiprecision/cpp_dec_float.hpp>
#include <boost/multiprecision/cpp_int.hpp>
*/#pragma GCC target("avx")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")    
using namespace std;
/*namespace mp = boost::multiprecision;
using Bint = mp::cpp_int;
using Real = mp::number<mp::cpp_dec_float<1024>>;
*/#define lli long long int
#define uli unsigned long long int
#define INF 999999999999
#define rep(i,m,n) for(lli i = m;i < n;i++)
#define rrep(i,m,n) for(lli i=m-1;i>=n;i--)
#define pb(n) push_back(n)
#define UE(N) N.erase(unique(N.begin(),N.end()),N.end());
#define Sort(n) sort(n.begin(), n.end())
#define Rev(n) reverse(n.begin(),n.end())
#define Out(S) cout << S << endl
#define NeOut(S) cout << S
#define HpOut(S) cout << setprecision(50) << S << endl
#define Vec(K,L,N,S) vector<L> K(N,S)
#define DV(K,L,N,M,S) vector<vector<L>> K(N,vector<L>(M,S))
#define TV(K,L,N,M,R,S) vector<vector<vector<L>>> K(N,vector<vector<L>>(M,vector<L>(R,S)))
#define pint pair<lli,lli>
#define paf(L,R) pair<L,R>
#define mod 1000000007
#define MAX 10000000
#define ALL(a)  a.begin(),a.end()
#define chmax(a, b) a = (((a)<(b)) ? (b) : (a))
#define chmin(a, b) a = (((a)>(b)) ? (b) : (a))
vector<lli> mtmul(vector<lli> &dp,vector<vector<lli>> &mt){
  int m=mt.size();
  Vec(ret,lli,m,0);
  rep(i,0,m)rep(j,0,m)ret[i]=(ret[i]+(mt[i][j]*dp[j])%mod)%mod;
  return ret;
}
vector<vector<lli>> update(vector<vector<lli>> &mt){
  int m=mt.size();
  DV(ret,lli,m,m,0);
  rep(i,0,m)rep(j,0,m)rep(k,0,m)ret[i][j]=(ret[i][j]+(mt[i][k]*mt[k][j])%mod)%mod;
  return ret;
}
int main(){
  ios::sync_with_stdio(false);
  cin.tie(0);
  lli A,B,C,H,W,M,K,L,R,N,num=0,sum=0,flag=0;string S,T;
  cin >> N;
  vector<vector<lli>> mt1(2);
  vector<lli> mt2(2);
  mt1[0].pb(-1);
  mt1[0].pb(1);
  mt1[1].pb(0);
  mt1[1].pb(3);
  mt2[0]=mt2[1]=1;
  while(N){
    if(N&1)mt2=mtmul(mt2,mt1);
    mt1=update(mt1);
    N/=2;
  }
  Out(mt2[0]);
}
0