• <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

            最新隨筆

            搜索

            積分與排名

            最新評論

            閱讀排行榜

            无码精品久久一区二区三区 | 91麻精品国产91久久久久| 99久久人妻无码精品系列| 欧美精品一本久久男人的天堂| 日韩精品国产自在久久现线拍 | 久久久精品无码专区不卡| 色婷婷久久久SWAG精品| 囯产极品美女高潮无套久久久 | 精品国产VA久久久久久久冰| 久久九九青青国产精品| 久久久亚洲欧洲日产国码是AV| 亚洲va久久久噜噜噜久久| 久久亚洲高清观看| 亚洲精品无码专区久久久| 日本精品久久久久中文字幕| 久久久久久综合网天天| 久久93精品国产91久久综合| 狠狠色婷婷久久一区二区 | 久久国产高清一区二区三区| 久久综合亚洲欧美成人| 久久精品成人免费观看97| 久久夜色精品国产欧美乱| 人妻系列无码专区久久五月天| 国产精品久久久久影院嫩草 | 久久国产免费观看精品| 国内精品久久久久久久久电影网 | 久久久久久综合一区中文字幕| 久久99国产精品久久99小说| 国产精品无码久久久久| 伊人久久免费视频| 亚洲综合婷婷久久| 9191精品国产免费久久| 国产一级持黄大片99久久| 国内精品伊人久久久久| 成人久久综合网| www.久久热| 国产日韩久久久精品影院首页| 国内精品久久久久国产盗摄| 精品无码久久久久久久久久| 久久久久久A亚洲欧洲AV冫| 亚洲国产成人精品无码久久久久久综合 |