結果
| 問題 | 
                            No.1340 おーじ君をさがせ
                             | 
                    
| コンテスト | |
| ユーザー | 
                             maspy
                         | 
                    
| 提出日時 | 2020-09-26 16:14:02 | 
| 言語 | C++14  (gcc 13.3.0 + boost 1.87.0)  | 
                    
| 結果 | 
                             
                                WA
                                 
                             
                            
                         | 
                    
| 実行時間 | - | 
| コード長 | 1,476 bytes | 
| コンパイル時間 | 861 ms | 
| コンパイル使用メモリ | 77,368 KB | 
| 実行使用メモリ | 15,244 KB | 
| 最終ジャッジ日時 | 2024-11-26 19:23:18 | 
| 合計ジャッジ時間 | 3,581 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge3 / judge1 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 2 WA * 1 | 
| other | AC * 43 WA * 13 RE * 3 | 
ソースコード
// 適当な時刻まで探索して、
// 最後 2 回の到着時刻で周期を決め打ってしまう
// 厳密解の意味では N^3M 解法?(自信なし)
#include<iostream>
#include<algorithm>
#include<vector>
using namespace std;
using ll = long long;
bool dp[1000000][100];
int frm[10000];
int to[10000];
ll N, M, T, ans;
int main(){
    cin >> N >> M >> T;
    vector<int> vec;
    for(int i=0;i<M;i++){
        int a, b;
        cin >> a >> b;
        vec.emplace_back(N*a+b);
    }
    sort(vec.begin(), vec.end());
    vec.erase(unique(vec.begin(), vec.end()), vec.end());
    M = vec.size();
    for(int i=0;i<M;i++){
        frm[i] = vec[i] / N;
        to[i] = vec[i] % N;
    }
    dp[0][0] = true;
    // 適当なところまで探索する
    int T0 = min(100000ll, 100000000/(N*M));
    for(int t=1;t<=T0;t++){
        for(int e=0;e<M;e++){
            dp[t][to[e]] |= dp[t-1][frm[e]];
        }
    }
    for(int v=0;v<N;v++){
        // last two
        ll t1 = -1;
        ll t2 = -1;
        for (int t=0;t<T0;t++){
            if(dp[t][v]){
                t1 = t2;
                t2 = t;
            }
        }
        if(T < T0){
            ans += (dp[T][v] ? 1 : 0);
        } else{
            if(t1 >= 0){
                ll p = t2 - t1;
                ll n = (T-T0)/p+1;
                ans += (dp[T-n*p][v] ? 1 : 0);
            }
        }
    }
    if(ans==0){
        ans-=1;
    }
    cout << ans << endl;
    return 0;
}
            
            
            
        
            
maspy