Problem:
MathJax.Hub.Config({ tex2jax: { inlineMath: [["$","$"], ["\\(","\\)"]], processEscapes: true }});

静かな町

アイヅ国では、毎年、駅伝大会が開催されています。アイヅ国には N 個の町が点在しており、それぞれ 1 から N までの番号が付いています。いくつかの町の間は、互いに直接行き来できる道でつながっています。また、どの町の間もいくつかの道を辿って行き来ができます。大会のコースは、次のようにして決めます。

２つの町のすべての組み合わせについて、最短経路の距離を求める。

それらの中で、最短経路の距離が最大になるような２つの町を、スタートの町とゴールの町とする。町の組み合わせが複数ある場合、そのうちのどれか一つを選ぶ。

選ばれたスタートの町とゴールの町の間の最短経路を大会のコースとする。最短経路が複数ある場合、そのうちのどれか一つを選ぶ。

ヤマト国から来たお坊さんのトクイツは、アイヅ国のできるだけ静かな町に新しいお寺を開きたいと考えています。そのため、駅伝大会のコースに使われる可能性がない町を知りたがっています。

アイヅ国の町を結ぶ道の情報が与えられたとき、駅伝大会のコースに使われる可能性がない町を求めるプログラムを作成せよ。

Input

入力は以下の形式で与えられる。

N R
s1 t1 d1
s2 t2 d2
:
sR tR dR

１行目に町の数 N (2 ≤ N ≤ 1500) と、町の間を直接つなぐ道の数 R (1 ≤ R ≤ 3000) が与えられる。続く R 行に２つの町を双方向に直接つなぐ道が与えられる。si と ti (1 ≤ si < ti ≤ N) は i 番目の道がつなぐ２つの町の番号を表し、di (1 ≤ di ≤ 1000) はその道の長さを表す。ただし、どの２つの町についても、それらを直接つなぐ道は高々１本とする。

Output

与えられた道の情報に対して、駅伝大会のコースになることがない町をすべて出力する。１行目に駅伝大会のコースになることがない町の数 M を出力する。続く M 行に、そのような町の番号を昇順で出力する。

Sample Input 1

4 5
1 2 2
1 3 2
2 3 1
2 4 2
3 4 1

Sample Output 1

1
2

Sample Input 2

7 11
1 2 2
1 3 1
2 4 2
2 5 2
2 6 1
3 4 1
3 5 1
3 6 1
4 7 2
5 7 2
6 7 1

Sample Output 2

3
2
4
5

Solution in C++ (judge verdict: Time Limit Exceeded)

#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
typedef pair<int, int> P;

#define pb push_back
#define rep(i, a, n) for(int i = (a); i < (n); i++)
#define dep(i, a, n) for(int i = (a); i >= (n); i--)
#define mod 1e9+7

__attribute__((constructor))
void initial() {
  cin.tie(0);
  ios::sync_with_stdio(false);
}

struct edge {int to, cost;};

int n, r;
vector<edge> g[1501];
int d[1501][1501];
map<int, int> mp;

void dijkstra(int s) {
  priority_queue<P, vector<P>, greater<P> > que;
  fill(d[s], d[s] + n, 1145141919);
  d[s][s] = 0;
  que.push(P(0, s));

  while(!que.empty()) {
    P p = que.top(); que.pop();
    int v = p.second;
    if(d[s][v] < p.first) continue;
    rep(i, 0, g[v].size()) {
      edge e = g[v][i];
      if(d[s][e.to] > d[s][v] + e.cost) {
        d[s][e.to] = d[s][v] + e.cost;
        que.push(P(d[s][e.to], e.to));
      }
    }
  }
}

int solve(int i, int j) {
  while(1) {
    mp[j] = 1;
    if(j == i) break;
    rep(k, 0, g[j].size()) {
      edge e = g[j][k];
      if(d[i][j] == d[i][e.to] + e.cost) {
        j = e.to;
        break;
      }
    }
  }
}

int main() {
  cin >> n >> r;
  rep(i, 0, r) {
    int a, b, c;
    cin >> a >> b >> c;
    edge e;
    e.to = b, e.cost = c;
    g[a].pb(e);
    e.to = a;
    g[b].pb(e);
  }
  rep(i, 1, n + 1) {
    dijkstra(i);
  }
  int mas = 0;
  rep(i, 1, n + 1) {
    rep(j, 1, n + 1) {
      mas = max(mas, d[i][j]);
    }
  }
  rep(i, 1, n + 1) {
    rep(j, 1, n + 1) {
      if(mas == d[i][j]) {
        solve(i, j);
      }
    }
  }
  int count = 0;
  rep(i, 1, n + 1 ) {
    if(!mp[i]) count++;
  }
  cout << count << endl;
  rep(i, 1, n + 1) {
    if(!mp[i]) cout << i << endl;
  }
}