結果

問題 No.1050 Zero (Maximum)
ユーザー kenta255kenta255
提出日時 2020-05-09 00:58:44
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 23 ms / 2,000 ms
コード長 1,859 bytes
コンパイル時間 1,835 ms
コンパイル使用メモリ 206,140 KB
実行使用メモリ 4,708 KB
最終ジャッジ日時 2023-09-17 04:55:45
合計ジャッジ時間 2,727 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,380 KB
testcase_01 AC 1 ms
4,380 KB
testcase_02 AC 6 ms
4,376 KB
testcase_03 AC 4 ms
4,376 KB
testcase_04 AC 15 ms
4,384 KB
testcase_05 AC 16 ms
4,380 KB
testcase_06 AC 8 ms
4,380 KB
testcase_07 AC 9 ms
4,376 KB
testcase_08 AC 2 ms
4,380 KB
testcase_09 AC 5 ms
4,376 KB
testcase_10 AC 20 ms
4,708 KB
testcase_11 AC 14 ms
4,632 KB
testcase_12 AC 2 ms
4,380 KB
testcase_13 AC 2 ms
4,384 KB
testcase_14 AC 1 ms
4,380 KB
testcase_15 AC 1 ms
4,376 KB
testcase_16 AC 21 ms
4,688 KB
testcase_17 AC 23 ms
4,584 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
//#define P pair<ll,ll>
#define FOR(I,A,B) for(ll I = ll(A); I < ll(B); ++I)
#define FORR(I,A,B) for(ll I = ll((B)-1); I >= ll(A); --I)
#define TO(x,t,f) ((x)?(t):(f))
#define SORT(x) (sort(x.begin(),x.end())) // 0 2 2 3 4 5 8 9
#define POSL(x,v) (lower_bound(x.begin(),x.end(),v)-x.begin()) //xi>=v  x is sorted
#define POSU(x,v) (upper_bound(x.begin(),x.end(),v)-x.begin()) //xi>v  x is sorted
#define NUM(x,v) (POSU(x,v)-POSL(x,v))  //x is sorted
#define REV(x) (reverse(x.begin(),x.end())) //reverse
ll gcd(ll a,ll b){if(a%b==0)return b;return gcd(b,a%b);}
ll lcm(ll a,ll b){ll c=gcd(a,b);return ((a/c)*(b/c)*c);}
#define NEXTP(x) next_permutation(x.begin(),x.end())
const ll INF=ll(1e16)+ll(7);
const ll MOD=1000000007LL;
#define out(a) cout<<fixed<<setprecision((a))
//tie(a,b,c) = make_tuple(10,9,87);
#define pop(a) __builtin_popcount((a))
ll keta(ll a){ll r=0;while(a){a/=10;r++;}return r;}



ll cnt[52][52];
ll B[52][52];
ll M;
vector< vector<ll> > pow_(vector<vector<ll>> A,ll k){
  vector< vector<ll> > res(M,vector<ll>(M,0));
  if(k==0){
    FOR(i,0,M)res[i][i] = 1;
    return res;
  }
  if(k==1){
    return A;
  }
  auto x = pow_(A,k/2);

  FOR(i,0,M){
    FOR(j,0,M){
      FOR(m,0,M){
        (res[i][j] += x[i][m]*x[m][j]) %=MOD;
      }
    }
  }

  if(k&1){
    FOR(i,0,M)FOR(j,0,M)x[i][j] = 0;
    FOR(i,0,M){
      FOR(j,0,M){
        FOR(m,0,M){
          (x[i][j] += res[i][m]*B[m][j]) %=MOD;
        }
      }
    }
    return x;
  }
  return res;
}




int main(){


  ll K;
  cin >> M >> K;

  FOR(i,0,M){
    FOR(j,0,M){
      cnt[i][(i*j)%M]++;
      cnt[i][(i+j)%M]++;
    }
  }

  vector< vector<ll> > A(M,vector<ll>(M,0));
  FOR(i,0,M)FOR(j,0,M){
    A[i][j] = cnt[j][i];
    B[i][j] = cnt[j][i];
  }

  auto ans = pow_(A,K);
  cout << ans[0][0] << endl;
}
0