Modified LCS

题目链接

• 题意：
  给两个等差数列的长度，起点和数列的增加值，求两个数列中有几个数相同
• 分析：
  将等差数列的通项公式化简后可以得到扩展欧几里得的结构，直接计算即可
• 反思：
  求出方程的一个解后，得到的是下标。此时如果继续用下标判断会比较麻烦，因为同一个数在两个序列中的下标是不一样的，所以需要两个数列的下标范围均需要判断是否合法。而如果采用值判断，因为两个数列的满足题意的值是相同的，所以直接采用值判断会简单很多
  
  但是，对于这个题目，上边的方法会带来一个问题，题目中的最后一项有可能超过long long的范围从而出错。正确的方法是，求出第一个数后转换成下标，x为在第一个数列中的下标，y为在第二个数列中的下标，dx表示x的变化量，dy表示y的变化量。那么ans = min((n1 - x) / dx, (n2 - y) / dy) + 1，这样转换成下标操作就可以避免计算最后一个数是多少（出错）

void ex_gcd(LL a, LL b, LL &x, LL &y, LL &d){    if (!b)        x = 1, y = 0, d = a;    else    {        ex_gcd(b, a % b, y, x, d);        y -= x * (a / b);    }}int main(){    //freopen("0.txt", "r", stdin);    int T;    RI(T);    FE(kase, 1, T)    {        LL n1, n2, f1, f2, d1, d2;        LL a, b, v, x, y, dx, dy, gcd, lcm;        cin >> n1 >> f1 >> d1 >> n2 >> f2 >> d2;        a = d1, b = -d2, v = d1 - d2 + f2 - f1;        ex_gcd(a, b, x, y, gcd);        lcm = abs(d1 * d2 / gcd);        if (v % gcd != 0)            puts("0");        else        {            x *= v / gcd;            x = f1 + (x - 1) * d1;            dx = abs(b / gcd * d1);            dy = abs(a / gcd);            LL l = max(f1, f2);            x = ((x - l) % dx + dx) % dx + l;            y = (x - f2) / d2 + 1;            x = (x - f1) / d1 + 1;            dx /= d1;            y = min((n1 - x) / dx, (n2 - y) / dy) + 1;            if (y >= 0)                cout << y << endl;            else                puts("0");        }    }    return 0;}