NOIP模拟题 2016.9.24 [贪心] [有依赖的背包问题] [图论] [spfa或tarjan缩点+DAGdp]

9 篇文章 0 订阅
9 篇文章 0 订阅

1.排座椅
(seat.pas/c/cpp)
【问题描述】
上课的时候总有一些同学和前后左右的人交头接耳,这是令小学班主任十分头疼的一件事情。不过,班主任小雪发现了一些有趣的现象,当同学们的座次确定下来之后,只有有限的D对同学上课时会交头接耳。同学们在教室中坐成了M行N列,坐在第i行第j列
的同学的位置是(i,j),为了方便同学们进出,在教室中设置了K条横向的通道,L条纵向的通道。于是,聪明的小雪想到了一个办法,或许可以减少上课时学生交头接耳的问题:她打算重新摆放桌椅,改变同学们桌椅间通道的位置,因为如果一条通道隔开了两个会交头接耳的同学,那么他们就不会交头接耳了。
请你帮忙给小雪编写一个程序,给出最好的通道划分方案。在该方案下,上课时交头接耳的学生对数最少。
【输入】
输入文件seat.in的第一行,有5各用空格隔开的整数,分别是M,N,K,L,D(2<=N,M<=1000,0<=K< M,0<=L< N,D<=2000)。
接下来D行,每行有4个用空格隔开的整数,第i行的4个整数Xi,Yi,Pi,Qi,表示坐在位置(Xi,Yi)与(Pi,Qi)的两个同学会交头接耳(输入保证他们前后相邻或者左右相邻)。
输入数据保证最优方案的唯一性。
【输出】
输出文件seat.out共两行。
第一行包含K个整数,a1a2„„aK,表示第a1行和a1+1行之间、第a2行和第a2+1行之间、„、第aK行和第aK+1行之间要开辟通道,其中ai< ai+1,每两个整数之间用空格隔开(行尾没有空格)。
第二行包含L个整数,b1b2„„bk,表示第b1列和b1+1列之间、第b2列和第b2+1列之间、„、第bL列和第bL+1列之间要开辟通道,其中bi< bi+1,每两个整数之间用空格隔开(行尾没有空格)。
【输入输出样例】
seat.in
seat.out
4 5 1 2 3
4 2 4 3
2 3 3 3
2 5 2 4
2
2 4
【输入输出样例解释】
* *


+ +
1 2 3 4 5
上图中用符号*、※、+ 标出了3 对会交头接耳的学生的位置,图中3 条粗线的位置表示通
道,图示的通道划分方案是唯一的最佳方案。


这道题还好,直接每次贪心找最大值即可,用排序的方法。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<vector>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<string>
#include<iomanip>
#include<ctime>
#include<climits>
#include<cctype>
#include<algorithm>
#ifdef WIN32
#define AUTO "%I64d"
#else
#define AUTO "%lld"
#endif
using namespace std;
#define smax(x,tmp) x=max((x),(tmp))
#define smin(x,tmp) x=min((x),(tmp))
#define maxx(x1,x2,x3) max(max(x1,x2),x3)
#define minn(x1,x2,x3) min(min(x1,x2),x3)
const int INF=0x3f3f3f3f;
const int maxn = 1005;
int row[maxn];
int column[maxn];
int M,N,K,L,D;
void init()
{
    scanf("%d%d%d%d%d",&M,&N,&K,&L,&D);
    for(int i=1;i<=D;i++)
    {
        int x,y,p,q;
        scanf("%d%d%d%d",&x,&y,&p,&q);
        if(x==p)
        {
            if(y>q) swap(y,q);
            column[y]++;
        }
        else
        {
            if(x>p) swap(x,p);
            row[x]++;
        }
    }
}
struct Node
{
    int val;
    int id;
    bool operator < (const Node t) const { return val < t.val; }
};
priority_queue <Node> c,r;
int ans[maxn],top;
void work()
{
    for(int i=1;i<M;i++) r.push((Node){row[i],i});
    if(K)
    {
        top=0;
        for(int i=1;i<=K;i++) ans[++top]=r.top().id,r.pop();
        sort(ans+1,ans+top+1);
        for(int i=1;i<top;i++) printf("%d ",ans[i]);
        printf("%d",ans[top]);
    }
    putchar('\n');
    for(int i=1;i<N;i++) c.push((Node){column[i],i});
    if(L)
    {
        top=0;
        for(int i=1;i<=L;i++) ans[++top]=c.top().id,c.pop();
        sort(ans+1,ans+top+1);
        for(int i=1;i<top;i++) printf("%d ",ans[i]);
        printf("%d",ans[top]);
    }
    putchar('\n');
}
int main()
{
    freopen("seat.in","r",stdin);
    freopen("seat.out","w",stdout);
    init();
    work();
    return 0;
}

