• <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>

            pku 3680 Intervals 最小費用流 有限制的區間覆蓋求最大權值和,經典題目

            簡明題意:
            給出一些開區間(s,e)以及權值c,要求選出一些區間,使得所有點的覆蓋次數小于k并且權值和最大。

            解法:
            對s1,e1,s2,e2..sn,en進行排序并去重,得到有序數列a1,a2,a3..ak,然后建立網絡,在ai與ai+1之間連一條容量為k,費用為0的邊,對于原有區間(s,e),在s',e'(s',e'指離散化后的值)之間連一條容量為1,費用為-c的邊,然后建立源點,向a1連一條容量為k,費用為0的邊;建立匯點,從ak向其連容量為k,費用為0的邊,最后求最小費用流即可
            代碼:
              1# include <cstdio>
              2# include <vector>
              3# include <cstring>
              4# include <algorithm>
              5using namespace std;
              6# define N 500                     //N+2
              7# define E 10000                    //2*E
              8#define typef int                     // type of flow
              9#define typec int                     // type of cost
             10const int inff = 0xfffffff;       // max of flow
             11const int infc = 0xfffffff;       // max of cost
             12struct network
             13{
             14  int nv, ne, pnt[E], nxt[E];
             15  int vis[N], que[N], head[N], pv[N], pe[N];
             16  typef flow, cap[E]; typec cost, dis[E], d[N];
             17  void addedge(int u, int v, typef c, typec w) 
             18  {
             19      pnt[ne] =  v; cap[ne] = c;
             20      dis[ne] = +w; nxt[ne] = head[u]; head[u] = (ne++);
             21      pnt[ne] =  u; cap[ne] = 0;
             22      dis[ne] = -w; nxt[ne] = head[v]; head[v] = (ne++);
             23  }

             24  int mincost(int src, int sink) 
             25  {
             26       int i, k, f, r; typef mxf;
             27       for (flow = 0, cost = 0; ; ) 
             28       {
             29           for(i=0;i<nv;i++) pv[i]=-1,d[i] = infc;
             30           memset(vis, 0sizeof(vis));
             31           d[src] = 0; pv[src] = src; vis[src] = 1;
             32           for (f = 0, r = 1, que[0= src; r != f; ) 
             33           {
             34                i = que[f++]; vis[i] = 0;
             35                if (N == f) f = 0;
             36                for (k = head[i]; k != -1; k = nxt[k])
             37                    if(cap[k] && dis[k]+d[i] < d[pnt[k]])
             38                    {
             39                        d[pnt[k]] = dis[k] + d[i];
             40                        if (0 == vis[pnt[k]]) 
             41                        {
             42                            vis[pnt[k]] = 1;
             43                            que[r++= pnt[k];
             44                            if (N == r) r = 0;
             45                        }

             46                        pv[pnt[k]]=i; pe[pnt[k]]=k;
             47                    }
              
             48            }

             49            if (-1 == pv[sink]) break;
             50            for (k = sink, mxf = inff; k != src; k = pv[k])
             51                if (cap[pe[k]] < mxf) mxf = cap[pe[k]];
             52            flow += mxf; cost += d[sink] * mxf;
             53            for (k = sink; k != src; k = pv[k]) 
             54            {
             55                cap[pe[k]] -= mxf; cap[pe[k] ^ 1+= mxf;
             56            }

             57       }

             58       return cost;
             59    }

             60  void build(int v) 
             61  {
             62      nv = v; ne = 0;
             63      for(int i=0;i<v;i++)
             64          head[i]=-1;
             65  }

             66  void add(int x,int y,int f,int w)//x->y up=f,cost=w
             67  {
             68      addedge(x, y, f, w);
             69  }

             70}
             net;
             71int data[400][3];
             72vector<int> map;
             73int main()
             74{
             75    int test;
             76    scanf("%d",&test);
             77    while(test--)
             78    {
             79       map.clear();
             80       int n,k;
             81       scanf("%d%d",&n,&k);
             82       for(int i=0;i<n;i++)
             83       {
             84         scanf("%d%d%d",&data[i][0],&data[i][1],&data[i][2]);
             85         map.push_back(data[i][0]);
             86         map.push_back(data[i][1]);
             87       }

             88       sort(map.begin(),map.end());
             89       vector<int>::iterator end=unique(map.begin(),map.end());
             90       while(map.end()!=end) map.pop_back();
             91       net.build(map.size()+2);
             92       for(int i=0;i<map.size()+1;i++)
             93          net.add(i,i+1,k,0);
             94       for(int i=0;i<n;i++)
             95           net.add(lower_bound(map.begin(),map.end(),data[i][0])-map.begin()+1,lower_bound(map.begin(),map.end(),data[i][1])-map.begin()+1,1,-data[i][2]);
             96       printf("%d\n",-net.mincost(0,map.size()+1));      
             97    }

             98    return 0;
             99}

            100

            posted on 2010-11-18 00:32 yzhw 閱讀(351) 評論(0)  編輯 收藏 引用 所屬分類: graph

            <2010年11月>
            31123456
            78910111213
            14151617181920
            21222324252627
            2829301234
            567891011

            導航

            統計

            公告

            統計系統

            留言簿(1)

            隨筆分類(227)

            文章分類(2)

            OJ

            最新隨筆

            搜索

            積分與排名

            最新評論

            閱讀排行榜

            久久99中文字幕久久| 欧洲人妻丰满av无码久久不卡| 77777亚洲午夜久久多喷| 国产麻豆精品久久一二三| 久久亚洲精品无码AV红樱桃| 狠狠色丁香婷婷久久综合不卡| 国产一区二区精品久久凹凸| 伊人久久一区二区三区无码| 97r久久精品国产99国产精| 久久久久久久综合日本| 久久香蕉超碰97国产精品| 伊人久久免费视频| 无遮挡粉嫩小泬久久久久久久| 精品国产乱码久久久久久浪潮 | 久久久久亚洲精品日久生情| 久久国产成人精品麻豆| 中文成人无码精品久久久不卡| 日产精品99久久久久久| 污污内射久久一区二区欧美日韩| 99久久免费国产特黄| 婷婷久久五月天| 无码精品久久一区二区三区| 精品久久久久久久无码| 一级做a爰片久久毛片免费陪| 66精品综合久久久久久久| 久久国产精品一国产精品金尊| 伊人色综合九久久天天蜜桃| 久久精品国产一区二区三区| 色综合久久精品中文字幕首页| 久久国产精品99国产精| 久久人人妻人人爽人人爽| 欧美成人免费观看久久| 亚洲美日韩Av中文字幕无码久久久妻妇 | 2021少妇久久久久久久久久| 色综合久久久久综合体桃花网 | 中文字幕久久精品| 麻豆久久| 久久久久人妻一区二区三区 | 无码精品久久一区二区三区| 香蕉99久久国产综合精品宅男自 | 色天使久久综合网天天|