結果
問題 | No.1448 和差算 |
ユーザー |
👑 ![]() |
提出日時 | 2021-03-31 14:17:43 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 3 ms / 2,000 ms |
コード長 | 751 bytes |
コンパイル時間 | 1,976 ms |
コンパイル使用メモリ | 191,852 KB |
最終ジャッジ日時 | 2025-01-20 01:25:37 |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 36 |
ソースコード
#include <bits/stdc++.h>using namespace std;using ll=long long;using ull=unsigned long long;#define rep(i,n) for(int i=0; i<(n); i++)const ull M=1000000007;ull powm(ull a,ull i){if(i==0) return 1;ull r=powm(a*a%M,i/2);if(i%2==1) r=r*a%M;return r;}ll A,B,C,D,N;int main(){cin>>A>>B>>C>>D>>N;A=(A+M*100)%M;B=(B+M*100)%M;C=(C+M*100)%M;D=(D+M*100)%M;ull Mul = powm(16,N/8); N%=8;if(N%8>=4){Mul=Mul*4%M;swap(A,B); A=(M-A)%M; B=(M-B)%M;swap(C,D); C=(M-C)%M; D=(M-D)%M;N-=4;}if(N%4>=2){Mul=Mul*2%M;swap(A,C); swap(B,D);swap(A,B); A=(M-A)%M; B=(M-B)%M;N-=2;}if(N%2==1){cout<<(B*2*Mul%M)<<"\n";}else{cout<<((B+D)*Mul%M)<<"\n";}return 0;}