USACO Section 4

前言

好久没更新这个系列了，最近闲的无聊写一下。有两题搜索懒得写了。

---

P2737 [USACO4.1]麦香牛块Beef McNuggets

https://www.luogu.com.cn/problem/P2737

解题思路

先只考虑\(a_1\)，假设我们拼出了\(w\)，那么一定能拼出\(w+ka_1\)，换句话我们考虑最大的不能拼出的\(ka_1+w(0<w<a_1)\)。

考虑到\(gcd\neq 1\)显然无解所以\(a\)中肯定存在\(gcd(a_i,a_1)=1\)，这样对于每个\(ka_i\% a_1(k<a_1)\)都会有不同的值，所以根据鸽笼原理最大不能拼出的\(ka_1+w\leq a_1a_i-a_i\)，这个范围不会很大，暴力背包就好了。

code

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int L=256*256;
int n,a[11],f[L];
int main()
{
	scanf("%d",&n);
	int d=0,flag=0;
	for(int i=1;i<=n;i++){
		scanf("%d",&a[i]);
		d=__gcd(d,a[i]);
		flag|=(a[i]==1);
	}
	if(d!=1||flag)return puts("0")&0;
	f[0]=1;
	for(int i=1;i<L;i++)
		for(int j=1;j<=n;j++)
			if(i>=a[j])f[i]|=f[i-a[j]];
	for(int i=L-1;i>=1;i--)
		if(!f[i])return printf("%d\n",i)&0;
	return 0;
}

---

P2738 [USACO4.1]篱笆回路Fence Loops

https://www.luogu.com.cn/problem/P2738

解题思路

每个篱笆开两个点，每个篱笆对之间记录他们用来连接的端点然后建图就变成最小环问题了。

枚举篱笆断边然后两端跑最短路就好了。

code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int N=210;
struct node{
	int to,next,w;
	bool ban;
}a[N*N];
int n,tot,ls[N],f[N],p[N][N],rev[N],pos[N],ans;
bool v[N];queue<int> q;
void addl(int x,int y,int w){
	a[++tot].to=y;
	a[tot].next=ls[x];
	ls[x]=tot;a[tot].w=w;
	return;
}
void SPFA(int s){
	memset(f,0x3f,sizeof(f));
	q.push(s);v[s]=1;f[s]=0;
	while(!q.empty()){
		int x=q.front();q.pop();v[x]=0;
		for(int i=ls[x];i;i=a[i].next){
			int y=a[i].to;
			if(a[i].ban)continue;
			if(f[x]+a[i].w<f[y]){
				f[y]=f[x]+a[i].w;
				if(!v[y])q.push(y);
			}
		}
	}
	return;
}
int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++){
		int x,l,s1,s2;
		scanf("%d%d%d%d",&pos[i],&l,&s1,&s2);
		rev[pos[i]]=i;
		addl(i*2-1,i*2,l);addl(i*2,i*2-1,l);
		for(int j=0;j<s1;j++)
			scanf("%d",&x),p[i][x]=i*2-1;
		for(int j=0;j<s2;j++)
			scanf("%d",&x),p[i][x]=i*2;
	}
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			if(p[i][j])addl(p[i][j],p[rev[j]][pos[i]],0);
	ans=2147483647;
	for(int i=1;i<=n;i++){
		int w=i*2;
		a[w].ban=a[w-1].ban=1;
		SPFA(w-1);ans=min(ans,a[w].w+f[w]);
		a[w].ban=a[w-1].ban=0;
	}
	printf("%d\n",ans);
	return 0;
}

---

P2740 [USACO4.2]草地排水Drainage Ditches

https://www.luogu.com.cn/problem/P2740

解题思路

最大流模板

