当前位置: 首页 > article >正文

[NOIP2009 提高组] 最优贸易

方法一 :分层图

#include<bits/stdc++.h>
using namespace std;
using ll = long long;
#define int long long
const int N = 1e5+10,M = 3e6+10;
const int inf = 0x3f3f3f3f;
const int mod = 1e9+7;int n,q,m;int e[M],ne[M],w[M],h[3*N],idx;
void add(int a,int b,int c){e[idx] = b,ne[idx] = h[a],w[idx] = c,h[a] = idx++;
}int cost[N];
int dist[N];
bool vis[N];void spfa()
{queue<int>q;q.push(1);memset(dist,-0x3f,sizeof dist);dist[1] = 0;vis[1] = true;while(q.size()){auto t = q.front();q.pop();vis[t] = false;for(int i=h[t];~i;i=ne[i]){int j = e[i];if(dist[j]<dist[t]+w[i]){dist[j] = dist[t]+w[i];if(!vis[j])q.push(j);}}}cout<<dist[n*3];
}void solve()
{cin>>n>>m;memset(h,-1,sizeof h);for(int i=1;i<=n;i++){cin>>cost[i];;add(i,i+n,-cost[i]);add(i+n,i+n+n,cost[i]);}while(m--){int a,b,op;cin>>a>>b>>op;add(a,b,0),add(a+n,b+n,0),add(a+n+n,b+n+n,0);if(op==2){add(b,a,0),add(b+n,a+n,0),add(b+n+n,a+n+n,0);			}}spfa();}signed main()
{ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);int _;//cin>>_;_ = 1;while(_--)solve();return 0;
}

方法二:spfa图上DP  用minn存储从1号点开始到达i点为之能搞到的最小花费,maxn表示从i号点到n号点的最大价值

#include<iostream>
#include<queue>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
using ll = long long;
#define int long long
const int N = 5e5+10,M = 2e6+10;
const int inf = 0x3f3f3f3f;
const int mod = 1e9+7;int n,q,m;int cost[N];
int dmin[N],dmax[N];
bool vis[N];
int e[M],ne[M],w[M],h[N],hs[N],idx;
void add(int h[],int a,int b){e[idx] = b,ne[idx] = h[a],h[a] = idx++;
}void spfa(int h[],int dist[],int type)
{queue<int>q;memset(vis,0,sizeof vis);if(!type){memset(dist,0x3f,4*N);dist[1] = cost[1];vis[1] = true;q.push(1);}else{memset(dist,-0x3f,4*N);dist[n] = cost[n];vis[n] = true;q.push(n);}while(q.size()){auto t = q.front();q.pop();vis[t] = false;for(int i=h[t];~i;i=ne[i]){int j = e[i];if((type==0&&dist[j]>min(dist[t],cost[j]))||(type==1&&dist[j]<max(dist[t],cost[j]))){if(!type)dist[j] = min(dist[t],cost[j]);if(type)dist[j] = max(dist[t],cost[j]);if(!vis[j]){q.push(j);vis[j] = false;}}}}}void solve()
{cin>>n>>m;for(int i=1;i<=n;i++)cin>>cost[i];memset(h,-1,sizeof h),memset(hs,-1,sizeof hs);for(int i=1;i<=m;i++){int a,b,c;cin>>a>>b>>c;add(h,a,b),add(hs,b,a);if(c==2)add(h,b,a),add(hs,a,b);}spfa(h,dmin,0);spfa(hs,dmax,1);int ans = 0;for(int i=1;i<=n;i++)ans = max(ans,dmax[i]-dmin[i]);cout<<ans;}signed main()
{ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);int _;//cin>>_;_ = 1;while(_--)solve();return 0;
}

http://www.lryc.cn/news/2415405.html

相关文章:

  • SQL连接MYSQL出现对象名无效_SQL对象名无效的解决方法
  • ASP.NET Core 初学者指南 - 基于 yoyomooc 的开源项目
  • AVI与NSX-T集成实验
  • Drupal 模块开发基本教程(二)
  • VINS-Mono代码阅读笔记(十四):posegraph的存储和加载
  • struts2.0(一)
  • 虚拟机 VMware Workstation-安装详细步骤
  • 修改grub 启动项
  • ERP选型之选型九步
  • 虚拟机VirtualBox安装MAC OS 10.12图文教程
  • VMware 虚拟机里连不上网的三种解决方案
  • windows提示“你可能是盗版软件受害者”的解决方法
  • 华创e路航固件_华创e路航地图升级工具 v1.0 官方版(图文)
  • qq音乐登录页面的html代码,musicQQ音乐协议登录源码
  • 服务器如何搭建
  • 区块链与大数据,打造智能经济(读书笔记)——井底望天
  • 儒豹搜索Android新版发布 新首页亮相
  • AOP切面用aspectjweaver.jar实现代码
  • CRM平台开发实战案例:客户关系管理的挑战与解决方案
  • 基于HTML+CSS+JavaScript制作学生网页——斗破苍穹动漫(6页) 排版整洁,内容丰富,主题鲜明...
  • 暴风影音2007全功能完美版和Symantec Norton的冲突
  • 使用Spring、Hibernate、Struts的一些错误总结(不断补充)
  • WordPress添加plugin
  • 一文读懂火山引擎A/B测试的实验类型(3)——多链接实验
  • SQL Server 2008中Analysis Services的新特性——深入SQL Server 2008
  • [渗透测试]—6.3 无线网络渗透测试工具
  • 水晶报表资源下载,不断更新
  • 几个地址。。备用
  • 使用虚拟光驱重装系统
  • 小黄鸡(小贱鸡)机器人