2.金明的预算方案
(budget.pas/c/cpp)
【问题描述】
金明今天很开心,家里购置的新房就要领钥匙了,新房里有一间金明自己专用的很宽
敞的房间。更让他高兴的是,妈妈昨天对他说:“你的房间需要购买哪些物品,怎么布置,
你说了算,只要不超过N 元钱就行”。今天一早,金明就开始做预算了,他把想买的物品分
为两类:主件与附件,附件是从属于某个主件的,下表就是一些主件与附件的例子:
主件 附件
电脑 打印机,扫描仪
书柜 图书
书桌 台灯,文具
工作椅 无
如果要买归类为附件的物品,必须先买该附件所属的主件。每个主件可以有0 个、1
个或2 个附件。附件不再有从属于自己的附件。金明想买的东西很多,肯定会超过妈妈限
定的N 元。于是,他把每件物品规定了一个重要度,分为5 等:用整数1~5 表示,第5 等
最重要。他还从因特网上查到了每件物品的价格(都是10 元的整数倍)。他希望在不超过N
元(可以等于N 元)的前提下,使每件物品的价格与重要度的乘积的总和最大。
设第j 件物品的价格为v[j],重要度为w[j],共选中了k 件物品,编号依次为j1,
j2,„„,jk,则所求的总和为:
v[j1]w[j1]+v[j2]*w[j2]+ „+v[jk]*w[jk]。(其中为乘号)
请你帮助金明设计一个满足要求的购物单。
【输入文件】
输入文件budget.in 的第1 行,为两个正整数,用一个空格隔开:
N m
(其中N(< 32000)表示总钱数,m(< 60)为希望购买物品的个数。)
从第2 行到第m+1 行,第j 行给出了编号为j-1 的物品的基本数据,每行有3 个非负整
4
3
4
2
1

v p q
(其中v表示该物品的价格(v<10000),p表示该物品的重要度(1~5),q表示该物品是主件还是附件。如果q=0,表示该物品为主件,如果q>0,表示该物品为附件,q是所属主件的编号)
【输出文件】
输出文件budget.out只有一个正整数,为不超过总钱数的物品的价格与重要度乘积的总和的最大值(<200000)。
【输入样例】
1000 5
800 2 0
400 5 1
300 5 1
400 3 0
500 2 0
【输出样例】
2200


有依赖的背包嘛。。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<vector>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<string>
#include<iomanip>
#include<ctime>
#include<climits>
#include<cctype>
#include<algorithm>
#ifdef WIN32
#define AUTO "%I64d"
#else
#define AUTO "%lld"
#endif
using namespace std;
#define smax(x,tmp) x=max((x),(tmp))
#define smin(x,tmp) x=min((x),(tmp))
#define maxx(x1,x2,x3) max(max(x1,x2),x3)
#define minn(x1,x2,x3) min(min(x1,x2),x3)
const int INF=0x3f3f3f3f;
const int maxn = 3205;
const int maxm = 65;
const int D = 5;
int val[maxm][D],cost[maxm][D];
int tot[maxm];
int dp[maxn];
int n,m;
inline void init()
{
    scanf("%d%d",&n,&m); n/=10;
    for(int i=1;i<=m;i++)
    {
        int v,p,q;
        scanf("%d%d%d",&v,&p,&q); v/=10;
        if(!q) val[i][++tot[i]]=p*v,cost[i][tot[i]]=v;
        else
        {
            int now=tot[q];
            for(int i=1;i<=now;i++) val[q][++tot[q]]=val[q][i]+v*p,cost[q][tot[q]]=cost[q][i]+v;
        }
    }
}
int dynamic()
{
    for(int i=1;i<=m;i++)
        for(int j=n;j>=1;j--)
            for(int k=1;k<=tot[i];k++)
                if(j-cost[i][k]>=0) smax(dp[j],dp[j-cost[i][k]]+val[i][k]);
    return dp[n];
}
int main()
{
    freopen("budget.in","r",stdin);
    freopen("budget.out","w",stdout);
    init();
    printf("%d",dynamic()*10);
    return 0;
}