code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int N=210;
struct node{
	int to,next,w;
}a[N<<1];
int n,m,tot=1,ls[N],dep[N],ans;
queue<int> q;
void addl(int x,int y,int w){
	a[++tot].to=y;a[tot].next=ls[x];ls[x]=tot;a[tot].w=w;
	a[++tot].to=x;a[tot].next=ls[y];ls[y]=tot;a[tot].w=0;
}
bool bfs(){
	memset(dep,0,sizeof(dep));dep[1]=1;
	while(!q.empty())q.pop();q.push(1);
	while(!q.empty()){
		int x=q.front();q.pop();
		for(int i=ls[x];i;i=a[i].next){
			int y=a[i].to;
			if(!a[i].w||dep[y])continue;
			dep[y]=dep[x]+1;
			if(y==m)return 1;
			q.push(y);
		}
	}
	return 0;
}
int dinic(int x,int flow){
	if(x==m)return flow;
	int rest=0,k;
	for(int i=ls[x];i;i=a[i].next){
		int y=a[i].to;
		if(dep[x]+1!=dep[y]||!a[i].w)continue;
		rest+=(k=dinic(y,min(flow-rest,a[i].w)));
		a[i].w-=k;a[i^1].w+=k;
		if(flow==rest)return flow;
	}
	if(!rest)dep[x]=0;
	return rest;
}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++){
		int x,y,w;
		scanf("%d%d%d",&x,&y,&w);
		addl(x,y,w);
	}
	while(bfs())
		ans+=dinic(1,2147483647);
	printf("%d\n",ans);
	return 0;
}

---

P1894 [USACO4.2]完美的牛栏The Perfect Stall

https://www.luogu.com.cn/problem/P1894

解题思路

最大匹配，把上面那个最大流改一下就好了。

code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int N=410;
struct node{
	int to,next,w;
}a[N*N];
int n,m,s,t,tot=1,ls[N],dep[N],ans;
queue<int> q;
void addl(int x,int y,int w){
	a[++tot].to=y;a[tot].next=ls[x];ls[x]=tot;a[tot].w=w;
	a[++tot].to=x;a[tot].next=ls[y];ls[y]=tot;a[tot].w=0;
}
bool bfs(){
	memset(dep,0,sizeof(dep));dep[s]=1;
	while(!q.empty())q.pop();q.push(s);
	while(!q.empty()){
		int x=q.front();q.pop();
		for(int i=ls[x];i;i=a[i].next){
			int y=a[i].to;
			if(!a[i].w||dep[y])continue;
			dep[y]=dep[x]+1;
			if(y==t)return 1;
			q.push(y);
		}
	}
	return 0;
}
int dinic(int x,int flow){
	if(x==t)return flow;
	int rest=0,k;
	for(int i=ls[x];i;i=a[i].next){
		int y=a[i].to;
		if(dep[x]+1!=dep[y]||!a[i].w)continue;
		rest+=(k=dinic(y,min(flow-rest,a[i].w)));
		a[i].w-=k;a[i^1].w+=k;
		if(flow==rest)return flow;
	}
	if(!rest)dep[x]=0;
	return rest;
}
int main()
{
	scanf("%d%d",&n,&m);s=n+m+1;t=s+1;
	for(int i=1;i<=n;i++){
		int k,x;scanf("%d",&k);
		for(int j=0;j<k;j++)
			scanf("%d",&x),addl(i,x+n,1);
	}
	for(int i=1;i<=n;i++)addl(s,i,1);
	for(int i=1;i<=m;i++)addl(i+n,t,1);
	while(bfs())
		ans+=dinic(s,2147483647);
	printf("%d\n",ans);
	return 0;
}

---

P2751 [USACO4.2]工序安排Job Processing

https://www.luogu.com.cn/problem/P2751

解题思路

考虑贪心，第一个正着做，开一个单调队列然后每次取最早做完的。

第二个和第一个一样反着做，然后两道工序时间长的配时间短的就好了。

code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#define mp(x,y) make_pair(x,y)
using namespace std;
const int N=1100;
int n,A,B,t[N],s[N],a[N],ans;
priority_queue<pair<int,int> > q;
int main()
{
	scanf("%d%d%d",&n,&A,&B);
	for(int i=1;i<=A;i++)scanf("%d",&a[i]),q.push(mp(-a[i],i));
	for(int i=1;i<=n;i++){
		pair<int,int> k=q.top();
		t[i]=-k.first;q.pop();
		q.push(mp(k.first-a[k.second],k.second));
	}
	printf("%d ",t[n]);
	while(!q.empty())q.pop();
	for(int i=1;i<=B;i++)scanf("%d",&a[i]),q.push(mp(-a[i],i));
	for(int i=1;i<=n;i++){
		pair<int,int> k=q.top();
		s[i]=-k.first;q.pop();
		q.push(mp(k.first-a[k.second],k.second));
	}
	for(int i=1;i<=n;i++)
		ans=max(ans,s[i]+t[n-i+1]);
	printf("%d\n",ans);
	return 0;
}

