程序设计思维与实践 Week13（必做） 作业

A - TT 的神秘任务1

这一天，TT 遇到了一个神秘人。
神秘人给了两个数字，分别表示 n 和 k，并要求 TT 给出 k 个奇偶性相同的正整数，使得其和等于 n。
例如 n = 10，k = 3，答案可以为 [4 2 4]。
TT 觉得这个任务太简单了，不愿意做，你能帮他完成吗？
本题是SPJ

Input

第一行一个整数 T，表示数据组数，不超过 1000。
之后 T 行，每一行给出两个正整数，分别表示 n（1 ≤ n ≤ 1e9）、k（1 ≤ k ≤ 100）。

Output

如果存在这样 k 个数字，则第一行输出 “YES”，第二行输出 k 个数字。
如果不存在，则输出 “NO”。

Example

Input
8
10 3
100 4
8 7
97 2
8 8
3 10
5 3
1000000000 9
Output 
YES
4 2 4
YES
55 5 5 35
NO
NO
YES
1 1 1 1 1 1 1 1
NO
YES
3 1 1
YES
111111110 111111110 111111110 111111110 111111110 111111110 111111110 111111110 111111120

思路

k个正整数之和为n,这些正整数奇偶性相同，要求输出一组这样的数，简化问题，一共只有这两种可能：
1 1 1 … 1 1   $n-(k-1)$
2 2 2 … 2 2   $n-2\times (k-1)$
举几个例子就会发现和n,k的奇偶性相关，有这四种情况：

• n奇 k偶
  奇+奇=偶，偶+偶=偶，k个数相加总是偶数，直接NO
• n奇 k奇
  前面k-1个数之和总是偶数，根据最后一个数选一个情况即可
• n偶 k奇
  同上
• n偶 k偶
  同上

代码

#include <iostream>
#include <stdio.h>
using namespace std;

int t,n,k;
int main(){
	scanf("%d",&t);
	for(int u=0;u<t;u++){
		scanf("%d %d",&n,&k);
		if(n<k){cout<<"NO"<<endl;continue;}
		if(n%2==1 && k%2==0){cout<<"NO"<<endl;continue;}
		
		int m1=n-(k-1),m2=n-2*(k-1);
		if(m1%2==1 && m1>0){
			cout<<"YES"<<endl;
			for(int i=1;i<=(k-1);i++) cout<<"1 ";
			cout<<m1<<endl;
		}else if(m2%2==0 && m2>0){
			cout<<"YES"<<endl;
			for(int i=1;i<=(k-1);i++) cout<<"2 ";
			cout<<m2<<endl;
		}else{
			cout<<"NO"<<endl;
		}
	}
}

B - TT 的神秘任务2

在你们的帮助下，TT 轻松地完成了上一个神秘任务。
但是令人没有想到的是，几天后，TT 再次遇到了那个神秘人。
而这一次，神秘人决定加大难度，并许诺 TT，如果能够完成便给他一个奖励。
任务依旧只给了两个数字，分别表示 n 和 k，不过这一次是要求 TT 给出无法被 n 整除的第 k 大的正整数。
例如 n = 3，k = 7，则前 7 个无法被 n 整除的正整数为 [1 2 4 5 7 8 10]，答案为 10。
好奇的 TT 想要知道奖励究竟是什么，你能帮帮他吗？

Input

第一行一个整数 T，表示数据组数，不超过 1000。
之后 T 行，每一行给出两个正整数，分别表示 n（2 ≤ n ≤ 1e9）、k（1 ≤ k ≤ 1e9）。

Output

对于每一组数据，输出无法被 n 整除的第 k 大的正整数。

Example

Input 
6
3 7
4 12
2 1000000000
7 97
1000000000 1000000000
2 1
Output 
10
15
1999999999
113
1000000001
1

思路

正整数序列可以被无法被n整除的数分成不同的段，每段n-1个数，这样可以分成$\frac{k}{n-1}$个段，要找的数在段的$k\mathrm{m}\mathrm{o}\mathrm{d}n-1$位，整除的情况看作是余数是除数

代码

#include <iostream>
#include <stdio.h>
using namespace std;

