• <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 閱讀(347) 評(píng)論(0)  編輯 收藏 引用 所屬分類: graph

            <2010年10月>
            262728293012
            3456789
            10111213141516
            17181920212223
            24252627282930
            31123456

            導(dǎo)航

            統(tǒng)計(jì)

            公告

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

            留言簿(1)

            隨筆分類(227)

            文章分類(2)

            OJ

            最新隨筆

            搜索

            積分與排名

            最新評(píng)論

            閱讀排行榜

            伊人久久大香线蕉av不卡| 久久精品国产亚洲AV香蕉| 久久精品午夜一区二区福利| 亚洲国产精品无码久久久久久曰| 99热精品久久只有精品| 人妻少妇久久中文字幕一区二区| 噜噜噜色噜噜噜久久| 久久亚洲国产最新网站| 欧美激情精品久久久久久久| 久久中文字幕视频、最近更新| 国产精品热久久毛片| 国产精品日韩深夜福利久久| 久久久久综合网久久| 国产精品成人精品久久久| 日本一区精品久久久久影院| 99久久精品免费看国产免费| 精品熟女少妇aⅴ免费久久| 国产亚洲精午夜久久久久久| 久久九色综合九色99伊人| 久久久久久国产a免费观看不卡| 日韩va亚洲va欧美va久久| 亚洲国产成人久久综合碰| 国内精品久久久久影院亚洲| 久久婷婷五月综合97色直播 | 亚洲日本va中文字幕久久| 亚洲精品午夜国产VA久久成人| 亚洲AV日韩AV永久无码久久| 久久久噜噜噜久久中文福利| 狠狠色噜噜狠狠狠狠狠色综合久久 | 少妇人妻综合久久中文字幕| 亚洲AV日韩精品久久久久| 99久久免费国产特黄| 94久久国产乱子伦精品免费| 日批日出水久久亚洲精品tv| 人妻少妇久久中文字幕一区二区 | 欧美精品一本久久男人的天堂| 久久精品国产精品亚洲艾草网美妙| 性做久久久久久久久| 久久精品国产网红主播| 久久er热视频在这里精品| 国产日韩欧美久久|