codeforces 147B Smile House

30 篇文章 0 订阅
26 篇文章 0 订阅

A smile house is created to raise the mood. It has n rooms. Some of the rooms are connected by doors. For each two rooms (number iand j), which are connected by a door, Petya knows their value cij — the value which is being added to his mood when he moves from room i to room j.

Petya wondered whether he can raise his mood infinitely, moving along some cycle? And if he can, then what minimum number of rooms he will need to visit during one period of a cycle?

Input

The first line contains two positive integers n and m (), where n is the number of rooms, and m is the number of doors in the Smile House. Then follows the description of the doors: m lines each containing four integers ijcij и cji(1 ≤ i, j ≤ n, i ≠ j,  - 104 ≤ cij, cji ≤ 104). It is guaranteed that no more than one door connects any two rooms. No door connects the room with itself.

Output

Print the minimum number of rooms that one needs to visit during one traverse of the cycle that can raise mood infinitely. If such cycle does not exist, print number 0.

Examples
input
4 4
1 2 -10 3
1 3 1 -10
2 4 -10 -1
3 4 0 -3
output
4
Note

Cycle is such a sequence of rooms a1a2, ..., ak, that a1 is connected with a2a2 is connected with a3, ..., ak - 1 is connected with akakis connected with a1. Some elements of the sequence can coincide, that is, the cycle should not necessarily be simple. The number of rooms in the cycle is considered as k, the sequence's length. Note that the minimum possible length equals two.

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

倍增~

用f[k][i][j]表示走2^k步以内,从i走到j的最大权值,再二分求解~


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

int n,m,x,y,v,w,f[10][301][301],ans[2][301][301];
bool b;

bool cal(int u)
{
	b=0;
	for(int i=1;i<=n;i++)
	  for(int j=1;j<=n;j++) ans[0][i][j]=i==j ? 0:-999999999;
	for(int k=0;(1<<k)<=u;k++)
	  if(u&(1<<k))
	  {
	  	b^=1;
	  	for(int i=1;i<=n;i++)
	  	  for(int j=1;j<=n;j++) ans[b][i][j]=-999999999;
	  	for(int i=1;i<=n;i++)
	  	  for(int j=1;j<=n;j++)
	  	    for(int z=1;z<=n;z++)
	  	      ans[b][i][j]=max(ans[b][i][j],ans[b^1][i][z]+f[k][z][j]);
	  }
	for(int i=1;i<=n;i++) if(ans[b][i][i]>0) return 1;
	return 0;
}

int main()
{
	scanf("%d%d",&n,&m);
	for(int k=0;(1<<k)<=n;k++)
	  for(int i=1;i<=n;i++)
	    for(int j=1;j<=n;j++) f[k][i][j]=i==j ? 0:-999999999;
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d%d%d",&x,&y,&v,&w);
		f[0][x][y]=v;
		f[0][y][x]=w;
	}
	for(int k=1;(1<<k)<=n;k++)
	  for(int i=1;i<=n;i++)
	    for(int j=1;j<=n;j++)
	      for(int z=1;z<=n;z++)
	       f[k][i][j]=max(f[k][i][j],f[k-1][i][z]+f[k-1][z][j]);
	int l=2,r=n+1;
	while(l<r)
	{
		int mid=(l+r)>>1;
		if(cal(mid)) r=mid;
		else l=mid+1;
	}
	printf("%d\n",l%(n+1));
	return 0;
}


  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值