10.7 数列 2391

题目

给定一个等差数列，第一项是a,从第二项开始看，每项与前一项的差都是一个定制b。如果用数学形式来表示，那么可以表示成a+b*x，其中x>=0，且都是整数。例如：a=1,b=2，那么这个等差数列就是:1,3,5,7,9……
再给定一个等比数列，第一项是c，从第二项开是，每项是前一项的d倍。如果用数学形式来表示等比数列，则是c*(d^y)。其中y>=0，且都是整数。例如:c=2,d=3，那么这个等比数列就是：2,6,18,54……
你的任务是计算在1至upperbound内的正整数，有多少正整数是合法的？
所谓“合法”是指：该整数属于上面给定的等差数列的某项或者属于等比数列的某项，或者既属于等差数列的项也属于等比数列的项。
1<=a,b,c,upperbound<=10^12,1<=d<=10^5

题解

因为等比数列所包含的数较少，所以考虑枚举符合条件的等比数列里的数，判断是否是等差数列里的数。数在范围内为x且（x-a)mod b<>0时，该数既合法也不与等比数列里的数重合。
然后(up-a) div b+1计算等差数列的量
注意合法的数是在范围内的

代码

var
  a,b,c,d,up,ans,t,i,j:int64;
  f:array[1..50]of int64;
begin
  readln(a,b,c,d,up);
  i:=1;f[1]:=1;
  if d=1 then
    begin
      if (((c-a)mod b<>0)or(c-a<0))and(c<=up) then inc(ans);
    end else
  while f[i]<up do
    begin
      if ((c*f[i]-a) mod b=0)and(c*f[i]-a>=0)and(c*f[i]<=up) then dec(ans);
      if (c*f[i]<=up) then inc(ans) else break;
      inc(i);
      f[i]:=f[i-1]*d;
    end;
  if up>=a then
  ans:=ans+((up-a) div b)+1;
  writeln(ans);
end.