3 5
2 1 608591421 720962564
3 2 608591421 641991574
3 2 718110706 353409763
2 3 482786015 858341723
3 2 825397574 357507279
 2 1
2 1 317425595 636156322
 50000 1
20845 28516 904146168 997637433

Yes
641991574 1
2 
 Yes
0 0

 Yes
0 0


#include <bits/stdc++.h>
using namespace std;
int IN() {
  int c, f, x;
  while (!isdigit(c = getchar()) && c != '-')
    ;
  c == '-' ? (f = 1, x = 0) : (f = 0, x = c - '0');
  while (isdigit(c = getchar())) x = (x << 1) + (x << 3) + c - '0';
  return !f ? x : -x;
}
const int N = 500000;
vector<int> V;
int Ans[N];
namespace TwoSAT {
const int N = 500000;
const int M = 1000000;
struct Edge {
  int y, nxt;
} E[M];
vector<int> res;
int las[N], dfn[N], low[N], scc[N], stk[N], seq[N], vis[N];
int cnt, id, num, n;
void Init() {
  memset(las, -1, sizeof(las));
  memset(scc, 0, sizeof(scc));
  memset(dfn, 0, sizeof(dfn));
  cnt = 0;
  id = 0;
  num = 0;
  *seq = 0;
}
void Add_Edge(int x, int y) {
  E[cnt] = (Edge){y, las[x]};
  las[x] = cnt++;
}
void Tarjan(int x) {
  low[x] = dfn[x] = ++id;
  stk[++*stk] = x;
  for (int i = las[x], y; ~i; i = E[i].nxt)
    if (!dfn[y = E[i].y]) {
      Tarjan(y);
      low[x] = min(low[x], low[y]);
    } else {
      if (!scc[y]) low[x] = min(low[x], dfn[y]);
    }
  if (low[x] == dfn[x]) {
    int u;
    num++;
    do {
      u = stk[(*stk)--];
      seq[++*seq] = u;
      scc[u] = num;
    } while (u != x);
  }
}
bool Main(int __n) {
  n = __n;
  for (int i = 0; i < n; i++) {
    if (!dfn[i]) Tarjan(i);
  }
  for (int i = 0; i < n; i++) {
    if (scc[i] == scc[i ^ 1]) return 0;
  }
  return 1;
}
void dfs(int x) {
  if (~vis[x]) return;
  vis[x] = 1;
  for (int i = las[x]; ~i; i = E[i].nxt) dfs(E[i].y);
}
void getsol() {
  memset(vis, -1, sizeof(vis));
  for (int i = *seq, x; i; i--)
    if (vis[x = seq[i]] == -1) {
      vis[x] = 0;
      dfs(x ^ 1);
    }
  for (int i = 0; i < n; i++) {
    if (vis[2 * i] == 1) Ans[i] = 1;
  }
}
}  // namespace TwoSAT
struct Edge {
  int x, y, c, t;
} E[N];
map<int, vector<int> > M[N];
int A[N], id[N], pre[N];
int n, m, l, r, res = -1, num;
bool check(int T) {
  TwoSAT::Init();
  for (int i = 1; i < m + 1; i++) id[i] = 2 * (i - 1);
  for (int i = 1; i < m + 1; i++) {
    if (E[i].t > T) TwoSAT::Add_Edge(id[i], id[i] ^ 1);
  }
  num = 2 * m;
  for (int x = 1; x < n + 1; x++) {
    *A = 0;
    for (map<int, vector<int> >::iterator it = M[x].begin(); it != M[x].end();
         it++) {
      if (it->second.size() == 1) {
        A[++*A] = it->second[0];
      } else {
        int x = it->second[0], y = it->second[1];
        TwoSAT::Add_Edge(id[x] ^ 1, id[y]);
        TwoSAT::Add_Edge(id[y] ^ 1, id[x]);
      }
    }
    pre[0] = num;
    num += 2;
    for (int i = 1; i < *A + 1; i++) {
      int x = A[i];
      pre[i] = num;
      num += 2;
      TwoSAT::Add_Edge(pre[i - 1], pre[i]);
      TwoSAT::Add_Edge(pre[i] ^ 1, pre[i - 1] ^ 1);
      TwoSAT::Add_Edge(id[x], pre[i]);
      TwoSAT::Add_Edge(pre[i] ^ 1, id[x] ^ 1);
      TwoSAT::Add_Edge(pre[i - 1], id[x] ^ 1);
      TwoSAT::Add_Edge(id[x], pre[i - 1] ^ 1);
    }
  }
  return TwoSAT::Main(num);
}
int main() {
  n = IN(), m = IN();
  for (int i = 1; i < m + 1; i++) {
    E[i] = (Edge){IN(), IN(), IN(), IN()};
    M[E[i].x][E[i].c].push_back(i);
    M[E[i].y][E[i].c].push_back(i);
    if (M[E[i].x][E[i].c].size() > 2) return puts("No"), 0;
    if (M[E[i].y][E[i].c].size() > 2) return puts("No"), 0;
  }
  l = 0, r = 1e9;
  while (l <= r) {
    int mid = (l + r >> 1);
    if (check(mid))
      res = mid, r = mid - 1;
    else
      l = mid + 1;
  }
  if (~res) {
    puts("Yes");
    printf("%d ", res);
    check(res);
    TwoSAT::getsol();
    for (int i = 1; i < m + 1; i++) {
      if (Ans[i]) V.push_back(i);
    }
    printf("%d\n", V.size());
    for (int i = 0; i < V.size(); i++) printf("%d ", V[i] + 1);
    puts("");
  } else {
    puts("No");
  }
}
