题目
给定一个等差数列,第一项是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.