HDU 5768 Lucky7 Theory of Chinese remainder theorem

The original topic see HDU 5768

In the [L,r] range is a multiple of 7 and does not satisfy the number of any given congruence. If the range is [1,100], does not meet the modulus of 3 to 2 or more than 3 of the modulus of 5 is 7,21,42,49,70,84,91, so the answer is 7.
Where the divisor is a prime number other than 7 (≤105) (\le 10^5), the product of the divisor is less than 1018 10^{18}. The same redundancy has a maximum of 15. Analysis

A multiple of 7 in all [L,r] range minus the number that satisfies any one congruence, that is, the answer is obtained. And when a number can satisfy multiple congruence, it is natural to think of the principle of tolerance to solve the problem. This collection is subtracted when the odd number of conditions are met, and vice versa.

The Chinese remainder theorem solves the solution that satisfies multiple congruence equations at the same time. Therefore 215 2^{15} within the scope of a variety of combinations of the same redundancy, you can get the answer.

Solution equation set X=ai (mod mi) mi between 22 coprime
int China (int r) {
    int M = 1, ans = 0;
    for (int i = 0; i < R; ++i)
        M *= m[i];
    for (int i = 0;i < r;i++) {
        int N = m/m[i];
        int x, y;
        Extend_euclid (N, M[i], x, y);
        x = (X%m[i]+m[i])%m[i];
        Ans = ((ans+a[i]*n%m*x%m)%M + M)%M;
    }
    return ans;
}

Pay special attention to the middle solution of each satisfied x, to take the smallest positive integer solution, or the subsequent multiplication will explode a long long (even if not explode, the answer is not, after all, M[i] is the approximate number of M, X after processing will become smaller). Ans is a constant iterative solution that is to be modeled M.
The person is still good, in the process of the ANS can also make a data explosion long long, so it is fast multiplication to engage. Note that the multiplier must be a positive number (small bin died in this hhh) code

/*--------------------------------------------* File name:hdu 5768 * Author:danliwoo * mail:danliwoo@outlook.com * Created time:2016-07-28 18:59:37--------------------------------------------*/#include <bits/stdc++.h> using
namespace Std;
#define N #define LL Long long LL a[n], m[n];
int s[n], N;
        ll Extend_euclid (ll A, ll B, LL &x, ll &y) {if (b==0) {x = 1; y = 0;
    return A;
    } LL r = extend_euclid (b, A%b, y, x);
    Y-= a/b*x;
return R;
} ll Gao (ll X, LL R, ll p) {return (x-r)/p;}
    ll Mult (ll A, LL K, ll m) {ll res = 0;
        while (k) {if (K & 1LL) res = (res + a)% m;
        K >>= 1;
    A = (a << 1)% m;
} return res;
    } ll China (ll L, ll R) {ll M = 1, ans = 0;
    for (int i = 0; I <= n; ++i) if (S[i]) {M *= m[i];
        } for (int i = 0;i <= n;i++) if (S[i]) {LL Nn = m/m[i];
        LL x, y;
        Extend_euclid (Nn, M[i], x, y); x = (X%m[i] + m[i])% m[i];
    Ans = ((Ans+mult (a[i]*nn%m, X, M))%M + m)% m;
    } LL ret = Gao (r+m, ans, m)-Gao (l-1+m, ans, M);
return ret;
    } int main () {int T, o = 0;
    scanf ("%d", &t);
        while (t--) {LL L, R;
        scanf ("%d%lld%lld", &n, &l, &r);
        memset (s, 0, sizeof (s)); M[n] = 7; A[n] = 0;
        S[n] = 1;
        for (int i = 0;i < n;i++) scanf ("%lld%lld", &m[i], &a[i]);
        LL ans = 0;
        int all = 1 << n;
            for (int i = 0;i < all;i++) {int t = i, k = 0;
                for (int j = 0;j < n;j++) {S[j] = t & 1;
                T >>= 1;
            K + = S[j]; } k = k & 1?
            -1:1; 
        Ans + = 1LL * k * China (L, R);
    } printf ("Case #%d:%lld\n", ++o, ans);
} return 0;
 }