计蒜客 2017 NOIP 提高组模拟赛（一）Day1

于 2017-09-12 00:35:47 发布
阅读量387 收藏
点赞数
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/w_h_h_/article/details/77940486
版权
T1:蒜头君打地鼠
当时用树状数组求和的，效率并不是很高，但也能过一部分数据，其实暴力用前缀和优化下似乎可以AC。

正解是矩阵旋转+二位前缀和。

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#include<vector>
#define MAXN 2005
#define pii pair<int,int>
using namespace std;
int n,k;
int dat[MAXN*2][MAXN*2];
vector<pii> a,b;
void add(int x,int y){
	for(int i=x;i<=n;i+=(i&-i)){
		for(int j=y;j<=n;j+=(j&-j)){
			dat[i][j]+=1;
		}
	}
}
int sum(int x,int y){
	int ret=0;
	for(int i=x;i>=1;i-=(i&-i)){
		for(int j=y;j>=1;j-=(j&-j)){
			ret+=dat[i][j];
		}
	}
	return ret;
}
int f(int x1,int y1,int x2,int y2){
	return sum(x2,y2)-sum(x1-1,y2)-sum(x2,y1-1)+sum(x1-1,y1-1);
}
int main()
{
//	freopen("data.in","r",stdin);
//	freopen("my.out","w",stdout);
	scanf("%d%d",&n,&k);
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			int temp;
			scanf("%d",&temp);
			if(temp){
				a.push_back(make_pair(i-j+n,i+j-1));
			}
			b.push_back(make_pair(i-j+n,i+j-1));
		}
	}
	n=n*2-1;
	k--;
	for(int i=0;i<a.size();i++){
		add(a[i].first,a[i].second);
	}
	int ans=0;
	for(int i=0;i<b.size();i++){
		int x=b[i].first,y=b[i].second;
		ans=max(ans,f(max(1,x-k),max(1,y-k),min(n,x+k),min(n,y+k)));
	}
	printf("%d\n",ans);
	return 0;
}
T2:蒜头君的树
对于一条边w=(u,v)
求和累加的次数可以这样算：
把u,v切成两个树Tu,Tv，
对于Tu中每个节点，要想求到Tu上另一个节点的距离不需要经过w（显然啊）
要想求到Tv上任何一个节点的距离都必须经过w
那么Tu的所有节点分别到Tv的所有节点必然经过w的次数为：size(Tu)*size(Tv)
然后水水就过了~
 
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#define MAXN 100005
#define ll long long
using namespace std;
ll d[MAXN];
int fa[MAXN],size[MAXN];
int cnt,first[MAXN],Next[MAXN],to[MAXN];
int n;
void Add(int u,int v){
	Next[++cnt]=first[u]; first[u]=cnt; to[cnt]=v;
	//single edge
}
void dfs(int u){
	size[u]=1;
	for(int e=first[u];e;e=Next[e]){
		int v=to[e];
		if(v==fa[u]) continue;
		dfs(v);
		size[u]+=size[v];
	}
}
int main()
{
	scanf("%d",&n);
	for(int i=2;i<=n;i++){
		scanf("%d%lld",&fa[i],&d[i]);
		Add(fa[i],i);
	}
	dfs(1);
	ll ans=0;
	for(int i=2;i<=n;i++){
		ans=(ans+d[i]*(size[1]-size[i])*size[i]);
	}
	printf("%lld\n",ans);
	int m;
	scanf("%d",&m);
	for(int i=1;i<=m;i++){
		int a; ll b;
		scanf("%d%lld",&a,&b);
		ans=(ans+(b-d[a])*(size[1]-size[a])*size[a]);
		d[a]=b;
		printf("%lld\n",ans);
	}
	return 0;
}
T3:蒜头君的坐骑
这题我用的方法是做成状态然后SPFA，考虑状态比较多，然后把其中“攻击力”直接做成值，和伤害构成一个pair，first是伤害，second是攻击力。因为伤害相同的情况下，攻击力越高越好。然后转移即可。
（其实这种写法感觉怪怪的ToT） 

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#include<queue>
#define MAXN 505
#define INF 0x7f7f7f7f
#define pii pair<int,int>
using namespace std;
int n,m;
int a[MAXN][MAXN];
int go[3]={0,1,0};
int atck,HP,mp,chi_xu;
struct Node{
	int first,second;
	Node(int p1=INF,int p2=0){
		first=p1,second=p2;
	}
	friend bool operator < (const Node &p1,const Node &p2){
		if(p1.first!=p2.first){
			return (p1.first<p2.first);
		}
		else{
			return (p1.second>p2.second);
		}
	}
	friend bool operator > (const Node &p1,const Node &p2){
		return !(p1<p2);
	}
};
Node d[MAXN][MAXN][11][6];
bool b[MAXN][MAXN][11][6];
//         number_use chixu_time
queue<pii > q1;
queue<int> q2,q3;
int getatc(int myatc,int youratc){
	if(!(HP%myatc)){
		return (HP/myatc-1)*youratc;
	}
	else{
		return (HP/myatc)*youratc;
	}
}
int SPFA(){
	d[1][1][mp][0]=Node(0,atck);
	b[1][1][mp][0]=1;
	q1.push(make_pair(1,1));
	q2.push(mp); q3.push(0);
	while(!q2.empty()){
		int x=q1.front().first,y=q1.front().second;
		int MP=q2.front(),chixu=q3.front();
		int dis=d[x][y][MP][chixu].first;
		int atc=d[x][y][MP][chixu].second;
		b[x][y][MP][chixu]=0;
		q1.pop(); q2.pop(); q3.pop();
		if(MP&&!chixu){
			int dMP=MP-1;
			int dchixu=chi_xu;
			if(d[x][y][dMP][dchixu]>Node(dis,atc)){
				d[x][y][dMP][dchixu]=Node(dis,atc);
				if(!b[x][y][dMP][dchixu]){
					b[x][y][dMP][dchixu]=1;
					q1.push(make_pair(x,y));
					q2.push(dMP); q3.push(dchixu);
				}
			}						
		}
		for(int k=0;k<2;k++){
			int dx=x+go[k],dy=y+go[k+1];
			if(1<=dx&&dx<=n&&1<=dy&&dy<=m){
				if(chixu){
					int dchixu=chixu-1;
					int datc=(dchixu==0)?atck:(atc+a[dx][dy]);
					int w=getatc(atc+a[dx][dy],a[dx][dy]);
					if(d[dx][dy][MP][dchixu]>Node(dis+w,datc)){
						d[dx][dy][MP][dchixu]=Node(dis+w,datc);
						if(!b[dx][dy][MP][dchixu]){
							b[dx][dy][MP][dchixu]=1;
							q1.push(make_pair(dx,dy));
							q2.push(MP); q3.push(dchixu);
						}
					}
				}
				else{
					int w=getatc(atc,a[dx][dy]);
					if(d[dx][dy][MP][chixu]>Node(dis+w,atc)){
						d[dx][dy][MP][chixu]=Node(dis+w,atc);
						if(!b[dx][dy][MP][chixu]){
							b[dx][dy][MP][chixu]=1;
							q1.push(make_pair(dx,dy));
							q2.push(MP); q3.push(chixu);
						}
					}
				}
			}
		}
	}
	int ret=INF;
	for(int i=0;i<=mp;i++){
		for(int j=0;j<=chi_xu;j++){
			ret=min(ret,d[n][m][i][j].first);
		}
	}
	return ret;
}
int main()
{
//	freopen("data.in","r",stdin);
	scanf("%d%d%d%d%d%d",&n,&m,&mp,&chi_xu,&HP,&atck);
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
			scanf("%d",&a[i][j]);
		}
	}
	printf("%d\n",SPFA());
	return 0;
}
然而正解并没有这么长，正解是直接三维状态f[i][j][k],表示(i,j)的位置还可以用k次技能

