codeforces632F Magic Matrix -- bitset

最新推荐文章于 2019-05-23 18:01:00 发布

gjghfd

最新推荐文章于 2019-05-23 18:01:00 发布

阅读量380

点赞数

分类专栏： bitset 文章标签： codeforces bitset

本文链接：https://blog.csdn.net/gjghfd/article/details/74838978

版权

bitset 专栏收录该内容

3 篇文章 0 订阅

订阅专栏

每一行维护一个 $bitset$ 记录已经出现过的列，将每个点按从小到大的顺序插入，插入 $a_{i,j}$ 时判断第 $i$ 行与第 $j$ 行 $and$ 起来是否大于 $0$ 就可以了。
手写 $bitset$ 时间复杂度 $O({n^3 \over 64})$

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<bitset>
#include<vector>
using namespace std;
inline char nc(){
    static char buf[100000],*p1=buf,*p2=buf;
    if(p1==p2){
        p2=(p1=buf)+fread(buf,1,100000,stdin);
        if(p1==p2)return EOF;
    }
    return *p1++;
}
inline void Read(int& x){
    char c=nc();
    for(;c<'0'||c>'9';c=nc());
    for(x=0;c>='0'&&c<='9';x=(x<<3)+(x<<1)+c-48,c=nc());
}
#define N 2510
#define ll long long
struct Node{
    int x,y;
    Node(){}
    Node(int x,int y):x(x),y(y){}
};
struct Bitset{
    ll bit[40];
    void set(int x){
        x--;
        bit[x/63]|=1ll<<x%63;
    }
}f[N];
struct L{
    int w,f1,f2;
}c[N*N];
vector<Node>g[N*N];
int a[N][N],i,j,k,n,m,M,x,y,Len,t;
inline int Max(int x,int y){
    return x<y?y:x;
}
inline bool Check(int x,int y){
    for(int i=0;i<=Len;i++)
    if(f[x].bit[i]&f[y].bit[i])return 1;
    return 0;
}
inline bool Cmp(L a,L b){
    return a.w<b.w;
}
int main(){
    Read(n);
    Len=(n-1)/63;
    for(i=1;i<=n;i++)
    for(j=1;j<=n;j++)
    Read(a[i][j]),c[++t].w=a[i][j],c[t].f1=i,c[t].f2=j;;
    for(i=1;i<=n;i++){
        if(a[i][i]){
            printf("NOT MAGIC\n");
            return 0;
        }
        for(j=1;j<=n;j++)
        if(a[i][j]^a[j][i]){
            printf("NOT MAGIC\n");
            return 0;
        }
    }
    sort(c+1,c+t+1,Cmp);
    M=a[c[1].f1][c[1].f2]=0;
    for(i=2;i<=t;i++)
    if(c[i].w==c[i-1].w)a[c[i].f1][c[i].f2]=M;else a[c[i].f1][c[i].f2]=++M;
    for(i=1;i<=n;i++)
    for(j=1;j<=n;j++)
    g[a[i][j]].push_back(Node(i,j));
    for(i=0;i<=M;i++){
        for(j=0;j<g[i].size();j++)
        if(Check(g[i][j].x,g[i][j].y)){
            printf("NOT MAGIC\n");
            return 0;
        }
        for(j=0;j<g[i].size();j++){
            x=g[i][j].x;y=g[i][j].y;
            f[x].set(y);f[y].set(x);
        }
    }
    printf("MAGIC\n");
    return 0;
}