#include <algorithm>
#include <iostream>
#include <climits>
#include <cstdio>
#include <queue>

using namespace std;

typedef unsigned int uint;

const int INF = INT_MAX >> 1;
const int N = 505;
const int W = 32;
const int B = N / W + 2;
const int M = 200005;

typedef pair<int , int> data;

#define mkp(a , b) make_pair(a , b)
#define ft first
#define sd second

queue<data> q;

uint done[N][B] , graph[N][B];
int tov[M] , nxt[M];
bool vis[N];
int f[N][N];
int last[N];
int n , m , tot , cirlen , ans;

inline void set(uint *bits , int x){bits[x / W] |= 1 << x % W;}

inline void insert(int x , int y){tov[++ tot] = y , nxt[tot] = last[x] , last[x] = tot;}

int dfs(int x , int dep = 0)
{
	vis[x] = 1;
	for (int i = last[x] , y , ret;i;i = nxt[i])
	{
		if ((y = tov[i]) == 1) return dep + 1;
		if (!vis[y])
		{
			ret = dfs(y , dep + 1);
			if (ret) return ret;
		}
	}
	return 0;
}

void bfs()
{
	ans = INF;
	if (!cirlen) return;
	for (int i = 0;i < cirlen;++ i)
		for (int x = 1;x <= n;++ x)
			f[i][x] = INF;
	f[0][1] = 0 , q.push(mkp(0 , 1)) , set(done[0] , 0);
	for (data cur;!q.empty();)
	{
		cur = q.front() , q.pop();
		for (int i = 0 , nxti = cur.ft == cirlen - 1 ? 0 : cur.ft + 1;i <= (n - 1) / W;++ i)
		{
			uint state = ~done[nxti][i] & graph[cur.sd][i];
			if (state)
				for (int j = 0 , y = i * W + 1;j < W;++ j , ++ y)
					if ((state >> j) & 1)
					{
						f[nxti][y] = f[cur.ft][cur.sd] + 1;
						q.push(mkp(nxti , y)) , set(done[nxti], y - 1);
					}
		}
	}
	for (int i = 0 , ret;i < cirlen;++ i)
	{
		ret = 0;
		for (int x = 1;x <= n;++ x) ret = max(ret , f[i][x]);
		ans = min(ans , ret);
	}
}

int main()
{
	//freopen("f.in" , "r" , stdin) , freopen("f.out" , "w" , stdout);
	scanf("%d%d" , &n , &m);
	for (int i = 1 , x , y;i <= m;++ i) scanf("%d%d" , &x , &y) , insert(x , y) , set(graph[x] , y - 1);
	cirlen = dfs(1) , bfs() , printf("%d\n" , ans == INF ? -1 : ans);
	//fclose(stdin) , fclose(stdout);
	return 0;
}

### Prompt

In cpp, your task is to solve the following problem:

Mr. Endo wanted to write the code that performs breadth-first search (BFS), which is a search algorithm to explore all vertices on a directed graph. An example of pseudo code of BFS is as follows:


1: $current \leftarrow \{start\_vertex\}$
2: $visited \leftarrow current$
3: while $visited \ne$ the set of all the vertices
4:   $found \leftarrow \{\}$
5:   for $u$ in $current$
6:     for each $v$ such that there is an edge from $u$ to $v$
7:       $found \leftarrow found \cup \{v\}$
8:   $current \leftarrow found \setminus visited$
9:   $visited \leftarrow visited \cup found$


However, Mr. Endo apparently forgot to manage visited vertices in his code. More precisely, he wrote the following code:


1: $current \leftarrow \{start\_vertex\}$
2: while $current \ne$ the set of all the vertices
3:   $found \leftarrow \{\}$
4:   for $u$ in $current$
5:     for each $v$ such that there is an edge from $u$ to $v$
6:       $found \leftarrow found \cup \{v\}$
7:   $current \leftarrow found$


You may notice that for some graphs, Mr. Endo's program will not stop because it keeps running infinitely. Notice that it does not necessarily mean the program cannot explore all the vertices within finite steps. Your task here is to make a program that determines whether Mr. Endo's program will stop within finite steps for a given directed graph in order to point out the bug to him. Also, calculate the minimum number of loop iterations required for the program to stop if it is finite. Since the answer might be huge, thus print the answer modulo $10^9 +7$, which is a prime number.



