小沙的算术题解报告

小沙的算术题解报告

标签：逆元 数学

题目链接

来源：牛客网

解题思路：

首先从示例开始分析，
给出算式1+1+3 * 1+1 * 6 * 7；
我们可以将式子分为1,1，3 * 1,1 * 6 * 7，四部分存在num 1~4中
设置指针数标记各个数所属的分组，记此时总和为ans；
当我们蒋第六位数字6改为7时，
首先从ans中减去num[4]，然令num[4] / 6 * 7，再令 ans += num[4]即可
然后我就发现自己挖了

逆元

问题出在题目数据太大而进行了取模运算，取模和除法的顺序不能换
比如我们将题目所给的mod缩小为5
根据上面的方法算
num[4]=((16%5)7)%5=2
变化后num[4]=2 / 6 * 7 % 5 =
2/6除不尽问题就暴露出来了 为防止出现这种情况我们要使用逆元
本人也是初学逆元，不是很懂，讲不明白
主要参考学习了下面这篇文章
https://zhuanlan.zhihu.com/p/378728642
根据费马小定理 当a与b互质时：

运用到上述例子中求出6的逆元为216；
num[4]=2216%57%5=4;
num[4]=7*7%5=4；
这就对了
我们可以大胆猜测1e9+7为质数，通过快速幂得到a模b的逆元

代码实现：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod = 1e9+7;
const int N = 1e6+8;
char s[N];
ll num[N],pos[N],n,q,a[N];
ll ksm(ll a, ll b) {
    ll ans = 1;
    while(b){
		if (b & 1) (ans *= a) %= mod;
		(a *= a) %= mod; 
		b >>= 1;
	}
    return ans;
}
int main()
{
    cin>>n>>q;
    scanf("%s",s+1);
    for(int i=1;i<=n;i++) cin>>a[i];
    int cnt=1;
    num[1]=a[1];
    pos[1]=1;
    for(int i=2;i<=n;i++){
        if(s[i-1]=='*'){
            pos[i]=cnt;
            num[cnt]=num[cnt]*a[i]%mod;
        }
        else{
            num[++cnt]=a[i];
            pos[i]=cnt;
        }
    }
    ll ans = 0,x=0,y=0;
    for(int i=1;i<=cnt;i++)  (ans+=num[i]) %= mod;
    while(q--){
        cin>>x>>y;
        ans=(ans-num[pos[x]]+mod)%mod;
        num[pos[x]]=num[pos[x]]*ksm(a[x],mod-2)%mod*y%mod;
        ans=(ans+num[pos[x]])%mod;
        a[x]=y;
        cout<<ans<<endl;
    }
}