AC日记——【清华集训2014】奇数国 uoj 38

#38. 【清华集训2014】奇数国

思路:

  题目中的number与product不想冲;

  即为number与product互素;

  所以,求phi(product)即可;

  除一个数等同于在模的意义下乘以一个数的逆元;

代码:

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm> using namespace std; #define maxn 100005
#define mod 19961993
#define ll long long struct TreeNodeType {
ll l,r,mid; ll dis1,dis2;
};
struct TreeNodeType tree[maxn<<]; ll n,m,cntp; ll bit[],ans1,ans2,pi_[],pi[]; bool if_p[]; inline void in(ll &now)
{
char Cget=getchar();now=;
while(Cget>''||Cget<'') Cget=getchar();
while(Cget>=''&&Cget<='')
{
now=now*+Cget-'';
Cget=getchar();
}
} ll poww(ll pos)
{
pos%=mod;
ll mi=mod-,res=;
while(mi)
{
if(mi&) res=(res*pos)%mod;
mi>>=,pos=(pos*pos)%mod;
}
return res;
} void tree_build(ll now,ll l,ll r)
{
tree[now].l=l,tree[now].r=r;
if(l==r)
{
tree[now].dis1=bit[],tree[now].dis2=;
return ;
}
tree[now].mid=l+r>>;
tree_build(now<<,l,tree[now].mid);
tree_build(now<<|,tree[now].mid+,r);
tree[now].dis1=tree[now<<].dis1|tree[now<<|].dis1;
tree[now].dis2=tree[now<<].dis2*tree[now<<|].dis2%mod;
} void tree_to(ll now,ll to,ll dis1,ll dis2)
{
if(tree[now].l==tree[now].r)
{
tree[now].dis1=dis1,tree[now].dis2=dis2;
return ;
}
if(to<=tree[now].mid) tree_to(now<<,to,dis1,dis2);
else tree_to(now<<|,to,dis1,dis2);
tree[now].dis1=tree[now<<].dis1|tree[now<<|].dis1;
tree[now].dis2=tree[now<<].dis2*tree[now<<|].dis2%mod;
} void tree_query(ll now,ll l,ll r)
{
if(tree[now].l==l&&tree[now].r==r)
{
ans1|=tree[now].dis1;
ans2=(ans2*tree[now].dis2)%mod;
return ;
}
if(l>tree[now].mid) tree_query(now<<|,l,r);
else if(r<=tree[now].mid) tree_query(now<<,l,r);
else tree_query(now<<,l,tree[now].mid),tree_query(now<<|,tree[now].mid+,r);
} int main()
{
for(ll i=;i<=;i++)
{
if(!if_p[i]) pi[++cntp]=i;
for(ll j=;pi[j]*i<=&&j<=cntp;j++)
{
if_p[i*pi[j]]=true;
if(i%pi[j]==) break;
}
}
for(ll i=;i<=;i++)
{
pi_[i]=poww(pi[i]);
if(i==) bit[i]=;
else bit[i]=bit[i-]<<;
}
in(n);tree_build(,,maxn-);
ll op,ai,bi;
for(;n--;)
{
in(op),in(ai),in(bi);
if(op==)
{
ll pos=;
for(ll i=;i<=;i++) if(bi&&(bi%pi[i]==)) pos+=bit[i];
bi%=mod;tree_to(,ai,pos,bi);
}
else
{
ll ans;ans1=,ans2=,tree_query(,ai,bi),ans=ans2;
for(ll i=;i<=;i++) if(bit[i]&ans1) ans=(ans*((pi[i]-)*pi_[i]%mod))%mod;
printf("%lld\n",ans);
}
}
return ;
}
上一篇:java.util (Collection接口和Map接口)


下一篇:清华集训2014 day2 task3 矩阵变换