結果

問題 No.1253 雀見椪
コンテスト
ユーザー shiomusubi496
提出日時 2020-10-09 23:45:31
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 105 ms / 2,000 ms
コード長 958 bytes
コンパイル時間 1,550 ms
コンパイル使用メモリ 167,316 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-07-20 14:35:46
合計ジャッジ時間 3,376 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 14
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int mod=1000000007;
int gcd(int a,int b){
  if(!b)return a;
  return gcd(b,a%b);
}
int lcm(int a,int b){
  return a/gcd(a,b)*b;
}
int modpow(int a,int b){
  if(!b)return 1;
  if(b&1)return modpow(a,b-1)*a%mod;
  int c=modpow(a,b/2);
  return c*c%mod;
}
signed main(){
  int Q;
  cin>>Q;
  while(Q--){
    int N,a,b,c,d,e,f;
    cin>>N>>a>>d>>b>>e>>c>>f;
    int n=gcd(a,d);a/=n,d/=n;
    n=gcd(b,e);b/=n,e/=n;
    n=gcd(c,f);c/=n,f/=n;
    n=lcm(d,lcm(e,f));
    a*=n/d,b*=n/e,c*=n/f;
    int m=(modpow(n,N)-(modpow(a+b,N-1)-modpow(a,N-1)+mod)%mod*a%mod
      -(modpow(a+c,N-1)-modpow(a,N-1)+mod)%mod*a%mod
      -(modpow(b+a,N-1)-modpow(b,N-1)+mod)%mod*b%mod
      -(modpow(b+c,N-1)-modpow(b,N-1)+mod)%mod*b%mod
      -(modpow(c+a,N-1)-modpow(c,N-1)+mod)%mod*c%mod
      -(modpow(c+b,N-1)-modpow(c,N-1)+mod)%mod*c%mod+mod*6)%mod;
    cout<<m*modpow(modpow(n,N),mod-2)%mod<<endl;
  }
}
0