GDKOI2006 纪念邮票 PASCAL

纪念邮票

问题描述：

邮局最近推出了一套特殊的纪念邮票，这套邮票共有N张，邮票面值各不相同，按编号顺序为1分，2分，．．．．．．，N分。

 

小杭是个集邮爱好者，他很喜欢这套邮票，可惜现在他身上只有M分，并不够把全套都买下。他希望尽量买，最好刚好花光所有钱。作为一个集邮爱好者，小杭也不想买的邮票编号断断续续。所以小杭打算买面值a分至b分的b-a+1张连续的邮票，且总价值刚好为M分。

 

你的任务是求出所有符合要求的方案，以[a,b]的形式输出。

 

输入格式(stamp.in)

输入文件只有一行，包含两个数N和M（1<=N，M<=10⁹）。

 

输出格式(stamp.out)

输出文件每行包含一个合法方案：[a,b]。按a值从小到大输出。

输出文件不含任何空格。

 

输入样例

20 15

 

输出样例

[1,5]

[4,6]

[7,8]

[15,15]

 

 

解题报告

由初等数学的知识可以知道：

a + (a+1) + ...... + (b-1) + b = (a+b)*(b-a+1)/2

（我去，这不是等差数列吗？当时没学怎么想到？)

所以 2*M = (a+b)*(b-a+1)

如果把 2*M 分解为两个数相乘: 2*M = X * Y

其中X = (a+b), Y = (b-a+1), (X>=Y)

可以推出：b=(X+Y-1)/2, a=(X-Y+1)/2, 数对(X,Y)跟(a,b)一一对应

枚举所有合法数对(X,Y)，就可以不重复地直接算出所有(a,b)

由于X>=Y，所以只需要对Y枚举1到sqr(2*M)，找到所有合法的(X,Y)

 

看起来好像不是很复杂，最主要是第一步没想到。

 

var
 m,n,x,y,a,b:longint;
begin
  read(n,m);
  for y:=trunc(sqrt(2*m))+1 downto 1 do
  if 2*m mod y=0 then
  begin
   x:=2*m div y;
   a:=trunc((x-y+1)/2);
   b:=trunc((x+y-1)/2);
   if (b<=n)and((a+b)*(b-a+1)=2*m) then
    writeln('[',a,',',b,']');
  end;

end.