Input

The input consists of a single test case formatted as follows.


$N$ $M$
$u_1$ $v_1$
:
$u_M$ $v_M$


The first line consists of two integers $N$ ($2 \leq N \leq 500$) and $M$ ($1 \leq M \leq 200,000$), where $N$ is the number of vertices and $M$ is the number of edges in a given directed graph, respectively. The $i$-th line of the following $M$ lines consists of two integers $u_i$ and $v_i$ ($1 \leq u_i, v_i \leq N$), which means there is an edge from $u_i$ to $v_i$ in the given graph. The vertex $1$ is the start vertex, i.e. $start\\_vertex$ in the pseudo codes. You can assume that the given graph also meets the following conditions.

* The graph has no self-loop, i.e., $u_i \ne v_i$ for all $1 \leq i \leq M$.
* The graph has no multi-edge, i.e., $(u_i, v_i) \le (u_j, v_j)$ for all $1 \leq i < j \leq M$.
* For each vertex $v$, there is at least one path from the start vertex $1$ to $v$.

Output

If Mr. Endo's wrong BFS code cannot stop within finite steps for the given input directed graph, print '-1' in a line. Otherwise, print the minimum number of loop iterations required to stop modulo $10^9+7$.

Examples

Input

4 4
1 2
2 3
3 4
4 1


Output

-1


Input

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


Output

7


Input

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


Output

3

### Response

```cpp
#include <algorithm>
#include <iostream>
#include <climits>
#include <cstdio>
#include <queue>

using namespace std;

typedef unsigned int uint;

const int INF = INT_MAX >> 1;
const int N = 505;
const int W = 32;
const int B = N / W + 2;
const int M = 200005;

typedef pair<int , int> data;

#define mkp(a , b) make_pair(a , b)
#define ft first
#define sd second

queue<data> q;

uint done[N][B] , graph[N][B];
int tov[M] , nxt[M];
bool vis[N];
int f[N][N];
int last[N];
int n , m , tot , cirlen , ans;

inline void set(uint *bits , int x){bits[x / W] |= 1 << x % W;}

inline void insert(int x , int y){tov[++ tot] = y , nxt[tot] = last[x] , last[x] = tot;}

int dfs(int x , int dep = 0)
{
	vis[x] = 1;
	for (int i = last[x] , y , ret;i;i = nxt[i])
	{
		if ((y = tov[i]) == 1) return dep + 1;
		if (!vis[y])
		{
			ret = dfs(y , dep + 1);
			if (ret) return ret;
		}
	}
	return 0;
}

void bfs()
{
	ans = INF;
	if (!cirlen) return;
	for (int i = 0;i < cirlen;++ i)
		for (int x = 1;x <= n;++ x)
			f[i][x] = INF;
	f[0][1] = 0 , q.push(mkp(0 , 1)) , set(done[0] , 0);
	for (data cur;!q.empty();)
	{
		cur = q.front() , q.pop();
		for (int i = 0 , nxti = cur.ft == cirlen - 1 ? 0 : cur.ft + 1;i <= (n - 1) / W;++ i)
		{
			uint state = ~done[nxti][i] & graph[cur.sd][i];
			if (state)
				for (int j = 0 , y = i * W + 1;j < W;++ j , ++ y)
					if ((state >> j) & 1)
					{
						f[nxti][y] = f[cur.ft][cur.sd] + 1;
						q.push(mkp(nxti , y)) , set(done[nxti], y - 1);
					}
		}
	}
	for (int i = 0 , ret;i < cirlen;++ i)
	{
		ret = 0;
		for (int x = 1;x <= n;++ x) ret = max(ret , f[i][x]);
		ans = min(ans , ret);
	}
}

int main()
{
	//freopen("f.in" , "r" , stdin) , freopen("f.out" , "w" , stdout);
	scanf("%d%d" , &n , &m);
	for (int i = 1 , x , y;i <= m;++ i) scanf("%d%d" , &x , &y) , insert(x , y) , set(graph[x] , y - 1);
	cirlen = dfs(1) , bfs() , printf("%d\n" , ans == INF ? -1 : ans);
	//fclose(stdin) , fclose(stdout);
	return 0;
}

```