---

P2687 [USACO4.3]逢低吸纳Buy Low, Buy Lower

https://www.luogu.com.cn/problem/P2687

解题思路

最长下降子序列，但是因为是股价序列不同，所以方案数转移的时候枚举到上一个与它相等的就要停止。

要开\(long\ double\)或者高精度

code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cctype>
#define ll long double
using namespace std;
const int N=5100;
long long n,a[N];
ll f[N],g[N];
signed main()
{
	scanf("%lld",&n);
	for(int i=1;i<=n;i++)
		scanf("%lld",&a[i]);
	a[++n]=0;a[0]=1e18;f[0]=g[0]=1;
	for(int i=1;i<=n;i++){
		for(int j=0;j<i;j++)
			if(a[j]>a[i])f[i]=max(f[i],f[j]+1);
		for(int j=i-1;j>=0;j--)
			if(a[j]==a[i])break;
			else if(a[j]>a[i]&&f[j]+1==f[i])
				g[i]+=g[j];
	}
	printf("%.0Lf %.0Lf\n",f[n]-2,g[n]);
	return 0;
}

---

P2752 [USACO4.3]街道赛跑Street Race

https://www.luogu.com.cn/problem/P2752

解题思路

先跑一遍\(Flody\)，然后枚举点，每次重新跑\(Flody\)判必经点。

然后如果是必经点用第一次的\(Flody\)判中间点。

code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
const int N=51;
int n;bool a[N][N],f[N][N],g[N][N];
vector<int>ans1,ans2;
int main()
{
	while(1){
		int x;scanf("%d",&x);
		if(x==-1)break;
		while(x!=-2){
			a[n][x]=g[n][x]=1;
			scanf("%d",&x);
		}
		n++;
	}
	for(int i=0;i<n;i++)a[i][i]=g[i][i]=1;
	for(int k=0;k<n;k++)
		for(int i=0;i<n;i++)
			for(int j=0;j<n;j++)
				g[i][j]|=g[i][k]&g[k][j];
	for(int p=1;p<n-1;p++){
		memset(f,0,sizeof(f));
		for(int i=0;i<n;i++)
			for(int j=0;j<n;j++)
				if(i!=p&&j!=p)f[i][j]=a[i][j];
		for(int k=0;k<n;k++)
			for(int i=0;i<n;i++)
				for(int j=0;j<n;j++)
					f[i][j]|=f[i][k]&f[k][j];
		if(!f[0][n-1]){
			ans1.push_back(p);
			bool flag=1;
			for(int i=0;i<n;i++){
				if(i==p)continue;
				if(f[0][i])flag&=(g[i][p]&&!g[p][i]);
				if(f[i][n-1])flag&=g[p][i];
			}
			if(flag)ans2.push_back(p);
		}
	}
	printf("%d",ans1.size());for(int i=0;i<ans1.size();i++)printf(" %d",ans1[i]);putchar('\n');
	printf("%d",ans2.size());for(int i=0;i<ans2.size();i++)printf(" %d",ans2[i]);putchar('\n');
	return 0;
}

---

P2753 [USACO4.3]字母游戏Letter Game

解题思路

先看懂题，然后发现虽然字符串很多，但是合法的字符串个数不会超过\(2^7\)所以只跑合法的不会超过\(O(n^2)\)

