8.10上课感悟

今天上课讲图论，感觉还行
2-sat问题我还没写，有空的时候我再写写吧
回家写了2题，感觉还行

poj 3613

题目传送门
上课已经讲过啦，就是简单的“矩阵快速幂”+floyd
考虑floyd的写法 f[i][j]=min{f[i][k]+f[k][j]}
和矩阵的 s[i][j]+=a[i][k]*b[k][j]是不是很像
只是把求和改成求Min,把乘改成加
我们暂时可以默认它有这个结论，只要跑出来能AC就行
实现可能比较重要，那就发一下代码吧

void floyd(LL c[][N],LL a[][N],LL b[][N]){
    for(int k=1;k<=cnt;k++)
        for(int i=1;i<=cnt;i++)
            for(int j=1;j<=cnt;j++)
                c[i][j]=min(c[i][j],a[i][k]+b[k][j]);
}
void copy(LL a[][N],LL b[][N]){
    for(int i=1;i<=cnt;i++)
        for(int j=1;j<=cnt;j++)
            a[i][j]=b[i][j],b[i][j]=inf;
}
void calc(int k){
    while(k){
        if(k&1){
            floyd(tmp,ans,a);
            copy(ans,tmp);//一开始我直接写floyd(ans,ans,a)像简短点代码，发现很有问题
        }
        floyd(tmp,a,a);copy(a,tmp);
        k>>=1;
    }
}

主要是求解的地方可能会写错

bzoj 1059

题目传送门
这题上课也讲过啦
把行和列看作是二分图两边的点，然后如果有黑点就连边
最后看最大匹配是不是等于N即可
代码还是发一下吧，怕一些人不会写匈牙利

#include<bits/stdc++.h>
#define N 210
using namespace std;
int n,vis[N],ans[N],f[N][N];
bool dfs(int x){
    for(int i=1;i<=n;i++)
        if(f[x][i]&&!vis[i]){
            vis[i]=1;
            if(!ans[i]||dfs(ans[i])){
                ans[i]=x;
                return true;
            }
        }
    return false;
}
int main(){
    int T;
    scanf("%d",&T);
    while(T--){
        scanf("%d",&n);
        memset(f,0,sizeof(f));
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++){
                int col;scanf("%d",&col);
                if(col) f[i][j]=1;
            }
        int cnt=0;
        memset(ans,0,sizeof(ans));
        for(int i=1;i<=n;i++){
            memset(vis,0,sizeof(vis));
            if(dfs(i)) cnt++;
        }
        if(cnt==n) puts("Yes"); else puts("No");
    }
    return 0;
}

题外话

表示现在上课认真听的人真的很少啊
除了我们这帮人全在颓，都不知道在干嘛
只能说后排是vip专享区……