动态规划考试2(test20170401)

前言

这次考试虽然称作dp考试，可怎么觉得第一题好像没什么关系。。。

神圣罗马帝国皇帝(emperor.pas/c/cpp)

注：Time Limit:2s Memory Limit:256MB

【题目描述】

题目描述只是借了历史的名义而已，纯属虚构娱乐，我亵渎历史我有罪。
中世纪的德意志，战火纷飞。近来ZJ、JS、HN、SD等地区为了扩大自己的版图面积，在内部分别举行了ZJOI、JSOI、HNOI以及SDOI(或者叫SDTSC)选拔了许多优秀的战士准备发动一场波及到全国各处的战争——NOI！
在HN举行HNOI过程中，金天成，绰号光下巴（或者没胡子）红鼻子，二试怒翻盘的成绩强势成为HN大军的一员并被封为公爵，率领整个HN大军。HN这个地方早已不能让金天成满足了！他在NOI之前发动了一场名为PKUSC的战役。金天成凭借着强大的军事才能，卓越的领导能力，强势统一了德意志的大部分地区。为他成为神圣罗马帝国的皇帝打下了坚实的基础！
现在，金天成需要加强自己统治区域内的运输，他现在有N块地盘，每块地盘之间都有些道路，总共有$M$条道路，第i条道路有两个属性：运输量$Ci$、危险指数$Di$。金天成希望有一种方案，使得维护最少的道路让$N$块地盘两两之间存在通路。并且选取出来的道路集合$S$的评估指数$Eval = \frac{\sum_{i\in S}C_i}{\sum_{i\in S}D_i}$最大！由于他日理万机，所以这件事情他就让你来做了……

【输入格式】

第一行两个数字$N$，$M$。
接下来M行，每行三个数字$A,B,C,D(1\le A,B\le N)$。表示A、B两块地盘之间有一条运输量为C，危险指数为D的可用通路。

【输出格式】

输出一行，表示所有道路选取方案中最大的评估指数。保留6位小数。

【样例输入输出】

emperor.in	emperor.out
3 3	2.000000
1 2 2 1
2 3 2 1
3 1 1 1

【样例解释】

总共三种选法{边1,边2},{边2,边3},{边3,边1}，然后评估值分别为2，1.5，1.5。最大评估值为2。

【数据范围】

10%的数据：$1≤N≤5$，$1≤M≤10$。
30%的数据：$1≤N≤100$，$1≤M≤1000$。
100%的数据：$1≤N≤10^4$，$1≤M≤5*10^4$。$1≤C≤10^4$，$1≤D≤10^4$。

【题解】

当 Eval 值最大时也会存在评估指数小于 Eval 的，但不可能大于 Eval ，所以我们就二分答案，枚举 Eval ，在图中寻找最大生成树，如果边权之和大于 0 就可能，否则不可能。其中边权为$vi=Ci−Di∗Eval$
由 Eval 的式子得出。
所以说这题目和dp有什么关系，难道仅仅和01分数规划有关吗？

【代码】

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;

const int N = 10100;
const int M = 50100;
const double eps = 1e-7/2.0;
struct Edge{int u,v,c,d; double w;}a[M];
bool cmp(const Edge& a1,const Edge& a2){return a1.w<a2.w;}
int p[N]; int n,m;

inline int read() {
    int x=0; char ch=getchar();
    while(ch<'0' || ch>'9') ch=getchar();
    while(ch>='0' && ch<='9'){
        x=x*10+ch-'0';
        ch=getchar();
    }
    return x;
}

int Find(int u) {
    if(p[u]!=u) p[u]=Find(p[u]);
    return p[u];
}

bool Union(int u,int v) {
    u=Find(u); v=Find(v);
    if(u==v) return false;
    p[u]=v; return true;
}

bool Kruskal(double Eval) {
    double MST=0; int k=0;
    for(int i=1;i<=m;i++)
        a[i].w=(double)(Eval*a[i].d-a[i].c);
    for(int i=1;i<=n;i++) p[i]=i;
    sort(a+1,a+m+1,cmp);
    for(int i=1;i<=m && k<n-1;i++){
        if(Union(a[i].u,a[i].v)){
            MST+=a[i].w; k++;
        }
    }
    return MST>=-eps;
}

double Solve(double l,double r) {
    double mid;
    while(r-l>eps){
        mid=(l+r)/2;
        if(Kruskal(mid)) r=mid;
        else l=mid;
    }
    return l;
}

int main() {
    double maxx=0;
    freopen("emperor.in","r",stdin);
    freopen("emperor.out","w",stdout);
    n=read(); m=read();
    for(int i=1;i<=m;i++){
        a[i].u=read(); a[i].v=read();
        a[i].c=read(); a[i].d=read();
        maxx=max(maxx,(double)(a[i].c/a[i].d));
    }
    printf("%.6lf\n",Solve(0,maxx));
    return 0;
}