code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
const int w[26]={2,5,4,4,1,6,5,5,1,7,6,3,5,2,3,5,7,2,1,2,4,6,6,7,5,7};
struct node{
	char s[7];
	int l,c;
}a[410000];
int n,ans,c[30],v[30];
char t[7];
vector<pair<int,int> >pr;
//bool operator==(node x,node y){
//	for(int i=0;i<7;i++)
//		if(x.s[i]!=y.s[i])return 0;
//	return 1;
//}
//bool operator<(node x,node y){
//	for(int i=0;i<7;i++)
//		if(x.s[i]!=y.s[i])return x.s[i]<y.s[i];
//	return 0;
//}
bool check(node x,node y){
	for(int i=0;i<26;i++)v[i]=c[i];
	for(int i=0;i<x.l;i++)
		if(v[x.s[i]-'a'])v[x.s[i]-'a']--;
		else return 0;
	for(int i=0;i<y.l;i++)
		if(v[y.s[i]-'a'])v[y.s[i]-'a']--;
		else return 0;
	return 1;
}
int main()
{
	scanf("%s",t);int m=strlen(t);
	for(int i=0;i<m;i++)c[t[i]-'a']++;
	while(1){
		++n;scanf("%s",a[n].s);
		if(a[n].s[0]=='.'){n--;break;}
		a[n].l=strlen(a[n].s);
		if(!check(a[n],a[0]))
		{n--;continue;}
		for(int i=0;i<a[n].l;i++)
			a[n].c+=w[a[n].s[i]-'a'];
	}
	ans=0;
	for(int i=1;i<=n;i++){
		if(a[i].c>ans)ans=a[i].c,pr.clear();
		if(a[i].c==ans)pr.push_back(make_pair(i,0));
		for(int j=i;j<=n;j++){
			if(check(a[i],a[j])){
				if(a[i].c+a[j].c>ans)ans=a[i].c+a[j].c,pr.clear();
				if(a[i].c+a[j].c==ans)pr.push_back(make_pair(i,j));
			}
		}
	}
	printf("%d\n",ans);
	for(int i=0;i<pr.size();i++)
		printf("%s %s\n",a[pr[i].first].s,a[pr[i].second].s);
	return 0;
}

---

P1344 [USACO4.4]追查坏牛奶Pollutant Control

https://www.luogu.com.cn/problem/P1344

解题思路

以前写的了，费用流模板，那个费用变成第一关键值乘上一个很大的值加上第二个关键值就好了。

code

// luogu-judger-enable-o2
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
#define ll long long
using namespace std;
const ll N=320,M=10010,inf=1e18,cs=2333;
struct node{
	ll to,next,w;
}a[M*2];
ll tot=1,n,s,t,m,ans;
ll dep[N],ls[N];
queue<int> q;
void addl(ll x,ll y,ll w)
{
	a[++tot].to=y;a[tot].next=ls[x];ls[x]=tot;a[tot].w=w;
	a[++tot].to=x;a[tot].next=ls[y];ls[y]=tot;a[tot].w=0;
}
bool bfs()
{
	memset(dep,0,sizeof(dep));
	while(!q.empty()) q.pop();
	q.push(s);dep[s]=1;
	while(!q.empty())
	{
		ll x=q.front();q.pop();
		for(ll i=ls[x];i;i=a[i].next)
		{
			ll y=a[i].to;
			if(!dep[y]&&a[i].w){
				dep[y]=dep[x]+1;
				if(y==t) return true;
				q.push(y);
			}
		}
	}
	return false;
}
ll dinic(ll x,ll flow)
{
	ll rest=0,k;
	if(x==t) return flow;
	for(ll i=ls[x];i;i=a[i].next)
	{
		ll y=a[i].to;
		if(dep[x]+1==dep[y]&&a[i].w)
		{
			rest+=(k=dinic(y,min(a[i].w,flow-rest)));
			a[i].w-=k;a[i^1].w+=k;
			if(rest==flow) return flow;
		}
	}
	if(!rest) dep[x]=0;
	return rest;
}
void netflow()
{
	while(bfs())
	  ans+=dinic(s,inf);
}
int main()
{
	scanf("%lld%lld",&n,&m);
	s=1;t=n;
	for(ll i=1;i<=m;i++)
	{
		ll x,y,w;
		scanf("%lld%lld%lld",&x,&y,&w);
		w=w*cs+1;addl(x,y,w);
	}
	netflow();
	printf("%lld %lld",ans/cs,ans%cs);
}

---

后面两道搜索懒得写了