結果
問題 | No.613 Solitude by the window |
ユーザー |
![]() |
提出日時 | 2017-11-22 01:40:38 |
言語 | Java (openjdk 23) |
結果 |
AC
|
実行時間 | 141 ms / 2,000 ms |
コード長 | 1,198 bytes |
コンパイル時間 | 2,436 ms |
コンパイル使用メモリ | 78,136 KB |
実行使用メモリ | 54,440 KB |
最終ジャッジ日時 | 2024-12-14 06:09:34 |
合計ジャッジ時間 | 6,717 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 21 |
ソースコード
import java.util.*;class Main {static void myAssert(boolean e){if(!e)throw new RuntimeException("assertion sippai");}static boolean isPrime(long x){if(x<=1)return false;for(long i=2;i*i<=x;++i)if(x%i==0)return false;return true;}static long powerMod(long x, long exponent,long m) {long prod = 1;for (int i = 63; i >= 0; --i) {prod = (prod * prod) % m;if ((exponent & 1L << i) != 0) {prod = (prod * x) % m;}}return prod;}public static void main(String[] args) {Scanner scan = new Scanner(System.in);long n=scan.nextLong();long m=scan.nextLong();myAssert(n>=1);myAssert(n<=1000000000000000000L);myAssert(m>=2);myAssert(m<=1100000000);myAssert(isPrime(m));long p=m==2||m%12==1||m%12==11?m-1:m+1;long e=powerMod(2,n,p);long x=1,y=0;for(int i=31;i>=0;--i){long z=(x*x+3*y*y)%m,w=2*x*y%m;x=z;y=w;if((e&1L<<i)!=0){z=(2*x+3*y)%m;w=(x+2*y)%m;x=z;y=w;}}System.out.println((2*x+m-2)%m);}}