然后直接配合dfs暴力转移，因为只能向右和向下，所以还算好写（我的代码没有抓住这一点）

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#define MAXN 505
#define pii pair<int,int>
#define INF 0x7f7f7f7f
using namespace std;
int f[MAXN][MAXN][12];
int a[MAXN][MAXN];
int n,m;
int lasts,HP,atk;
int injure(int myatc,int youratc){
	if(!(HP%myatc)){
		return ((HP/myatc)-1)*youratc;
	}
	else{
		return (HP/myatc)*youratc;
	}
}
void dfs(int x,int y,int p,int last,int atc,int w){
	atc+=a[x][y];
	w+=injure(atc,a[x][y]);
	last--;
	if(!last){
		f[x][y][p+1]=min(f[x][y][p+1],w);
		return;
	}
	if(x+1<=n){
		dfs(x+1,y,p,last,atc,w);
	}
	if(y+1<=m){
		dfs(x,y+1,p,last,atc,w);
	}
}
void debug(int x){
	printf("======debug===========\n");
	printf("%d:\n",x);
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
			printf("%d ",f[i][j][x]);
		}
		printf("\n");
	}
	printf("\n\n");
}
int main()
{
//	freopen("data.in","r",stdin);
	memset(f,0x7f,sizeof(f));
	int mp;
	scanf("%d%d%d%d%d%d",&n,&m,&mp,&lasts,&HP,&atk);
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
			scanf("%d",&a[i][j]);
		}
	}
	f[1][1][0]=0;
	for(int p=0;p<=mp;p++){
		for(int i=1;i<=n;i++){
			for(int j=1;j<=m;j++){
				f[i][j][p]=min(f[i][j][p],injure(atk,a[i][j])+min(f[i-1][j][p],f[i][j-1][p]));
			}
		}
//		debug(p);
		for(int i=1;i<=n;i++){
			for(int j=1;j<=m;j++){
				if(i+1<=n){
					dfs(i+1,j,p,lasts,atk,f[i][j][p]);
				}
				if(j+1<=m){
					dfs(i,j+1,p,lasts,atk,f[i][j][p]);
				}
			}
		}	
	}
	printf("%d\n",f[n][m][mp]);
	return 0;
}


w_h_h_
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值