[NOIP模拟赛]等差数列

题目描述
麦克在第一张纸上写下了N个十进制实数。然后，对于每个实数，他都写下了一个从0开始，以该数为公差的等差数列。比如当前数为x，则他写的等差数列为0,x,2x,3x,……接下来，他把第一张纸上的所有在区间[A,B]中的数挑选出来，剔除掉重复的，按升序写在第二张纸上。但第二天，他把第一张纸弄丢了。现在，他请你根据第二张纸还原出第一张纸上的内容。

输入

第一行包含一个自然数K,(K<=50),表示在第二张纸上有K个数。这些数都是在区间[A,B]之间的。第二行两个整数A和B（1<=A<B<=10 6)接下来K行，每行一个实数，按升序排列。每个数最多5位小数。

输出

N行，每行一个数。表示第一张纸上最开始的N个数。按升序排列。如果有多种可能，则输出个数最少的一种,如果还有多种，任意输出一种即可。

样例输入

4
1 2
1
1.4
1.5
2

样例输出

0.5

0.7

题解：

因为最多5位小数，就可以先把每个数乘上100000，输出答案的时候再除回来，方便计算。

首先寻找所有可能合法的差值。然后逐个枚举，看在当前差值下是否可以组成合法数列。

最后检验每个差值是否多余，删去多余的差值。

详见代码。

#include<cstdio>
#include<algorithm>
#include<map>
using namespace std;
typedef long long LL;
const int N=55;
const int M=1500;

int n;
LL A, B, num[N];

void Getin( LL &shu ) {
	char c; int f=5; shu=0;
	for( c=getchar(); c<'0' || c>'9'; c=getchar() );
	for( ; c>='0' && c<='9'; c=getchar() ) shu=shu*10+c-'0';
	if( c=='.' ) for( c=getchar(); c>='0' && c<='9'; c=getchar() ) shu=shu*10+c-'0', f--;
	while( f-- ) shu*=10;
}

map< LL, int >id;
LL d[M], dcnt, ans[N], acnt;
int way[N], cnt[N];

void Find_way( LL d ) {
	LL l=A, wcnt=0;
	if( A%d ) l=(A/d+1)*d;
	for( LL v=l; v<=B; v+=d )//寻找可能的数列
		if( !id[v] ) return;//数列不合法
		else way[++wcnt]=id[v];
	
	bool flg=0;
	for( int i=1; i<=wcnt; i++ )//只要该数列的一个数未在之前的数列中出现过, 该数列有存在的意义
		if( !cnt[ way[i] ] ) { flg=1; break; }
	if( !flg ) return;
	
	ans[++acnt]=d;
	for( int i=1; i<=wcnt; i++ ) cnt[ way[i] ]++;
}

void Judge( int i, LL d ) {
	LL l=A;
	if( A%d ) l=(A/d+1)*d;
	
	bool flg=0;
	for( LL v=l; v<=B; v+=d )//若该数列的所有数都在之前出现过不止一次, 该数列没有存在的意义
		if( cnt[ id[v] ]==1 ) { flg=1; break; }
	if( !flg ) ans[i]=0;
}

int main() {
	scanf( "%d", &n );
	Getin(A); Getin(B);
	
	for( int i=1; i<=n; i++ ) {
		Getin( num[i] ); id[ num[i] ]=i;
		for( int j=0; j<i; j++ )
			d[++dcnt]=num[i]-num[j];//找可能的差值
	}
	sort( d+1, d+dcnt+1 );
	dcnt=unique( d+1, d+dcnt+1 )-d-1;
	
	for( int i=1; i<=dcnt; i++ ) {
		bool flg=0;
		for( int j=1; j<=acnt; j++ )//该差值是之前某一合法差值的的倍数, 肯定不优
			if( !( d[i]%ans[j] ) ) { flg=1; break; }
		if( !flg ) Find_way( d[i] );
	}
	
	for( int i=1; i<=acnt; i++ ) Judge( i, ans[i] );
	
	for( int i=1; i<=acnt; i++ ) if( ans[i] )
		printf( "%.5lf\n", 1.0*ans[i]/1e5 );
	return 0;
}