最小环-Floyd

floyd求最小环

在Floyd的同时，顺便算出最小环。

Floyd算法

for(k=1；k<=n；k++)
  { for(i=1；i<k；i++)
     for(j=i+1；j<k；j++)
      if(d[i][j]+m[i][k]+m[k][j]<min)
       min=d[i][j]+m[i][k]+m[k][j]；
    for(i=1；i<=n；i++)
     for(j=1；j<=n；j++)
      if(d[i][k]+d[k][j]<d[i][j])
       d[i][j]=d[i][k]+d[k][j]；
  }

 

保证了最外层循环到 k 时所有顶点间已求得以 0...k-1 为中间点的最短路径。一
个环至少有 3 个顶点,设某环编号最大的顶点为 L ,在环中直接与之相连的两个顶点编号
分别为 M 和 N (M,N < L),则最大编号为 L 的最小环长度即为 Graph(M,L) + Graph(N,L) +
Dist(M,N) ,其中 Dist(M,N) 表示以 0...L-1 号顶点为中间点时的最短路径,刚好符合 Floyd
算法最外层循环到 k=L 时的情况,则此时对 M 和 N 循环所有编号小于 L 的顶点组合即
可找到最大编号为 L 的最小环。再经过最外层 k 的循环,即可找到整个图的最小环。

上面是对无向图的情况

若是有向图,只需稍作改动。注意考虑有向图中 2 顶点即可组成环的情况

 

题目：

D. Shortest Cycle

time limit per test

1 second

memory limit per test

256 megabytes

input

standard input

output

standard output

You are given nn integer numbers a1,a2,…,ana1,a2,…,an. Consider graph on nn nodes, in which nodes ii, jj (i≠ji≠j) are connected if and only if, aiaiAND aj≠0aj≠0, where AND denotes the bitwise AND operation.

Find the length of the shortest cycle in this graph or determine that it doesn't have cycles at all.

Input

The first line contains one integer nn (1≤n≤105)(1≤n≤105) — number of numbers.

The second line contains nn integer numbers a1,a2,…,ana1,a2,…,an (0≤ai≤10180≤ai≤1018).

Output

If the graph doesn't have any cycles, output −1−1. Else output the length of the shortest cycle.

Examples

input

Copy

4
3 6 28 9

output

Copy

4

input

Copy

5
5 12 9 16 48

output

Copy

3

input

Copy

4
1 2 4 8

output

Copy

-1

Note

In the first example, the shortest cycle is (9,3,6,28)(9,3,6,28).

In the second example, the shortest cycle is (5,12,9)(5,12,9).

The graph has no cycles in the third example.

分析：对于大于2 * 64个正数的情况直接输出3；其余的情况怼入vector跑floyd求最小环。

代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int inf = 0x3f3f3f3f;
const int maxn = 1e5 + 10;
vector<ll> vec;
ll g[210][210];
ll dis[210][210];
int num;

int floyd()
{
    ll res = 1e9;
    for (int i = 1; i <= num; i++)
        for (int j = 1; j <= num; j++)
            dis[i][j] = g[i][j];
    for (int k = 1; k <= num; k++)
    {
        for (int i = 1; i < k; i++)
            for (int j = i + 1; j < k; j++)
                res = min(res, dis[i][j] + g[i][k] + g[k][j]);
        for (int i = 1; i <= num; i++)
            for (int j = 1; j <= num; j++)
                dis[i][j] = min(dis[i][j], dis[i][k] + dis[k][j]);
    }
    return res != 1e9 ? res : -1;
}

int main()
{
    int n; cin >> n;
    int cnt = 0;
    ll x;

    for (int i = 1; i <= n; i++)
    {
        scanf("%lld", &x);
        if (x > 0) cnt++, vec.push_back(x);
    }
    if (cnt >= 200)
    {
        cout << "3" << endl;
        return 0;
    }
    num = vec.size();
    for (int i = 0; i < num; i++)
        for (int j = i + 1; j < num; j++)
            g[i + 1][j + 1] = g[j + 1][i + 1] = (vec[i] & vec[j]) ? 1 : inf;
    cout << floyd() << endl;
}