codeforces 181 div2 C&E

比赛没打，B题估计写不熟练，C题虽然赛后搞定了，但不是很熟练，比赛做坑定不行

这题正面思考会很吐，尤其是YYdp的话，就及其不可取了。从一个条件入手，就是数字之和为good数，由于n=10^6,所以good数，也就是所有excellent数的数字之后也不过就是2^1+2^2+2^3+..2^7<=2^8，才200多种，很少。枚举结果，由于数字不是a就是b,再加上题目限定条件,必须是n位，好了，方程组ax+by=sum(digit),x+y=n，x和y分别代表各种数字的个数。接下来求组合数C(n,x)%mod,这个要用到逆元。。。。其实这题不过就是个意识的问题，意识快，白送之。。。附代码：

#include <iostream>#include <cstdio>#define u 127#define mod 1000000007  #define ll __int64using namespace std;ll sum;void extend(ll a0,ll b0,ll &x,ll &y)  {      if (!b0)      {          x=1;y=0;      }      else      {          extend(b0,a0%b0,x,y);          ll temp=x;x=y;          y=temp-a0/b0*y;      }  }  int com(int n,int k){    if (k>n-k) k=n-k;    ll p=1,q=1,x,y;    for (int i=1;i<=k;i++)    {        p=p*(n+1-i)%mod;        q=(q*i)%mod;    }    extend(q,mod,x,y);    x=(x+mod)%mod;    ll res=(p*x)%mod;    return (int)res;}void dfs(int value,int a,int b,int n){    if (value>b*n) return;    int s=value-a*n;    if (s>=0&&s%(b-a)==0) sum=(sum+com(n,s/(b-a)))%mod;    dfs(value*10+a,a,b,n);    dfs(value*10+b,a,b,n);}int main(){    int n,a,b;    cin>>a>>b>>n;    sum=0;    dfs(0,a,b,n);    cout<<sum<<endl;    return 0;}

E.请教cxlove爱神的，主要是如何处理分母的各个质因子。普通的作法坑定超时，于是就把a数组中所有的绑定在一起，不一一计算a[i]!的质因子情况，用cnt[i]表示数组中有多少个数大于等于i。然后枚举每个10^7的质因数，按照a[i],a[i]^2,a[i]^m..逐次筛选。预处理完了之后，再二分。最小值是max(a[i]),最大值是sum(a[i]),特别提醒，再用max(a[i])时，最好不要直接写max(x,y)函数，这个函数不支持int64，附代码：

#include <iostream>#include <cstdio>#include <vector>#define N 10000006#define ll __int64#define pb push_backusing namespace std;ll cnt[N],a[N];vector<ll>b,c;int check(ll n){    ll i,j,k,z;    for (i=0;i<b.size();i++)    {        j=b[i];        z=0;        k=j;        while (true)        {            z+=n/k;            if (n/k>=j) k*=j;            else break;        }        if (z<c[i]) return 0;     }    return 1;}ll search(ll l,ll r){    if (r-l<=1)    {        if (check(l)) return l;        else return r;    }    ll m =(l+r)/2;    if (check(m)) return search(l,m);    else return search(m+1,r);}int main(){    ll n,x,t,i,j,k,mmax=-1,zs=0;    cin>>n;    for (i=1;i<=n;i++)     {        cin>>x;        cnt[x]++;        if (x>mmax) mmax = x;        zs+=x;    }    for (i=N-2;i>=1;i--) cnt[i]+=cnt[i+1];    for (i=2;i<=mmax;i++)        if (!a[i])        {            t=0;            j=i;            while (j<=mmax)            {                               for (k=j;k<=mmax;k+=j)                 {                    if (k!=i) a[k]=1;                    t+=cnt[k];                }                if (mmax/j>=i) j*=i;                else j=mmax+1;            }            b.pb(i);            c.pb(t);        }    cout<<search(mmax,zs)<<endl;    //system("pause");    return 0;}