标签：tmp floyd ons 最小 一个 开始 去除 style gre
给你n个数字，对于ai和aj，如果ai&aj不等于0，那么就在i点和j点之间连上一条边，问你图中的最小环长度，没有环输出-1
最开始看n=1e5，觉得基本跟floyd绝缘，但是jxc给了一个牛逼结论：
如果n足够大，那么一定存在三个点，它们互相有边相连，最小环长度为3
进一步说，n有多大才满足条件呢？
考虑ai的二进制形式，由于ai<=1e18，最长不会超过64位
如果有至少64*2+1个非零数字，不管你如何安排，至少有一位上会有三个1，那么他们三者就可以构成最小环。
这样思路就明确了，只要去除0的影响（任何数字&0=0，这意味着这个点不会有边相连），再考虑剩下的数字，大于200直接输出3，小的时候n方连边再判环
注意一个细节，位运算优先级比不等于还低，所以判断有没有边的时候记得加括号
代码：
#include <bits/stdc++.h> #define int long long #define sc(a) scanf("%lld",&a) #define scc(a,b) scanf("%lld %lld",&a,&b) #define sccc(a,b,c) scanf("%lld %lld %lld",&a,&b,&c) #define schar(a) scanf("%c",&a) #define pr(a) printf("%lld",a) #define fo(i,a,b) for(int i=a;i<b;++i) #define re(i,a,b) for(int i=a;i<=b;++i) #define rfo(i,a,b) for(int i=a;i>b;--i) #define rre(i,a,b) for(int i=a;i>=b;--i) #define prn() printf("\n") #define prs() printf(" ") #define mkp make_pair #define pii pair<int,int> #define pub(a) push_back(a) #define pob() pop_back() #define puf(a) push_front(a) #define pof() pop_front() #define fst first #define snd second #define frt front() #define bak back() #define mem0(a) memset(a,0,sizeof(a)) #define memmx(a) memset(a,0x3f3f,sizeof(a)) #define memmn(a) memset(a,-0x3f3f,sizeof(a)) #define debug #define db double #define yyes cout<<"YES"<<endl; #define nno cout<<"NO"<<endl; using namespace std; typedef vector<int> vei; typedef vector<pii> vep; typedef map<int,int> mpii; typedef map<char,int> mpci; typedef map<string,int> mpsi; typedef deque<int> deqi; typedef deque<char> deqc; typedef priority_queue<int> mxpq; typedef priority_queue<int,vector<int>,greater<int> > mnpq; typedef priority_queue<pii> mxpqii; typedef priority_queue<pii,vector<pii>,greater<pii> > mnpqii; const int maxn=500005; const int inf=0x3f3f3f3f3f3f3f3f; const int MOD=100000007; const db eps=1e-10; int qpow(int a,int b){int tmp=a%MOD,ans=1;while(b){if(b&1){ans*=tmp,ans%=MOD;}tmp*=tmp,tmp%=MOD,b>>=1;}return ans;} int lowbit(int x){return x&-x;} int max(int a,int b){return a>b?a:b;} int min(int a,int b){return a<b?a:b;} int mmax(int a,int b,int c){return max(a,max(b,c));} int mmin(int a,int b,int c){return min(a,min(b,c));} void mod(int &a){a+=MOD;a%=MOD;} bool chk(int now){} int half(int l,int r){while(l<=r){int m=(l+r)/2;if(chk(m))r=m-1;else l=m+1;}return l;} int ll(int p){return p<<1;} int rr(int p){return p<<1|1;} int mm(int l,int r){return (l+r)/2;} int lg(int x){if(x==0) return 1;return (int)log2(x)+1;} bool smleql(db a,db b){if(a<b||fabs(a-b)<=eps)return true;return false;} db len(db a,db b,db c,db d){return sqrt((a-c)*(a-c)+(b-d)*(b-d));} bool isp(int x){if(x==1)return false;if(x==2)return true;for(int i=2;i*i<=x;++i)if(x%i==0)return false;return true;} int n,a,d[205][205],e[205][205]; vei v; int floyd(){ int mn=inf; re(k,1,n){ re(i,1,k-1){ re(j,i+1,k-1){ mn=min(mn,d[i][j]+e[j][k]+e[k][i]); } } re(i,1,n) re(j,1,n) d[i][j]=min(d[i][j],d[i][k]+d[k][j]); } return mn; } signed main(){ ios_base::sync_with_stdio(0); cin.tie(0),cout.tie(0); cin>>n; re(i,1,n){ cin>>a; if(a!=0) v.pub(a); } if(v.size()>=200) cout<<3; else if(v.size()==0) cout<<-1; else{ n=v.size(); re(i,1,n) re(j,1,n) d[i][j]=e[i][j]=1000000000; fo(i,0,v.size()){ fo(j,0,v.size()){ if(i!=j&&((v[i]&v[j])!=0)) d[i+1][j+1]=e[i+1][j+1]=1; } } // re(i,1,n) re(j,1,n) cout<<e[i][j]<<" \n"[j==n]; int ans=floyd(); if(ans>=1000000000) cout<<-1; else cout<<ans; } return 0; }
1205B Shortest Cycle（floyd找最小环）
标签：tmp floyd ons 最小 一个 开始 去除 style gre
原文地址：https://www.cnblogs.com/oneman233/p/11478800.html
