結果

問題 No.1112 冥界の音楽
コンテスト
ユーザー fastmath
提出日時 2020-07-10 21:59:41
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 22 ms / 2,000 ms
コード長 2,104 bytes
コンパイル時間 1,600 ms
コンパイル使用メモリ 167,728 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-11 09:10:10
合計ジャッジ時間 2,967 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 34
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;
#define int long long
#define ii pair <int, int>
#define app push_back
#define all(a) a.begin(), a.end()
#define bp __builtin_popcountll
#define ll long long
#define mp make_pair
#define f first
#define s second
#define Time (double)clock()/CLOCKS_PER_SEC
#define debug(x) std::cout << #x << ": " << x << '\n';

const int MOD = 1000 * 1000 * 1000 + 7;

const int N =36;
void add(int &a, int b) {
    a += b;
    a %= MOD;
}   

void clear(int a[N][N], bool t) {
    for (int i = 0; i < N; ++i)
        for (int j = 0; j < N; ++j)
            a[i][j] = t * (i == j);
}   
void cop(int a[N][N], int b[N][N]) {
    for (int i = 0; i < N; ++i)
        for (int j = 0; j < N; ++j)
            a[i][j] = b[i][j];
}   
int t[N][N];
void add(int a[N][N], int b[N][N]) {
    clear(t, 0);
    for (int i = 0; i < N; ++i)
        for (int j = 0; j < N; ++j)
            for (int k = 0; k < N; ++k)
                add(t[i][j], a[i][k] * b[k][j]);
    cop(a, t);
}   
void print(int a[N][N]) {
    for (int i = 0; i < N; ++i) {
        for (int j = 0; j < N; ++j) 
            cout << a[i][j] << ' ';
        cout << '\n';
    }   
    cout << '\n';
}
int ans[N][N], cur[N][N];
void fp(int a[N][N], int p) {
    cop(cur, a);
    clear(ans, 1);
    for (int i = 0; (1ll << i) <= p; ++i) {
        if ((p >> i) & 1) 
            add(ans, cur);
        add(cur, cur);
    }       
    cop(a, ans);
}   

int a[N][N];

signed main() {
    #ifdef HOME
    freopen("input.txt", "r", stdin);
    #else
    #define endl '\n'
    ios_base::sync_with_stdio(0); cin.tie(0);
    #endif

    int k, m, n;
    cin >> k >> m >> n;
                        
    for (int i = 0; i < m; ++i) {
        int aa, b, c;
        cin >> aa >> b >> c;

        --aa; --b; --c;

        int u = k * aa + b;
        int v = k * b + c;

        a[u][v] = 1;
    }   

    fp(a, n - 2);

    int ans = 0;
    for (int i = 0; i < k * k; ++i)
        for (int j = 0; j < k * k; ++j) {
            if (i<k && j%k==0) {
                add(ans, a[i][j]);
            }
        }
    cout << ans << endl;
}
0