poj 3343 Against Mammoths

题目超长，难在读题，理解题以后就可以看出是一道二分查找+最大匹配的题目了

题意：有n个人类星球，m个外星球，每个星球上开始sh艘飞船，之后每年会生产p艘飞船，人类想要战胜所有的外星球，每个人类星球和每个外星球的距离已知（需要耗费年），比如在a年初，人类星球H向外星球A宣战，那么人类会带上sh+p*a艘飞船进攻，经过k年后到达外星球，这时外星球飞船的数量为sh'+p'*(a+k)，如果这时人类飞船的数量>=外星球飞船的数量，那么就会战胜这个外星球，每个人类星球只能攻击一个外星球，一个外星球只能被一个人类星球攻击，问最少需要多少年人类才能全部战胜所有的外星球。

这题要注意的就是一些初始化的问题，我在这上面WA了很久。

#include<iostream>
using namespace std;
const int MAXN = 40001;
struct edge
{
	int v,next;
}vetex[300*300];
__int64 map[300][300];
int head[600],k;
bool vis[300];
int result[300];
struct Space_Ship
{
	__int64 inital;
	__int64 raise;
}Human[251],Alien[251];
int H,A;__int64 ans;int tmp;
void add(int u,int v)
{
	vetex[k].v = v;
    vetex[k].next = head[u];
    head[u] = k;
    ++k;
}
bool check(__int64 limit,__int64 dis,Space_Ship h,Space_Ship a)
{
	if(limit < dis)
        return false;
	if(h.inital >= a.inital+dis*a.raise)
        return true;
	if(h.inital+(limit-dis)*h.raise >= a.inital+limit*a.raise)
        return true;
    return false;
}
void build(__int64 time)
{
	k=1;
	memset(result,0,sizeof(result));
	memset(head,0,sizeof(head));
	for(int i=1;i!=H+1;i++)
		for(int j=1;j!=A+1;j++)
		{
			if(check(time,map[i][j],Human[i],Alien[j]))
				add(i,j);
		}
}
bool refresh(int i)
{
	for(int k=head[i];k;k=vetex[k].next)
	{
		int j = vetex[k].v;
		if(!vis[j])
		{
			vis[j] = true;
			if(result[j]==0 || refresh(result[j]))
			{				
				result[j] = i;
				return true;
			}	
		}
	}
	return false;
}
void solve()
{
	__int64 l_time(0);
	__int64 h_time((1<<31)-1);
	__int64 m_time;
	tmp = -1;
	while(l_time<=h_time)
	{
		m_time =(__int64)(l_time+h_time)>>1;
		build(m_time);
		ans = 0;
		for(int i=1;i!=H+1;i++)
		{
			memset(vis,false,sizeof(vis));
			if(refresh(i))
				ans++;
		}
		if(ans==A)
		{
			tmp = m_time;
			
            h_time = m_time - 1;
		}
		else
			l_time = m_time + 1;
	}
	if(tmp == -1)
		cout<<"IMPOSSIBLE"<<endl;
	else
		printf("%d\n",tmp);
}
int main()
{
	while(cin>>H>>A&&(H!=0&&A!=0))
	{
		memset(map,0,sizeof(map));
		memset(vetex,0,sizeof(vetex));
		for(int i=1;i!=H+1;i++)
			cin>>Human[i].inital>>Human[i].raise;
		for(int i=1;i!=A+1;i++)
			cin>>Alien[i].inital>>Alien[i].raise;
		for(int i=1;i!=H+1;i++)
			for(int j=1;j!=A+1;j++)
				cin>>map[i][j];
		solve();
	}
	return 0;
}