zoukankan      html  css  js  c++  java
  • 传染病防治

    传染病防治(爆搜+剪枝)

    研究表明,这种传染病的传播具有两种很特殊的性质;第一是它的传播途径是树型的,一个人X只可能被某个特定的人Y感染,只要Y不得病,或者是XY之间的传播途径被切断,则X就不会得病。第二是,这种疾病的传播有周期性,在一个疾病传播周期之内,传染病将只会感染一代患者,而不会再传播给下一代。这些性质大大减轻了蓬莱国疾病防控的压力,并且他们已经得到了国内部分易感人群的潜在传播途径图(一棵树)。但是,麻烦还没有结束。由于蓬莱国疾控中心人手不够,同时也缺乏强大的技术,以致他们在一个疾病传播周期内,只能设法切断一条传播途径,而没有被控制的传播途径就会引起更多的易感人群被感染(也就是与当前已经被感染的人有传播途径相连,且连接途径没有被切断的人群)。当不可能有健康人被感染时,疾病就中止传播。所以,蓬莱国疾控中心要制定出一个切断传播途径的顺序,以使尽量少的人被感染。你的程序要针对给定的树,找出合适的切断顺序。输入格式的第一行是两个整数n(1≤n≤300)和p。接下来p行,每一行有两个整数i和j,表示节点i和j间有边相连(意即,第i人和第j人之间有传播途径相连)。其中节点1是已经被感染的患者。

    ​ 这道题就是个爆搜+剪枝的题目,每次搜每一行有几个点,如果人数已经超过ans就剪掉。做题的时候我还在疑惑为什么不会爆栈,事实上可以证明搜索层数不大于logn。感觉写的有点乱。。

    #include <cstdio>
    #include <cstring>
    using namespace std;
    
    const int maxn=305, INF=1e9;
    
    class Graph{
    public:
        struct Edge{
            int to, next; Graph *belong;
            void set(int x, int y, Graph *g){
                to=x; next=y; belong=g; }
            inline int operator*(){ return to; }
            Edge& operator ++(){ return *this=belong->edge[next]; }
        };
        void addedge(int x, int y){
            edge[++cntedge].set(y, fir[x], this);
            fir[x]=cntedge;
        }
        Edge& getlink(int x){ return edge[fir[x]]; }
    private:
        int cntedge, fir[maxn];
        Edge edge[maxn*2];
    };
    
    Graph g;
    //now数组储存当前层被传染的结点
    int n, m, fa[maxn], cntnow;
    int now[maxn], ans;
    
    void dfs(int depth, int cnt){
        if (cnt>ans) return;
        Graph::Edge e;
        int tmp[maxn], t=0;
        //先把下一层的所有结点搞出来
        for (int i=1; i<=cntnow; ++i){
            e=g.getlink(now[i]);
            for (; *e; ++e){
                if (*e==fa[now[i]]) continue;
                fa[*e]=now[i];
                tmp[++t]=*e;
            }
        }
        if (t==0){
            if (cnt<ans) ans=cnt;
            return;
        }
        //屏蔽下一层的某一个结点
        for (int i=1; i<=t; ++i){
            cntnow=0;
            for (int j=1; j<=t; ++j)
                if (j!=i) now[++cntnow]=tmp[j];
            dfs(depth+1, cnt+cntnow);
        }
    }
    
    int main(){
        scanf("%d%d", &n, &m);
        int x, y;
        for (int i=0; i<m; ++i){
            scanf("%d%d", &x, &y);
            g.addedge(x, y); g.addedge(y, x);
        }
        now[++cntnow]=1;
        ans=INF;
        dfs(1, 1);
        printf("%d
    ", ans);
        return 0;
    }
    
  • 相关阅读:
    宏定义中的#
    HDU1506 Largest Rectangle in a Histogram 动态规划
    HDU1864 最大报销额 DP
    POJ2771 Guardian of Decency 最大独立子集
    POJ1698 Alice's Chance 最大流
    HDU1003 Max Sum 动态规划
    eval格式化事件类型的字符串
    C#虚方法virtual详解
    c# 利用反射获得某个类或者对象的所有属性
    windows服务的通常写法
  • 原文地址:https://www.cnblogs.com/MyNameIsPc/p/7714421.html
Copyright © 2011-2022 走看看