猛狮亨利(lion.pas/c/cpp)

注：Time Limit:1s Memory Limit:256MB

【题目描述】

题目描述只是借了历史的名义而已，纯属虚构娱乐，我亵渎历史我有罪。
自金天成在德意志各个公国树立威信后。他准备开始攻打波兰（好像中间少了什么……），后勤工作由萨克森和巴伐利亚两个公国负责（JS、ZJ、SD都酱油去了）。
而在金天成攻打波兰期间。德意志最强大的王公“饭桌猛狮亨利”，别名“二胖”，宣布他才是德意志的国王！发动了叛乱。
为了食物，金天成国王表示“一定要胖揍这个（哔~~~~~~~~~~~~~~~~~~~~~~）”于是，他留了部分人在战斗前线，带领一队人马前往镇压叛乱。
这时，他发现，他能带走的粮食有限，而士兵都是有饭量的。比如某位叫“巨胖”的军官可以使用碾压战术杀死一条线上的所有敌人（无视兵种），也可以当城墙使用，但是他一顿可以吃一吨。又比如某位叫“中二现”的军官，他饭量不多，但是杀伤力没“巨胖”军官强。这就让金天成犯难了。因为金天成一顿也可以吃一吨……所以，金天成希望能够在减去他的饭量后所剩下粮食数量有限的情况下，带领的士兵杀伤力最大，且不至于断粮而饿死沙场。
现在，他手下有N类可以带走的士兵或者军官，数量无限。每类士兵或军官都有个饭量Ai和杀伤力Bi，行军时每顿饭可以提供的粮食数量为M（不包括金天成的饭量）。现在请告诉金天成他能带走的杀伤力最强的部队杀伤力为多少。

【输入格式】

第1行两个数字N，M。
第2行到第n+1行，第i+1行两个数字Ai，Bi，表示第i个士兵（或军官）的饭量和杀伤力。

【输出格式】

输出一行，即最强部队杀伤力数值。

【样例输入输出】

lion.in	lion.out
4 10	16
2 3
3 5
4 6
5 7

【样例解释】

带2个2号士兵（或军官）和1个3号士兵（或军官）。饭量10，杀伤力16。

【数据范围】

30%的数据：$1≤N≤10$，$1≤M≤20$。
100%的数据：$1≤N≤10000$，$1≤M≤10000$，$1≤Ai≤10000$，$1≤Bi≤10000$。

【题解】

这个题目如果不懂再看一遍题，就会发现就是就是就是完全背包。。。。。。

【代码】

#include <cstdio>
#include <algorithm>
using namespace std;

const int size = 10000+10;
int n,m;
int w[size],c[size];
int f[size];

int main() {
    freopen("lion.in","r",stdin);
    freopen("lion.out","w",stdout);
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
    scanf("%d%d",&w[i],&c[i]);
    for(int i=1;i<=n;i++)
    for(int j=w[i];j<=m;j++)
        f[j]=max(f[j],f[j-w[i]]+c[i]);
    printf("%d\n",f[m]);
    return 0;
}

教皇与僭称教皇(pope.pas/c/cpp)

【题目描述】

题目描述只是借了历史的名义而已，纯属虚构娱乐，我亵渎历史我有罪。
稳固了后方之后，光下巴开始寻求加冕为神圣罗马帝国皇帝，1153年发动了第一次远征意大利。在进军途中，他接到了教皇尤金三世的求救信，原来是教士阿诺德领导人民夺取了城市政权，选举了罗马元老院成员和执政官，主张教会放弃领地过使徒的简朴生活。这真是天赐的机会，金天成加速进军镇压了这次起义，以拯救者的身份进入罗马城。
但新任教皇哈德良四世是个不懂变通的倔老头，他在为金天成加冕时，坚持要金天成按照惯例为教皇牵马、扶镫，激怒了金天成，结果使庆典会场变成了战场，近千人被杀。但不管过程如何狼狈，金天成已通过加冕获得了“神圣罗马帝国皇帝”的称号（在他之前的德皇，一直只称为“罗马帝国皇帝”，从他开始又冠上“神圣”二字）。
好吧，现在出题人要管管这不知道如何狼狈过程了，教堂的区域很少，只有N块区域。然后某些区域之间有通路。这N块区域中每块区域都有敌军，N个区域之间有M条通道。金天成大帝现在需要以闪电般的速度将这N个区域中敌军迅速歼灭。金天成大帝的军队比教皇军队要强很多，他们到达某个区域就能瞬间歼灭该区域的军队。但是走第i条通道花的时间为Ti。由于某种特殊的作者编造不出来的原因导致金天成大帝不能回到以前曾经过的区域，所以金天成大帝希望有一种方案，使得他能最快歼灭教皇的军队，如果不存在方案，则输出’-1’。一开始金天成大帝所在区域为1号区域。