int n,k,t;
int main(){
	scanf("%d",&t);
	for(int u=0;u<t;u++){
		scanf("%d %d",&n,&k);
		int part=k/(n-1);
		int mod=k%(n-1);
		if(mod==0) cout<<k+part-1<<endl;
		else cout<<k+part<<endl;
		
	}
	return 0;
}

C - TT 的奖励

在大家不辞辛劳的帮助下，TT 顺利地完成了所有的神秘任务。
神秘人很高兴，决定给 TT 一个奖励，即白日做梦之捡猫咪游戏。
捡猫咪游戏是这样的，猫咪从天上往下掉，且只会掉在 [0, 10] 范围内，具体的坐标范围如下图所示。

TT 初始站在位置五上，且每秒只能在移动不超过一米的范围内接住掉落的猫咪，如果没有接住，猫咪就会跑掉。例如，在刚开始的一秒内，TT 只能接到四、五、六这三个位置其中一个位置的猫咪。
喜爱猫咪的 TT 想要接住尽可能多的猫咪，你能帮帮他吗？

Input

多组样例。每组样例输入一个 m (0 < m < 100000)，表示有 m 只猫咪。
在接下来的 m 行中，每行有两个整数 a b (0 < b < 100000)，表示在第 b 秒的时候有一只猫咪掉落在 a 点上。
注意，同一个点上同一秒可能掉落多只猫咪。m = 0 时输入结束。

Output

输出一个整数 x，表示 TT 可能接住的最多的猫咪数。

Sample Input

6
5 1
4 1
6 1
7 2
7 2
8 3
0

思路

下一个结果取决于上一个结果，很典型的动态规划，每一次要么向左走，要么向右走，要么不动。
$dp[i][p]$：时间为i，位置为p的结果
$cnt[i][p]$：时间为i，位置为p的猫的数量
d p [ i ] [ p ] = m a x { d p [ i − 1 ] [ p − 1 ] , d p [ i − 1 ] [ p ] , d p [ i − 1 ] [ p + 1 ] } + c n t [ i ] [ p ] dp[i][p] = max\{dp[i-1][p-1],dp[i-1][p],dp[i-1][p+1] \}+cnt[i][p] dp[i][p]=max{dp[i−1][p−1],dp[i−1][p],dp[i−1][p+1]}+cnt[i][p]
可以使用滚动数组来记录每一次的上一次的结果，这样就可以将滚动数组的第一维优化掉，极大的减小了空间

	int begin=pos-1,end=pos+1;
	for(int i=1;i<=maxtime;i++){
		for(int p=begin;p<=end;p++){
			if(p==0) dp[p]=max(pre[p],pre[p+1])+cnt[i][p];
			else if(p==10) dp[p]=max(pre[p],pre[p-1])+cnt[i][p];
			else dp[p]=max(pre[p] , max(pre[p-1],pre[p+1]))+cnt[i][p];
		}
		strcpy(pre,dp);
		if(begin>0) begin--;
		if(end<10) end++;
	}

代码

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;

const int mmax=1e5+5;
int m,ans,pos,a,b;
int cnt[mmax][11],dp[mmax][11];

int main(){
	scanf("%d",&m);
	pos=5;
	int maxtime=0;
	while(m!=0){
		memset(cnt,0,sizeof cnt);
		memset(dp,0,sizeof dp);
		for(int i=0;i<m;i++){
			scanf("%d %d",&a,&b);
			cnt[b][a]++;
			maxtime=max(maxtime,b);
		}
		
		int begin=pos-1,end=pos+1;
		for(int i=1;i<=maxtime;i++){
			for(int p=begin;p<=end;p++){
				if(p==0) dp[i][p]=max(dp[i-1][p],dp[i-1][p+1])+cnt[i][p];
				else if(p==10) dp[i][p]=max(dp[i-1][p],dp[i-1][p-1])+cnt[i][p];
				else dp[i][p]=max(dp[i-1][p] , max(dp[i-1][p-1],dp[i-1][p+1]))+cnt[i][p];
			}
			if(begin>0) begin--;
			if(end<10) end++;
		}
		
		
		ans=0;
		for(int i=0;i<=10;i++){
			ans=max(ans,dp[maxtime][i]);
		}
		cout<<ans<<endl;
		scanf("%d",&m);
	}
	
	return 0;
}