题1

A  A×B

Time Limit:1000MS  Memory Limit:65535K

题型: 编程题   语言: 无限制

 

描述

给出两个数列a、b，构造一个矩阵mat，使得mat[i][j] = a[i] * b[j]。
首先，需要求出一个数列cnt，cnt[x]表示能被x整除的mat[i][j]的个数，x∈[1, max{mat[i][j]}]。
然后，请对数列cnt求和，并输出结果。

 

 

输入格式

输入T，表示以下有T组测试数据：

对于每一组测试数据，输入一行9个整数 n, m, a0, b0, p, q, x, y, mx。
其中，a[0] = a0，a[i] = a[i - 1] * p + x (1 <= i < n)；
      b[0] = b0, b[j] = b[j - 1] * q + y (1 <= j < m)。
由于生成出来的数字比较大，所以对于a和b每一个元素对其模mx。

数据范围：
T <= 5
1 <= n, m <= 2000
1 <= a0, b0, p, q, x, y, mx <= 1000

 

输出格式

每组数据一行，输出对应的结果。

 

输入样例

1
2 2 1 2 3 4 5 7 10

 

输出样例

17

 

Hint

解释sample：
a[] = { 1, 8 };
b[] = { 2, 5 };
mat[][] = {
    { 2, 16 },
    { 5, 40 }
};
除数： 1 2 4 5 8 10 16 20 40
个数： 4 3 2 2 2  1  1  1  1

 

#include<stdio.h>
#include<string.h>
int a[2000],b[2000];
long long count[1000000],cc[4000000];
int main()
{
    int i,j,amax,bmax;
    int n,m,a0,b0,p,q,x,y,mx;
    for(i=1;i<1000000;i++)\\分别计算前1000000的因子个数
        for(j=i;j<1000000;j+=i)\\i的倍数，那么i也是这个倍数的一个因子
        cc[j]++;
    int T,t;
    scanf("%d",&T);
    for(t=0;t<T;t++)
    {
        scanf("%d%d%d%d%d%d%d%d%d",&n,&m,&a0,&b0,&p,&q,&x,&y,&mx);
        amax=a[0]=a0%mx;
        bmax=b[0]=b0%mx;
        for(i=1;i<n;i++)
            {a[i]=(a[i-1]*p+x)%mx;if(a[i]>amax)amax=a[i];}
        for(i=1;i<m;i++)
            {b[i]=(b[i-1]*q+y)%mx;if(b[i]>bmax)bmax=b[i];}
        cc[0]=amax*bmax;
        for(i=0;i<n;i++)
        for(j=0;j<m;j++)
        {
            count[t]+=cc[a[i]*b[j]];
        }
    }
    for(i=0;i<T;i++)
        printf("%lld\n",count[i]);
    return 0;
}