【输入格式】

第一行两个数字N，M。
接下来M行，每行三个数字A,B,C，表示A到B有一条通路，耗时为C。

【输出格式】

输出一行，一个数字，表示金天成大帝干掉教皇军队的最少时间。

【样例输入输出】

pope.in	pope.out
5 4	18
1 2 3
2 3 4
3 4 5
4 5 6

【样例解释】

从1开始一条路走到黑就OK。

【数据范围】

30%的数据：$1≤N≤8$，$1≤M≤28$。
100%的数据：$1≤N≤16$，$1≤M≤120$。$1≤C≤10^8$

【题解】

刚开始打搜索，竟然能够过70分
注意到这里的 N 很小，可以采用状态压缩，设 $dis_{i,j}$ 表示从 1 走到 i 时且状态为 j 的最短时间，其中 j 是一个二进制串，剩下求 $dis_{i,j}$

【代码】

70分代码

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

typedef long long LL;
const int M = 150;
const int N = 20;
const LL INF = 9223372036854775807LL;
struct node {
    int u,v,w;
}e[M];

struct edge {
    int u,v,w,nxt;
}a[M<<1];

bool vis[N];
int head[N],n,m,cnt=0;
LL ans=INF;

inline bool cmp(const node &e1,const node &e2) {
    return e1.w>e2.w;
}

inline int read() {
    int x=0; char ch=getchar();
    while(ch<'0' || ch>'9') ch=getchar();
    while(ch>='0' && ch<='9'){
        x=x*10+ch-'0';
        ch=getchar();
    }
    return x;
}

inline void add_edge(int u,int v,int w) {
    a[++cnt]=(edge){u,v,w,head[u]};head[u]=cnt;
}

inline void dfs(int u,int x,LL now) {
    if(now>ans) return;
    if(x==n) {
    ans=now; return;
    }
    vis[u]=true;
    for(int i=head[u];i;i=a[i].nxt) {
    int v=a[i].v;
    if(vis[v]) continue;
    dfs(v,x+1,now+a[i].w);
    }
    vis[u]=false;
}

inline void init() {
    n=read();m=read();
    for(int i=1;i<=m;i++) {
    e[i].u=read();e[i].v=read();e[i].w=read();
    }
    sort(e+1,e+m+1,cmp);
    for(int i=1;i<=m;i++) {
    add_edge(e[i].u,e[i].v,e[i].w);
    add_edge(e[i].v,e[i].u,e[i].w);
    }
}

inline void output() {
    if(ans==INF) puts("-1");
    else printf("%lld\n",ans);
}

int main() {
    freopen("pope.in","r",stdin);
    freopen("pope.out","w",stdout);
    init();
    dfs(1,1,0);
    output();
    return 0;
}

AC代码

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cctype>
using namespace std;
typedef long long LL;
#define Min(a,b) ((a)<(b)?(a):(b))
#define now node[point].data

inline int read() {
    int in=0,f=1;
    char ch=getchar();
    for(;ch<'0'||ch>'9';ch=getchar())
    if(ch=='-')
        f=-1;
    for(;ch>='0'&&ch<='9';ch=getchar())
    in=in*10+ch-'0';
    return in*f;
}

const int maxn=20;
const LL INF=9*1e13;
struct Node{
    int data,next,weight;
}node[310];
int n,m,cnt;
int head[maxn];
LL ans;
LL f[(1<<16)+10][maxn];

inline void add_edge(int u,int v,int w) {
    node[cnt].data=v;node[cnt].next=head[u];node[cnt].weight=w;head[u]=cnt++;
    node[cnt].data=u;node[cnt].next=head[v];node[cnt].weight=w;head[v]=cnt++;
}

int main() {
    freopen("pope.in","r",stdin);
    freopen("pope.out","w",stdout);
    n=read();m=read();
    for(int i=1;i<=n;i++) head[i]=-1;
    for(int i=1,u,v,w;i<=m;i++) u=read(),v=read(),w=read(),add_edge(u,v,w);
    int Lim=(1<<n);
    memset(f,0x3f,sizeof(f));
    f[1][1]=0;
    for(int i=1;i<Lim;i++){
    for(int j=1;j<=n;j++)
        if((1<<(j-1))&i)
        for(int point=head[j];point!=-1;point=node[point].next)
            if((i&(1<<(now-1)))==0)
            f[i|(1<<(now-1))][now]=Min(f[i|(1<<(now-1))][now],f[i][j]+node[point].weight);
    }
    ans=INF;
    for(int i=1;i<=n;i++)
    ans=Min(ans,f[Lim-1][i]);
    if(ans==INF) printf("-1");
    else printf("%lld\n",ans);
    return 0;
}

总结

还是考炸了啊