2718: [Violet 4]毕业旅行/1143: [CTSC2008]祭祀river 最长反链=最小边覆盖

这个证明我真的不会啊……
Vfleaking证明传送门
Dilworth定理:最小链覆盖数 = 最长反链长度。
其对偶定理:最长链长度 = 最小反链覆盖数。

知道这个定理后,先floyd求出点对间的关系再跑dinic就好了。

#include<iostream>
#include<cstdio>
#include<cstring>
#define inf 1000000007
#define N 500
using namespace std;
int n,m,T,ans,cnt=1,a[N>>1][N>>1],dis[N],q[N],head[N];
int next[N*N],list[N*N],key[N*N];
inline int read()
{
    int a=0,f=1; char c=getchar();
    while (c<'0'||c>'9') {if (c=='-') f=-1; c=getchar();}
    while (c>='0'&&c<='9') {a=a*10+c-'0'; c=getchar();}
    return a*f;
}
inline void insert(int x,int y,int z)
{
    next[++cnt]=head[x];
    head[x]=cnt;
    list[cnt]=y;
    key[cnt]=z;
}
inline void floyd()
{
    for (int k=1;k<=n;k++)
        for (int i=1;i<=n;i++)
            for (int j=1;j<=n;j++)
                a[i][j]|=(a[i][k]&a[k][j]);
}
inline bool BFS()
{
    memset(dis,-1,sizeof(dis));
    dis[0]=1; q[1]=0;
    int t=0,w=1,x;
    while (t!=w)
    {
        x=q[++t];
        for (int i=head[x];i;i=next[i])
            if (key[i]&&dis[list[i]]==-1)
                dis[list[i]]=dis[x]+1,q[++w]=list[i];
    }
    return dis[T]!=-1;
}
int find(int x,int flow)
{
    if (x==T) return flow;
    int w,used=0;
    for (int i=head[x];i;i=next[i])
        if (key[i]&&dis[list[i]]==dis[x]+1)
        {
            w=find(list[i],min(key[i],flow-used));
            key[i]-=w; key[i^1]+=w; used+=w;
            if (used==flow) return flow;
        }
    if (!used) dis[x]=-1;
    return used;
}
inline void build()
{
    for (int i=1;i<=n;i++)
        for (int j=1;j<=n;j++)
            if (a[i][j]) insert(i,j+n,inf),insert(j+n,i,0);
    for (int i=1;i<=n;i++)
        insert(0,i,1),insert(i,0,0),insert(i+n,T,1),insert(T,i+n,0);
}
int main()
{
    n=read(); m=read(); T=2*n+1;
    for (int i=1;i<=m;i++)
    {
        int u=read(),v=read();
        a[u][v]=1;
    }
    floyd();
    build();
    while (BFS()) ans+=find(0,inf);
    cout << n-ans;
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
import matplotlib.pyplot as plt import np as np import numpy as np from scipy import signal from scipy import fftpack import matplotlib.font_manager as fm t = np.linspace(-1, 1, 200, endpoint=False) x = (np.cos(2,np.pi5t) + np.sin(2np.pi20t) * np.exp(-t**3/0.4)) X = fftpack.fft(x) fig, axs = plt.subplots(2, 2, figsize=(16, 8)) axs[0, 0].plot(t, x, color='pink') axs[0, 0].set_title('原信号', fontproperties=fm.FontProperties(fname='C:/Windows/Fonts/simsun.ttc'), color='plum') axs[0, 0].tick_params(axis='x', colors='red') axs[0, 0].tick_params(axis='y', colors='blue') axs[0, 1].plot(t, np.abs(X), color='brown') axs[0, 1].set_title('傅里叶变换', fontproperties=fm.FontProperties(fname='C:/Windows/Fonts/simsun.ttc'), color='violet') axs[0, 1].set_ylim([0, 25]) axs[0, 1].tick_params(axis='x', colors='red') axs[0, 1].tick_params(axis='y', colors='blue') b1, a1 = signal.butter(16, 0.2) y = signal.filtfilt(b1, a1, x) axs[1, 0].plot(t, y, color='grey') axs[1, 0].set_title('高通滤波', fontproperties=fm.FontProperties(fname='C:/Windows/Fonts/simsun.ttc'), color='indigo') axs[1, 0].tick_params(axis='x', colors='red') axs[1, 0].tick_params(axis='y', colors='blue') b2, a2 = signal.butter(4, 0.3) z = signal.filtfilt(b2, a2, x) axs[1, 1].plot(t, z, color='orange') axs[1, 1].set_title('低通滤波', fontproperties=fm.FontProperties(fname='C:/Windows/Fonts/simsun.ttc'), color='navy') axs[1, 1].tick_params(axis='x', colors='red') axs[1, 1].tick_params(axis='y', colors='blue') plt.tight_layout() plt.show()有错误
06-02

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值