• <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 最小費(fèi)用流 有限制的區(qū)間覆蓋求最大權(quán)值和,經(jīng)典題目

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

            解法:
            對(duì)s1,e1,s2,e2..sn,en進(jìn)行排序并去重,得到有序數(shù)列a1,a2,a3..ak,然后建立網(wǎng)絡(luò),在ai與ai+1之間連一條容量為k,費(fèi)用為0的邊,對(duì)于原有區(qū)間(s,e),在s',e'(s',e'指離散化后的值)之間連一條容量為1,費(fèi)用為-c的邊,然后建立源點(diǎn),向a1連一條容量為k,費(fèi)用為0的邊;建立匯點(diǎn),從ak向其連容量為k,費(fèi)用為0的邊,最后求最小費(fèi)用流即可
            代碼:
              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 閱讀(359) 評(píng)論(0)  編輯 收藏 引用 所屬分類: graph

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

            導(dǎo)航

            統(tǒng)計(jì)

            公告

            統(tǒng)計(jì)系統(tǒng)

            留言簿(1)

            隨筆分類(227)

            文章分類(2)

            OJ

            最新隨筆

            搜索

            積分與排名

            最新評(píng)論

            閱讀排行榜

            久久久中文字幕| www久久久天天com| 久久久久久久综合狠狠综合| 日本亚洲色大成网站WWW久久| 亚洲欧美另类日本久久国产真实乱对白 | 91精品日韩人妻无码久久不卡| 国产精品欧美久久久久天天影视| 欧美午夜精品久久久久久浪潮| 久久精品国产久精国产一老狼| 婷婷综合久久中文字幕| 伊人久久大香线蕉av不变影院 | 国产69精品久久久久9999APGF| 久久精品人人做人人妻人人玩| 狠狠人妻久久久久久综合蜜桃| 亚洲色大成网站WWW久久九九| 久久久精品久久久久久 | 国产精品久久久久久福利漫画| 精品99久久aaa一级毛片| 久久精品黄AA片一区二区三区| 久久久久亚洲av成人无码电影| 99久久超碰中文字幕伊人| 久久狠狠爱亚洲综合影院| 亚洲国产精品综合久久一线| 久久这里只有精品久久| 日本精品久久久久中文字幕8 | 色综合久久久久网| 久久国产精品久久国产精品| www久久久天天com| 潮喷大喷水系列无码久久精品| 亚洲国产精品无码久久久蜜芽 | www性久久久com| 精品国产VA久久久久久久冰| 无码精品久久久久久人妻中字| 伊人色综合久久天天人守人婷| 久久久久久极精品久久久| 精品欧美一区二区三区久久久| 久久精品嫩草影院| 亚洲国产精品一区二区久久| 国内精品久久久久久久影视麻豆| 亚洲国产精品一区二区久久| 麻豆久久|