神圣罗马帝国皇帝(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 条道路有两个属性:运输量 C i 、危险指数 D i 。金天成希望有一种方案,使得维护最少的道路让 N 块地盘两两之间存在通路。并且选取出来的道路集合 S 的评估指数 Eval=
∑i∈SCi∑i∈SDi
最大!由于他日理万机,所以这件事情他就让你来做了……
输入
第一行两个数字 N,M。
接下来 M 行,每行三个数字 A,B,C,D(1≤A,B≤N)。表示 A、B 两块地盘之间有一条运输
量为 C,危险指数为 D 的可用通路。
输出
输出一行,表示所有道路选取方案中最大的评估指数。保留 6 位小数。
输入样例
3 3
1 2 2 1
2 3 2 1
3 1 1 1
输出样例
2.000000
样例解释
总共三种选法{边 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∗104
。1≤C≤
104
,1≤D≤
104
。
分析
分数规划+最大生成树
代码
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int sm = 10000*5+10;
const double EPS=1e-7;
int n,m,f[sm],u,v;
long double l,r,mid,res;
struct edge{
int from,to,c,d;
}e[sm];
struct tmp{
int from,to;
long double val;
}tt[sm];
int find(int x)
{
if(x!=f[x])f[x]=find(f[x]);
return f[x];
}
bool cmp(tmp x,tmp y){return x.val>y.val;}
bool judge(long double Q)
{
res=0;
for(int i=1;i<=m;++i)
{
tt[i].from=e[i].from;
tt[i].to=e[i].to;
tt[i].val=e[i].c-Q*e[i].d;
}
sort(tt+1,tt+m+1,cmp);
for(int i=1;i<=n;++i)f[i]=i;
for(int i=1;i<=m;++i)
{
u=find(tt[i].from);
v=find(tt[i].to);
if(v!=u)
{
f[u]=v;
res+=tt[i].val;
}
}
if(res>=0)return 1;
return 0;
}
int main()
{
freopen("emperor.in","r",stdin);
freopen("emperor.out","w",stdout);
scanf("%d%d",&n,&m);
for(int i=1;i<=m;++i)
scanf("%d%d%d%d",&e[i].from,&e[i].to,&e[i].c,&e[i].d);
l=0;r=10000;
while(l+EPS<r)
{
mid=(l+r)/2.0;
if(judge(mid))l=mid;
else r=mid;
}
printf("%.6Lf\n",l);
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 行两个数字 A i ,B i ,表示第 i 个士兵(或军官)的饭量和杀伤力。
输出
输出一行,即最强部队杀伤力数值。
输入样例
4 10
2 3
3 5
4 6
5 7
输出样例
16
样例解释
带 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<cstring>
#include<algorithm>
using namespace std;
const int sm = 10000+20;
int n,m,w[sm],c[sm];
long long f[sm];
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-w[i]]+c[i],f[j]);
printf("%lld\n",f[m]);
return 0;
}
教皇与僭称教皇(pope.pas/c/cpp Time Limit:0.3s Memory Limit:256MB)
问题描述
题目描述只是借了历史的名义而已,纯属虚构娱乐,我亵渎历史我有罪。
稳固了后方之后,光下巴开始寻求加冕为神圣罗马帝国皇帝,1153 年发动了第一次远征意大利。在进军途中,他接到了教皇尤金三世的求救信,原来是教士阿诺德领导人民夺取了城市政权,选举了罗马元老院成员和执政官,主张教会放弃领地过使徒的简朴生活。这真是天赐的机会,金天成加速进军镇压了这次起义,以拯救者的身份进入罗马城。
但新任教皇哈德良四世是个不懂变通的倔老头,他在为金天成加冕时,坚持要金天成按照惯例为教皇牵马、扶镫,激怒了金天成,结果使庆典会场变成了战场,近千人被杀。但不管过程如何狼狈,金天成已通过加冕获得了“神圣罗马帝国皇帝”的称号(在他之前的德皇,一直只称为“罗马帝国皇帝”,从他开始又冠上“神圣”二字)。
好吧,现在出题人要管管这不知道如何狼狈过程了,教堂的区域很少,只有 N 块区域。
然后某些区域之间有通路。这 N 块区域中每块区域都有敌军,N 个区域之间有 M 条通道。金天成大帝现在需要以闪电般的速度将这 N 个区域中敌军迅速歼灭。金天成大帝的军队比教皇军队要强很多,他们到达某个区域就能瞬间歼灭该区域的军队。但是走第 i 条通道花的时间为 T i 。由于某种特殊的作者编造不出来的原因导致金天成大帝不能回到以前曾经过的区域,所以金天成大帝希望有一种方案,使得他能最快歼灭教皇的军队,如果不存在方案,则输出’-1’。一开始金天成大帝所在区域为 1 号区域。
输入
第一行两个数字 N,M。
接下来 M 行,每行三个数字 A,B,C,表示 A 到 B 有一条通路,耗时为 C。
输出
输出一行,一个数字,表示金天成大帝干掉教皇军队的最少时间。
输入样例
5 4
1 2 3
2 3 4
3 4 5
4 5 6
输出样例
18
样例解释
从 1 开始一条路走到黑就 OK
数据范围
30%的数据:1≤N≤8,1≤M≤28。
100%的数据:1≤N≤16,1≤M≤120。1≤C≤
108
分析
n非常小,状压DP,dis[i][j]表示走到i点且状态为j的最短时间,j用16位的二进制数表示第i个点是否被消灭,结合spfa即可
代码
#include<queue>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
struct node{int x,y;node(int X=0,int Y=0){x=X;y=Y;}};
const int sz=120+10;
int ans,r,s,n,m,a,b,c,h[sz],dis[20][1<<17],vis[20][1<<17];
struct edge{
int to,nxt,len;
}e[sz];
void build(int from,int to,int len)
{
e[++s].to=to;
e[s].nxt=h[from];
e[s].len=len;
h[from]=s;
}
void spfa()
{
queue<node> q;
node now;
q.push(node(0,1));
memset(dis,0x3f,sizeof(dis));
dis[0][1]=0;
while(!q.empty())
{
now=q.front();q.pop();
vis[now.x][now.y]=0;
for(int i=h[now.x];i;i=e[i].nxt)
{
if(!(now.y&(1<<e[i].to))&&e[i].len+dis[now.x][now.y]<dis[e[i].to][now.y|1<<e[i].to])
{
dis[e[i].to][now.y|1<<e[i].to]=e[i].len+dis[now.x][now.y];
if(!vis[e[i].to][now.y|1<<e[i].to])
{
vis[e[i].to][now.y|1<<e[i].to]=1;
q.push(node(e[i].to,now.y|1<<e[i].to));
}
}
}
}
}
int main()
{
freopen("pope.in","r",stdin);
freopen("pope.out","w",stdout);
scanf("%d%d",&n,&m);
for(int i=1;i<=m;++i)
{
scanf("%d%d%d",&a,&b,&c);
a--;b--;
build(a,b,c);build(b,a,c);
}
spfa();
r=(1<<n)-1;ans=0x3f3f3f3f;
for(int i=1;i<=n;++i)
ans=min(ans,dis[i][r]);
if(ans!=0x3f3f3f3f)printf("%d\n",ans);
else printf("-1\n");
return 0;
}