zoukankan      html  css  js  c++  java
  • BZOJ 3930: [CQOI2015]选数

    3930: [CQOI2015]选数

    Time Limit: 10 Sec  Memory Limit: 512 MB

    Description

     我们知道,从区间[L,H](L和H为整数)中选取N个整数,总共有(H-L+1)^N种方案。小z很好奇这样选出的数的最大公约数的规律,他决定对每种方案选出的N个整数都求一次最大公约数,以便进一步研究。然而他很快发现工作量太大了,于是向你寻求帮助。你的任务很简单,小z会告诉你一个整数K,你需要回答他最大公约数刚好为K的选取方案有多少个。由于方案数较大,你只需要输出其除以1000000007的余数即可。

     

    Input

    输入一行,包含4个空格分开的正整数,依次为N,K,L和H。

     

    Output

    输出一个整数,为所求方案数。

     

    Sample Input

    2 2 2 4

    Sample Output

    3

    HINT

     样例解释


    所有可能的选择方案:(2, 2), (2, 3), (2, 4), (3, 2), (3, 3), (3, 4), (4, 2), (4, 3), (4, 4)

    其中最大公约数等于2的只有3组:(2, 2), (2, 4), (4, 2)

    对于100%的数据,1≤N,K≤10^9,1≤L≤H≤10^9,H-L≤10^5

     题解大法好:

    http://blog.csdn.net/L_0_Forever_LF/article/details/53031515

     1 #include<iostream>
     2 #include<algorithm>
     3 #include<cstdio>
     4 #include<cstring>
     5 #include<cmath>
     6 #include<cstdlib>
     7 #include<vector>
     8 using namespace std;
     9 typedef long long ll;
    10 typedef long double ld;
    11 typedef pair<int,int> pr;
    12 const double pi=acos(-1);
    13 #define rep(i,a,n) for(int i=a;i<=n;i++)
    14 #define per(i,n,a) for(int i=n;i>=a;i--)
    15 #define Rep(i,u) for(int i=head[u];i;i=Next[i])
    16 #define clr(a) memset(a,0,sizeof(a))
    17 #define pb push_back
    18 #define mp make_pair
    19 #define fi first
    20 #define sc second
    21 #define pq priority_queue
    22 #define pqb priority_queue <int, vector<int>, less<int> >
    23 #define pqs priority_queue <int, vector<int>, greater<int> >
    24 #define vec vector
    25 ld eps=1e-9;
    26 ll pp=1000000007;
    27 ll mo(ll a,ll pp){if(a>=0 && a<pp)return a;a%=pp;if(a<0)a+=pp;return a;}
    28 ll powmod(ll a,ll b,ll pp){ll ans=1;for(;b;b>>=1,a=mo(a*a,pp))if(b&1)ans=mo(ans*a,pp);return ans;}
    29 void fre() { freopen("c://test//input.in", "r", stdin); freopen("c://test//output.out", "w", stdout); }
    30 //void add(int x,int y,int z){ v[++e]=y; next[e]=head[x]; head[x]=e; cost[e]=z; }
    31 int dx[5]={0,-1,1,0,0},dy[5]={0,0,0,-1,1};
    32 ll read(){ ll ans=0; char last=' ',ch=getchar();
    33 while(ch<'0' || ch>'9')last=ch,ch=getchar();
    34 while(ch>='0' && ch<='9')ans=ans*10+ch-'0',ch=getchar();
    35 if(last=='-')ans=-ans; return ans;
    36 }
    37 const int N=1e05+5;
    38 ll f[N];
    39 int main(){
    40     int n=read(),k=read(),l=read(),h=read();
    41     int l_=l/k+((l%k)!=0),r_=h/k;
    42     for (int i=h-l+1;i>=1;i--){
    43         int L=l_/i+((l_%i)!=0),R=r_/i;
    44         if (L<=R){
    45             f[i]=mo(powmod(R-L+1,n,pp)-(R-L+1),pp);
    46                     for(int j=2*i;j<=N;j+=i) f[i]=mo(f[i]-f[j],pp);
    47         }
    48     }
    49     if(l_<=1&&l_<=r_) f[1]=mo(f[1]+1,pp);
    50         printf("%lld
    ",f[1]);
    51     return 0;
    52 }
    View Code
  • 相关阅读:
    专职DBA-MySQL体系结构与基本管理
    JSON
    MIME类型
    文件上传下载
    response常用的方法
    2020.11.27小记
    HTTP请求状态码
    1561. Maximum Number of Coins You Can Get
    1558. Minimum Numbers of Function Calls to Make Target Array
    1557. Minimum Number of Vertices to Reach All Nodes
  • 原文地址:https://www.cnblogs.com/SXia/p/7701363.html
Copyright © 2011-2022 走看看