420. [SDOI2009] 晨跑
★★☆ 输入文件:run.in
输出文件:run.out
简单对比 时间限制:2 s 内存限制:128 MB
Elaxia最近迷恋上了空手道,他为自己设定了一套健身计划,比如俯卧撑、仰卧起坐等等,不过到目前为止,他坚持下来的只有晨跑。 现在给出一张学校附近的地图,这张地图中包含N个十字路口和M条街道,Elaxia只能从一个十字路口跑向另外一个十字路口,街道之间只在十字路口处相交。Elaxia每天从寝室出发跑到学校,保证寝室编号为1,学校编号为N。 Elaxia的晨跑计划是按周期(包含若干天)进行的,由于他不喜欢走重复的路线,所以在一个周期内,每天的晨跑路线都不会相交(在十字路口处),寝室和学校不算十字路口。Elaxia耐力不太好,他希望在一个周期内跑的路程尽量短,但是又希望训练周期包含的天数尽量长。 除了练空手道,Elaxia其他时间都花在了学习和找MM上面,所有他想请你帮忙为他设计一套满足他要求的晨跑计划。
输入 第一行2个数n,m。表示十字路口数和街道数。 接下来m行,每行3个数a,b,c,表示十字路口a和十字路口b之间有条长度为c的街道(单向)。
输出 2个数,第一个数为最长周期的天数,第二个数为满足最长天数的条件下最短的路程长度。
输入样例 7 10 1 2 1 1 3 1 2 4 1 3 4 1 4 5 1 4 6 1 2 5 5 3 6 6 5 7 1 6 7 1
输出样例 2 11
N<=200 M<=20000
拆点 终点不动
i->(i+n)
2~(n-1) 容量为1 (1)为mod
#include<iostream> #include<cmath> #include<cstring> #include<cstdio> #include<queue> #include<map> #include<cstdlib> #include<algorithm> #define V 505 #define mod 1000000007 #define LL long long using namespace std; int n,m,sd,T,S,rt; int ff,ss,nn; int a[V],pre[V],dep[V],q[V],dis[V],pr[V],p; int sb[V][V]; int ins[V]; int b[V],s[V]; int fd,csd; struct da {int f,to,next,dis,cast; }Edge[V*V*2]; int head[V],tot; inline void add(int x,int y,int zz,int dd) {Edge[tot].f=x;Edge[tot].to=y;Edge[tot].dis=zz;Edge[tot].cast=dd;Edge[tot].next=head[x];head[x]=tot++; Edge[tot].f=y;Edge[tot].to=x;Edge[tot].dis=0;Edge[tot].cast=-dd;Edge[tot].next=head[y];head[y]=tot++; } bool FIND(int st,int ed) {//cout<<st<<endl;memset(dis,0x3f,sizeof(dis));memset(ins,false,sizeof(ins));memset(pr,0,sizeof(pr));memset(s,0x3f,sizeof(s));int inf;inf = dis[0];queue <int> q;q.push(st), dis[st] = 0, s[st] = inf;while(!q.empty()) {int op = q.front(); q.pop();for(int i = head[op] ; i != -1 ; i = Edge[i].next) {if(Edge[i].dis>0&&dis[Edge[i].to]>dis[op]+Edge[i].cast) {dis[Edge[i].to] = dis[op]+Edge[i].cast;pr[Edge[i].to] = i;s[Edge[i].to] = min(s[op],Edge[i].dis);//printf("i &&& %d to %d %d\n",i,Edge[i].to,dis[Edge[i].to]);//cout<<dis[Edge[i].to]<<" "<<Edge[i].to<<" "<<ed<<endl;if(!ins[Edge[i].to]) {ins[Edge[i].to] = true;q.push(Edge[i].to);}}}ins[op] = false;} //cout<<Cost<<endl;if(dis[ed]==inf) return false;fd += s[ed];csd += dis[ed]*s[ed];int w = ed;while(w!=st) {Edge[pr[w]].dis -= s[ed];Edge[pr[w]^1].dis += s[ed];w = Edge[pr[w]].f;}//cout<<Cost<<" $%"<<endl;return true;} int st[V][V]; inline int haha() {// freopen("in.txt","r",stdin);freopen("out.txt","w",stdout);freopen("run.in","r",stdin); freopen("run.out","w",stdout);memset(head,-1,sizeof(head));cin>>n>>m;S=1;T=n;int x,y,z;for(int i=1;i<=m;i++){cin>>x>>y>>z;add(x+n,y,1,z); }add(1,1+n,mod,0);for(int i=2;i<n;i++) add(i,i+n,1,0);while(FIND(S,T));cout<<fd<<" "<<csd<<endl; return 0; } int gg=haha(); int main() {;}