3.最优贸易
(trade.pas/c/cpp)
【问题描述】
C国有n个大城市和m条道路,每条道路连接这n个城市中的某两个城市。任意两个城市之间最多只有一条道路直接相连。这m条道路中有一部分为单向通行的道路,一部分为双向通行的道路,双向通行的道路在统计条数时也计为1条。
C国幅员辽阔,各地的资源分布情况各不相同,这就导致了同一种商品在不同城市的价格不一定相同。但是,同一种商品在同一个城市的买入价和卖出价始终是相同的。
商人阿龙来到C国旅游。当他得知同一种商品在不同城市的价格可能会不同这一信息之后,便决定在旅游的同时,利用商品在不同城市中的差价赚回一点旅费。设C国n个城市的标号从1-n,阿龙决定从1号城市出发,并最终在n号城市结束自己的旅行。在旅游的过程中,任何城市可以重复经过多次,但不要求经过所有n个城市。阿龙通过这样的贸易方式赚取旅费:他会选择一个经过的城市迈入他最喜欢的商品——水晶球,并在之后经过的另一个城市卖出这个水晶球。用赚取的差价当作旅费。由于阿龙主要是来C国旅游,他决定这个贸易只进行最多一次。当然,在赚不到差价的情况下它就无需进行贸易。
假设C国有5个大城市,城市的编号和道路连接情况如下图,单向箭头表示这条道路为单向通行。双向箭头表示这条道路为双向通行。
3
1
2
4
5
假设1~n号城市的水晶球价格分别为4,3,5,6,1。
阿龙可以选择如下一条线路:1->2->3->5,并在2号城市以3的价格买入水晶球,在3号城市以5的价格卖出水晶球,赚取的旅费数为2。
阿龙也可以选择如下一条线路:1->4->5->4->5,并在第1次到达5号城市时以1的价格买入水晶球,在第2次到达4号城市时以6的价格卖出水晶球,赚取的旅费数为5。
现在给出n个城市的水晶球价格,m条道路的信息(每条道路所连接的两个城市的编号以及该条道路的通行情况)。请你告诉阿龙,他最多能赚钱多少旅费。
【输入】
输入文件名为trade.in。
第一行包含2个正整数n和m,中间用一个空格隔开,分别表示城市的数目和道路的数目。
第二行n个正整数,每两个正整数之间用一个空格隔开,按标号顺序分别表示这n个城市的商品价格。
接下来m行,每行有3个正整数,x,y,z,每两个整数之间用一个空格隔开。如果z=1,表示这条道路是城市x到城市y之间的单向道路;如果z=2,表示这条道路为城市x和城市y之间的双向道路。
【输出】
输出文件trade.out共1行,包含1个整数,表示最多能赚取的旅费。如果没有进行贸易,则输出0。
【输出输出样例】
trade.in
trade.out
5 5
4 3 5 6 1
1 2 1
1 4 1
2 3 2
3 5 1
4 5 2
4
【数据范围】
输入数据保证1号城市可以到达n号城市。
对于10%的数据,1≤n≤6。
对于30%的数据,1≤n≤100。
对于50%的数据,不存在一条旅游路线,可以从一个城市出发,再回到这个城市。
对于100%的数据,1≤n≤100000,1≤m≤500000,1≤x,y≤n,1≤z≤2,1≤各城市水晶球价格≤100。


