hdu 6053 TrickGCD [莫比乌斯函数]
来源:互联网 发布:ubuntu 16.04改中文 编辑:程序博客网 时间:2024/06/05 09:47
题意:给出长度为n的A数列,求满足条件的B数组的个数,条件:①1<=b[i]<=a[i] ②对于任意区间【L,R】,区间gcd>=2
题解:对于某个数及其倍数能产生的B数组方案个数为:
但其中会产生重复的数列,这时候我们需要判断x的质因子数的奇偶性,若为奇数则加上该方案数,否则减去该方案数(容斥原理)
可以用莫比乌斯函数判断当前枚举的GCD对答案的贡献。
AC代码:
#include<stdio.h>#include<iostream>#include<algorithm>#include<string.h>#define mod 1000000007using namespace std;typedef long long ll;ll mu[100005],cas=1;ll num[200005];void mobius(ll mn){ mu[1]=1; for(ll i=1;i<=mn;i++){ for(ll j=i+i;j<=mn;j+=i){ mu[j]-=mu[i]; } }}ll a[100005];ll qmi(ll a,ll b){ll ans=1;while(b){if(b%2==1)ans=ans*a%mod;a=(a*a)%mod;b/=2;}return ans;}int main(){mobius(100000);ll T;scanf("%lld",&T);while(T--){memset(num,0,sizeof(num));ll n;scanf("%lld",&n);ll mi=100005;for(ll i=1;i<=n;i++){scanf("%lld",&a[i]);mi=min(a[i],mi);}for(int i=1;i<=n;i++)num[a[i]]++;for(int i=1;i<=200000;i++)num[i]+=num[i-1];ll sum=0;for(ll i=2;i<=mi;i++){ll gg=1;for(ll j=1;j*i<=100000;j++){gg=(gg*qmi(j,num[(j+1)*i-1]-num[(j)*i-1]))%mod;}sum=(sum-gg*mu[i]%mod+mod)%mod;}printf("Case #%lld: %lld\n",cas++,sum);}}
阅读全文
1 0
- hdu 6053 TrickGCD [莫比乌斯函数]
- Hdu-6053 TrickGCD(莫比乌斯函数)
- HDU 6053 TrickGCD (莫比乌斯函数)
- HDU 6053 TrickGCD(莫比乌斯函数)
- hdu 6053 TrickGCD 容斥 莫比乌斯
- HDU 6053 TrickGCD(莫比乌斯反演+分块)
- HDU 6053 TrickGCD 莫比乌斯反演||筛法
- 【HDU 6053 TrickGCD】 + 莫比乌斯反演
- hdu 6053 TrickGCD 【容斥&筛法|莫比乌斯】
- HDU 6053 TrickGCD(莫比乌斯反演)
- 【数论-莫比乌斯】hdu 6053 TrickGCD
- hdu 6053 TrickGCD(莫比乌斯反演)
- HDU 6053 TrickGCD(莫比乌斯反演+前缀和)
- HDU 6053 TrickGCD(莫比乌斯反演)
- hdu 6053 TrickGCD 筛法+莫比乌斯函数+分块处理
- hdu 6053 TrickGCD(容斥,分段,莫比乌斯函数)
- HDU 6053 TrickGCD (莫比乌斯函数+分块筛)
- HDU-6053 TrickGCD 前缀和数组,莫比乌斯函数
- JSP其他主题
- WKWebView禁止捏合的方法
- POJ 2502 Subway(将各种数据转化成图+最短路+迪杰斯特拉算法)
- (五)View的测量 MeasureSpec类
- git版本控制
- hdu 6053 TrickGCD [莫比乌斯函数]
- Centos7 udev
- 模拟购物车的实现过程(详细讲解)
- Java_Day13
- 使用Generic Webhook Trigger插件实现Jenkins+WebHooks(码云)持续集成--指定具体分支集成
- Java的流程控制
- ZOJ 3699 (贪心模拟)【加油站类模板】
- (六)VIew的绘制过程-measure、layout、draw
- Java知识(标识符、数据的存储和类型)