3032. 是坚挺数吗？

单点时限: 2.0 sec

内存限制: 256 MB

考虑如下自然数序列：

1 2 3 4 5 6 7 8 9 10 11 12 …
隔一个数划去一个数，形成如下序列：

1 3 5 7 9 11 13 15 17 …
再隔两个数划去一个数，形成新的序列：

1 3 7 9 13 15 19 21 …
再隔三个数划去一个数，形成新的序列：

1 3 7 13 15 19 25 27 …
这个过程无限地进行下去。无限次地进行后，依然未被划去的数被称为 坚挺数。

编写程序：判断一个给定的整数，是否是坚挺数。

输入格式
第 1 行：一个整数 T (1≤T≤10) 为问题数。

接下来共 T 行整数，每行输入一个正整数 N（1≤N≤10000）。

输出格式
对于每个问题，输出一行问题的编号（0 开始编号，格式：case #0: 等）。

然后对应每个问题在一行中输出整数 N 是否为坚挺数（Yes或者 No）。若为 Yes，则输出该整数在所有坚挺数中的位次（按从小到大的顺序排列）。注意：中间要输出一个空格。

样例
input
3
19
148
4783
output
case #0:
Yes 5
case #1:
No
case #2:
Yes 78

/*
思想：用两个数组模拟。
vector中无find函数，使用find(first,last,value)在algorithm库中，如果寻找不到返回v.end(),查找成功返回迭代器。
*/
#include<iostream>
#include<vector>
#include<algorithm>
using namespace std;
int main() {
	vector<int>v1,v2;
	int index=2;
	for(int i = 1; i < 10000; i+=2)
		v1.push_back(i);
	while(1) {
		v2.clear();
		for(int i = 0; i < v1.size(); i++) {
			if((i+1)%(index+1)) {
				v2.push_back(v1[i]);
			}
		}
		index++;
		if(v1==v2)
			break;
		else
			v1=v2;
	}
	int t;
	cin>>t;
	for(int i = 0; i < t; i++)
	{
		int m;
		cin>>m;
		printf("case #%d:\n",i);
		vector<int>::iterator ans=find(v2.begin(),v2.end(),m);
		if(ans!=v2.end())
		cout<<"Yes"<<" "<<++ans-v2.begin()<<endl;
		else
		cout<<"No"<<endl;
	}
	return 0;
}

思路二：
设数字当前的位置是 n，当前是第 k 轮淘汰，那么有如下三种情况：
如果 n<k，那么 n 再也不可能被淘汰了，所以它是一个坚挺数。
如果 n%k=0，那么 n 就在这一轮被淘汰了，所以它不够坚挺。
如果前两种情况都不能判断，说明 n 暂时幸存了下来，它的新位置是 n−⌊n/k⌋，这是因为在它前面有 ⌊n/k⌋ 个数在这一轮中被淘汰掉了。继续进行下一轮，直到得出结论。

#include<iostream>
#include<vector>
#include<algorithm>
using namespace std;
int is_firm(int n,int k) {
	if(n<k)
		return n;
	else if(n%k==0)
		return -1;
	else {
		return is_firm(n-n/k,k+1);
	}
}
int main() {
	int t;
	cin>>t;
	for(int i = 0; i < t; i++)
	{
		int m;
		cin>>m;
		printf("case #%d:\n",i);
		if(is_firm(m,2)!=-1)
		cout<<"Yes"<<" "<<is_firm(m,2)<<endl;
		else
		cout<<"No"<<endl;
	}
	return 0;

}