第一种解法,枚举买入卖出分界点,前面找最小值买入后面找最大值卖出,spfa自动处理环了。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<vector>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<string>
#include<iomanip>
#include<ctime>
#include<climits>
#include<cctype>
#include<algorithm>
#ifdef WIN32
#define AUTO "%I64d"
#else
#define AUTO "%lld"
#endif
using namespace std;
#define smax(x,tmp) x=max((x),(tmp))
#define smin(x,tmp) x=min((x),(tmp))
#define maxx(x1,x2,x3) max(max(x1,x2),x3)
#define minn(x1,x2,x3) min(min(x1,x2),x3)
const int INF=0x3f3f3f3f;
const int maxn = 100005;
const int maxm = 500005;
struct Edge
{
    int to,next;
}edgeD[maxm<<1],edgeR[maxm<<1];
int headD[maxn],headR[maxn];
int maxedgeD,maxedgeR;
inline void addedge(int u,int v,bool bi_directed)
{
    edgeD[++maxedgeD] = (Edge) { v,headD[u] };
    headD[u] = maxedgeD;
    edgeR[++maxedgeR] = (Edge) { u,headR[v] };
    headR[v] = maxedgeR;
    if(bi_directed)
    {
        edgeD[++maxedgeD] = (Edge) { u,headD[v] };
        headD[v] = maxedgeD;
        edgeR[++maxedgeR] = (Edge) { v,headR[u] };
        headR[u] = maxedgeR;
    }
}
int val[maxn];
int n,m;
inline void init()
{
    scanf("%d%d",&n,&m);
    memset(headD,-1,sizeof(headD)); maxedgeD=-1;
    memset(headR,-1,sizeof(headR)); maxedgeR=-1;
    for(int i=1;i<=n;i++) scanf("%d",val+i);
    for(int i=1;i<=m;i++)
    {
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        addedge(x,y,z-1);
    }
}
bool inque[maxn];
int mn[maxn],mx[maxn];
void spfa(int S,bool flag,int dis[maxn]) // flag==true indicates maximum
{
    queue <int> que;
    for(int i=0;i<=n+1;i++) dis[i]=flag?-INF:INF;
    que.push(S); inque[S]=true; dis[S]=val[S];
    while(!que.empty())
    {
        int u=que.front(); que.pop(); inque[u]=false;
        for(int i=flag?headR[u]:headD[u];~i;i=(flag?edgeR[i]:edgeD[i]).next)
        {
            int v=(flag?edgeR[i]:edgeD[i]).to;
            if((dis[v]>min(dis[u],val[v]) && !flag) || (dis[v]<max(dis[u],val[v]) && flag))
            {
                dis[v]=flag?max(dis[u],val[v]):min(dis[u],val[v]);
                if(inque[v]) continue;
                que.push(v); inque[v]=true;
            }
        }
    }
}
int work()
{
    spfa(1,0,mn);
    spfa(n,1,mx);
    int ret = -INF;
    for(int i=1;i<=n;i++) smax(ret,mx[i]-mn[i]);
    return ret;
}
int main()
{
    freopen("trade.in","r",stdin);
    freopen("trade.out","w",stdout);
    init();
    int ans = work();
    printf("%d",max(ans,0));
    return 0;
}

