bzoj1021

这个题写了三个小时。。。。

太洋务了最近。。。

估计是因为那个计算几何闹得。。。

这个题大概就是2*3一共6中情况， 分别是一个人给另外两个或者是两个人给一个人，那么用ijk分别表示对于面值是i，a有j，b有k的方案数

然后用队列暴力转移就可以

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<queue>
#include<cmath>
#define inf 0x7fffffff/2
 
using namespace std;
 
int c[5],a1,a2,a3,b[5],a[5][100],tot,f[8][1009][1009];
 
struct wbysr
{
	int num,a,b;
};
 
int GCD(int x,int y)
{
	if(!y)
		return x;
	return GCD(y,x%y);
}
 
bool check(int a,int b,int c,int d)
{
	a=abs(a);
	b=abs(b);
	c=abs(c);
	if(a%d||b%d||c%d)
		return 0;
	return 1;
}
 
int main()
{
	int money[]={0,1,5,10,20,50,100};
	int gcd[7];
	gcd[6]=100;
	for(int i=5;i>=1;i--)
		gcd[i]=GCD(gcd[i+1],money[i]);
 
 
 
	scanf("%d%d%d",&a1,&a2,&a3);
	for(int i=1;i<=3;i++)
		for(int j=6;j>=1;j--)
			scanf("%d",&a[i][j]),b[i]+=a[i][j]*money[j];
	tot=b[1]+b[2]+b[3];
	c[1]=b[1]-a1+a3;
	c[2]=b[2]-a2+a1;
	c[3]=b[3]-a3+a2;
 
	for(int i=1;i<=3;i++)
		if(c[i]<0)
		{
			printf("impossible\n");
			return 0;
		}
	for(int i=1;i<=6;i++)
		for(int j=0;j<=1000;j++)
			for(int k=0;k<=1000;k++)
				f[i][j][k]=inf;
 
 
	f[6][b[1]][b[2]]=0;
	queue<wbysr>q;
	q.push((wbysr){0,b[1],b[2]});
	while(!q.empty())
	{
		//printf("------------------------------\n");
		wbysr now=q.front();
		q.pop();
		int num=now.num+1;
		if(num==7)
			continue;
		int next[5],current[5];
		current[0]=now.a;
		current[1]=now.b;
		current[2]=tot-now.a-now.b;
		//printf("%d %d %d %d    %d\n",now.num,now.a,now.b,tot-now.a-now.b,f[now.num][now.a][now.b]);
		
		//1->2
		for(int i=0;i<=2;i++)
			for(int j=0;j<=a[i+1][num];j++)
				for(int k=0;k<=a[i+1][num];k++)
				{
					if(j+k>a[i+1][num]||(!j&&!k))
						continue;
 
					next[i]=current[i]-money[num]*(j+k);
					next[(i+1)%3]=current[(i+1)%3]+money[num]*j;
					next[(i+2)%3]=current[(i+2)%3]+money[num]*k;
 
					if(check(next[0]-c[1],next[1]-c[2],next[2]-c[3],gcd[num+1]))
					{
						q.push((wbysr){num,next[0],next[1]});
						f[num][next[0]][next[1]]=min(f[num][next[0]][next[1]],f[now.num][now.a][now.b]+j+k);
					}
				}
 
		//2->1
		for(int i=0;i<=2;i++)
			for(int j=0;j<=a[(i+1)%3+1][num];j++)
				for(int k=0;k<=a[(i+2)%3+1][num];k++)
				{
					if(!j&&!k)
						continue;
 
					next[i]=current[i]+money[num]*(j+k);
					next[(i+1)%3]=current[(i+1)%3]-money[num]*j;
					next[(i+2)%3]=current[(i+2)%3]-money[num]*k;
 
					if(check(next[1]-c[2],next[0]-c[1],next[2]-c[3],gcd[num+1]))
					{
						q.push((wbysr){num,next[0],next[1]});
						f[num][next[0]][next[1]]=min(f[num][next[0]][next[1]],f[now.num][now.a][now.b]+j+k);
					}
				}
		
		f[num][current[0]][current[1]]=min(f[num][current[0]][current[1]],f[now.num][now.a][now.b]);
		q.push((wbysr){num,now.a,now.b});
	}
	if(f[6][c[1]][c[2]]!=inf)
		printf("%d\n",f[6][c[1]][c[2]]);
	else
		printf("impossible\n");
	return 0;
}