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

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

            解法:
            對s1,e1,s2,e2..sn,en進行排序并去重,得到有序數(shù)列a1,a2,a3..ak,然后建立網(wǎng)絡(luò),在ai與ai+1之間連一條容量為k,費用為0的邊,對于原有區(qū)間(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 閱讀(347) 評論(0)  編輯 收藏 引用 所屬分類: graph

            <2011年3月>
            272812345
            6789101112
            13141516171819
            20212223242526
            272829303112
            3456789

            導(dǎo)航

            統(tǒng)計

            公告

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

            留言簿(1)

            隨筆分類(227)

            文章分類(2)

            OJ

            最新隨筆

            搜索

            積分與排名

            最新評論

            閱讀排行榜

            香蕉久久夜色精品国产尤物| 精品久久久久久国产牛牛app| 中文字幕久久波多野结衣av| 熟妇人妻久久中文字幕| 色综合合久久天天综合绕视看| 久久www免费人成看国产片| 久久受www免费人成_看片中文| 国产精品视频久久| 久久久久久伊人高潮影院| 欧美激情精品久久久久| 日本五月天婷久久网站| 狠狠色综合网站久久久久久久| 国产精品久久久香蕉| 99久久婷婷国产综合精品草原| 人妻无码久久精品| 久久这里只有精品久久| 亚洲中文字幕无码久久2017| 久久人人爽人人爽AV片| 久久精品国产免费一区| 久久99精品久久久久久久久久 | 国内精品人妻无码久久久影院 | 久久久久久久久久久久中文字幕 | 亚洲中文字幕久久精品无码APP| 精品免费久久久久国产一区| 青青草国产精品久久| 97久久国产亚洲精品超碰热| 无码人妻久久一区二区三区免费| 久久久久国色AV免费观看| 中文字幕亚洲综合久久2| 99久久er这里只有精品18| 久久精品蜜芽亚洲国产AV| 天天爽天天狠久久久综合麻豆| 青青草原综合久久大伊人| 无码人妻久久一区二区三区蜜桃 | 日韩久久久久中文字幕人妻| 久久久久亚洲AV无码去区首| 久久996热精品xxxx| 欧美日韩成人精品久久久免费看 | 囯产极品美女高潮无套久久久| 久久亚洲国产成人影院| 国产69精品久久久久9999APGF|