• BZOJ1433: [ZJOI2009]假期的宿舍


    1433: [ZJOI2009]假期的宿舍

    Time Limit: 10 Sec  Memory Limit: 162 MB
    Submit: 1141  Solved: 491
    [Submit][Status]

    Description

    Input

    Output

    Sample Input

    1
    3
    1 1 0
    0 1 0
    0 1 1
    1 0 0
    1 0 0

    Sample Output

    ˆ ˆ

    HINT

    对于30% 的数据满足1 ≤ n ≤ 12。
    对于100% 的数据满足1 ≤ n ≤ 50,1 ≤ T ≤ 20。

    题解:
    左边是人,右边是床,dinic即可
    注意逻辑关系即可
    代码:
     1 #include<cstdio>
     2 #include<cstdlib>
     3 #include<cmath>
     4 #include<cstring>
     5 #include<algorithm>
     6 #include<iostream>
     7 #include<vector>
     8 #include<map>
     9 #include<set>
    10 #include<queue>
    11 #define inf 1000000000
    12 #define maxn 100+100
    13 #define maxm 2500+100
    14 #define eps 1e-10
    15 #define ll long long
    16 using namespace std;
    17 inline int read()
    18 {
    19     int x=0,f=1;char ch=getchar();
    20     while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    21     while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}
    22     return x*f;
    23 }
    24 struct edge{int go,next,v;}e[2*maxm];
    25 int n,m,s,t,tot=1,maxflow,head[maxn],q[maxn],h[maxn],cur[maxn];
    26 bool a[maxn],b[maxn];
    27 void ins(int x,int y,int z){e[++tot].go=y;e[tot].v=z;e[tot].next=head[x];head[x]=tot;}
    28 void insert(int x,int y,int z){ins(x,y,z);ins(y,x,0);}
    29 bool bfs()
    30 {
    31     for(int i=s;i<=t;i++)h[i]=-1;
    32     int l=0,r=1;q[1]=s;h[s]=0;
    33     while(l<r) 
    34     {
    35         int x=q[++l];
    36         for(int i=head[x];i;i=e[i].next)
    37          if(e[i].v&&h[e[i].go]==-1)
    38          {
    39              h[e[i].go]=h[x]+1;q[++r]=e[i].go;
    40          }
    41     }
    42     return h[t]!=-1;
    43 }
    44 int dfs(int x,int f)
    45 {
    46     if(x==t) return f;
    47     int tmp,used=0,i;
    48     for(int i=head[x];i;i=e[i].next)
    49      if(e[i].v&&h[e[i].go]==h[x]+1) 
    50     {
    51         tmp=dfs(e[i].go,min(e[i].v,f-used));
    52         e[i].v-=tmp;if(e[i].v)cur[x]=i; 
    53         e[i^1].v+=tmp;used+=tmp;
    54         if(used==f)return f;         
    55     }
    56     if(!used) h[x]=-1;
    57     return used; 
    58 }
    59 void dinic()
    60 {
    61     maxflow=0;
    62     while(bfs())
    63     {
    64         for (int i=s;i<=t;i++)cur[i]=head[i];maxflow+=dfs(s,inf); 
    65     } 
    66 }
    67 int main()
    68 {
    69     freopen("input.txt","r",stdin);
    70     freopen("output.txt","w",stdout);
    71     int cs=read();
    72     while(cs--)
    73     {    
    74      memset(a,0,sizeof(a));
    75      memset(b,0,sizeof(b));
    76      memset(head,0,sizeof(head));tot=1;
    77      n=read();int x;
    78      for(int i=1;i<=n;i++)a[i]=read();
    79      for(int i=1;i<=n;i++){x=read();if(x&&a[i])b[i]=1;}
    80      //for(int i=1;i<=n;i++)cout<<i<<' '<<a[i]<<' '<<b[i]<<endl;
    81      s=0;t=2*n+1;
    82      for(int i=1;i<=n;i++)
    83      {
    84          if(!b[i]&&a[i])insert(i,i+n,1);
    85         for(int j=1;j<=n;j++){x=read();if(!b[i]&&x&&a[j])insert(i,j+n,1);}
    86      }
    87      int ans=0;
    88      for(int i=1;i<=n;i++)if(!b[i])insert(s,i,1),ans++;
    89      for(int i=1;i<=n;i++)if(a[i])insert(i+n,t,1);
    90      dinic();
    91      //cout<<ans<<' '<<maxflow<<endl;
    92      if(maxflow==ans)printf("%c%c%c
    ",94,95,94);else printf("%c%c%c
    ",84,95,84);
    93     }
    94     return 0;
    95 }
    View Code
  • 相关阅读:
    『实践』Yalmip建模+Cplex类求解(文末附程序、文章和算例)
    『实践』Matlab实现Flyod求最短距离及存储最优路径
    『转载』Matlab中fmincon函数获取乘子
    『实践』Yalmip获取对偶函数乘子
    『转载』hadoop 1.X到2.X的变化
    『转载』hadoop2.x常用端口、定义方法及默认端口
    『实践』VirtualBox 5.1.18+Centos 6.8+hadoop 2.7.3搭建hadoop完全分布式集群及基于HDFS的网盘实现
    C#设计模式学习笔记:(6)适配器模式
    ASP.NET 开源导入导出库Magicodes.IE 导出Pdf教程
    Rx基础
  • 原文地址:https://www.cnblogs.com/zyfzyf/p/3930998.html
Copyright © 2020-2023  润新知