第二种解法,Tarjan缩点之后求DAG最短路,用topsort,注意判断点的可达性。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<vector>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<string>
#include<iomanip>
#include<ctime>
#include<climits>
#include<cctype>
#include<algorithm>
#ifdef WIN32
#define AUTO "%I64d"
#else
#define AUTO "%lld"
#endif
using namespace std;
#define smax(x,tmp) x=max((x),(tmp))
#define smin(x,tmp) x=min((x),(tmp))
#define maxx(x1,x2,x3) max(max(x1,x2),x3)
#define minn(x1,x2,x3) min(min(x1,x2),x3)
const int INF=0x3f3f3f3f;
const int maxn = 100005;
const int maxm = 500005;
struct Edge
{
    int to,next;
}edge[maxm<<1],edgeR[maxm<<1],edgeRR[maxm<<1];
int head[maxn],headR[maxn],headRR[maxn];
int maxedge,maxedgeR,maxedgeRR;
inline void addedge(int u,int v,int bi_directed)
{
    edge[++maxedge] = (Edge) { v,head[u] };
    head[u] = maxedge;
    if(bi_directed)
    {
        edge[++maxedge] = (Edge) { u,head[v] };
        head[v] = maxedge;
    }
}
inline void addedgeR(int u,int v)
{
    edgeR[++maxedgeR] = (Edge) { v,headR[u] };
    headR[u] = maxedgeR;
    edgeRR[++maxedgeRR] = (Edge) { u,headRR[v] };
    headRR[v] = maxedgeRR;
}
int val[maxn];
int n,m;
namespace Scc
{
    int cnt;
    int sccno[maxn];
    int in[maxn];
    bool marked[maxn];
    int mn[maxn],mx[maxn];
    int dfs_clock,dfn[maxn];
    bool insta[maxn];
    stack <int> sta;
    int dfs(int u)
    {
        sta.push(u); insta[u]=true;
        int lowu=dfn[u]=++dfs_clock;
        for(int i=head[u];~i;i=edge[i].next)
        {
            int v=edge[i].to;
            if(!dfn[v])
            {
                int lowv=dfs(v);
                smin(lowu,lowv);
            }
            else if(insta[v]) smin(lowu,dfn[v]);
        }
        if(lowu>=dfn[u])
        {
            cnt++;
            int tmp;
            do
            {
                tmp=sta.top(); sta.pop(); insta[tmp]=false;
                sccno[tmp]=cnt;
                smin(mn[cnt],val[tmp]); smax(mx[cnt],val[tmp]);
            }while(tmp^u);
        }
        return lowu;
    }
    void tarjan()
    {
        memset(mn,0x3f,sizeof(mn));
        dfs(1);
    }
}
void mark(int u)
{
    Scc::marked[u]=true;
    for(int i=headRR[u];~i;i=edgeRR[i].next)
    {
        int v=edgeRR[i].to;
        if(!Scc::marked[v]) mark(v);
    }
}
bool vis[maxn];
void get_deg(int u)
{
    vis[u]=true;
    for(int i=headR[u];~i;i=edgeR[i].next)
    {
        int v=edgeR[i].to;
        Scc::in[v]++;
        if(!vis[v]) get_deg(v);
    }
}
void init()
{
    memset(head,-1,sizeof(head)); maxedge=-1;
    memset(headR,-1,sizeof(headR)); maxedgeR=-1;
    memset(headRR,-1,sizeof(headRR)); maxedgeRR=-1;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",val+i);
    for(int i=1;i<=m;i++)
    {
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        addedge(x,y,z-1);
    }
    using namespace Scc;
    tarjan();
    for(int u=1;u<=n;u++)
        for(int i=head[u];~i;i=edge[i].next)
        {
            int v=edge[i].to;
            if(sccno[u]^sccno[v]) addedgeR(sccno[u],sccno[v]);
        }
    mark(sccno[n]);
    get_deg(sccno[1]);
}
int dis[maxn];
void spfa(int S)
{
    using namespace Scc;
    memset(dis,0x3f,sizeof(dis));
    queue <int> que;
    que.push(S); dis[S]=mn[S];
    while(!que.empty())
    {
        int u=que.front(); que.pop();
        for(int i=headR[u];~i;i=edgeR[i].next)
        {
            int v = edgeR[i].to;
            if(dis[v]>min(dis[u],mn[v]))
                dis[v]=min(dis[u],mn[v]);
            if(--in[v]==0) que.push(v);
        }
    }
}
int main()
{
    freopen("trade.in","r",stdin);
    freopen("trade.out","w",stdout);
    init();
    using namespace Scc;
    spfa(sccno[1]);
    int ans=0;
    for(int i=1;i<=cnt;i++) if(marked[i])
        smax(ans,mx[i]-dis[i]);
    printf("%d",ans);
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值