[树状数组]uva 5798 Jupiter Atacks!

H(f_i,..., f_j) = B^kf_j-k    ( mod P)

/**uva 5798 
动态更新fi 的值，查询H(fi,...fj)写出f11,f12,...f15;f22,f23,,,f25;....
以最长的f15 作参照, 设 a1...a5 分别为f15的第5....1项，sum 为a的和， 
则 f(i,j) = (sum(j) - sum(i-1) ) / b^(n-j) % p
求和用树状数组，模p用欧拉定理(题目关键2≤B <P≤109 and P prime)。trick：注意减法会产生负数
*/
#include <stdio.h>
#include <string.h>

#define N 100001
#define lowbit(i) (i) & (-i)
__int64 B[N],f[N],A[N];
int n;
void add(int i,__int64 c)
{
    for(; i <= n;i += lowbit(i