BZOJ1902: Zju2116 Christopher

$n \leq 10^{100}$，问$C_n^m,0<=m<=n$有多少是质数$p \leq 1e7$的倍数。

一样，套高精度的题，只有战胜他才能鄙视他。

但是我TM被他鄙视了一上午！！！

好先冷静分析。用Lucas的观点看组合数，这里就是个明显的数位DP了，统计每一位时大于当前数、小于等于当前数的合法和不合法方案数，很简单的转移，详见代码。

被鄙视*1：方程抄错了。。

被鄙视*2：高精度乘单精度乘法写错了。。

当然这也不能怪我鬼知道他有乘零！

好吧怪我

 //#include<iostream>
 #include<cstring>
 #include<cstdlib>
 #include<cstdio>
 //#include<vector>
 //#include<queue>
 //#include<time.h>
 //#include<complex>
 #include<algorithm>
 #include<stdlib.h>
 using namespace std;

 int n;
 #define maxn 1011
 int a[maxn],mod,b[maxn],lb; char s[maxn]; int len;

 #define LL long long
 struct LLL
 {
     int a[],len;
     LLL() {memset(a,,sizeof(a)); len=;}
     void operator = (int x)
     {
         len=;
         while (x) a[++len]=x%,x/=;
     }
     void operator = (const LLL &b)
     {
         len=b.len;
         for (int i=;i<=len;i++) a[i]=b.a[i];
     }
     LLL operator * (int x)
     {
         LLL ans;
         for (int i=;i<=len;i++)
         {
             LL tmp=a[i]*1ll*x;
             ans.a[i+]+=(ans.a[i]+tmp)/;
             ans.a[i]=(ans.a[i]+tmp)%;
         }
         ans.len=len;
         while (ans.a[ans.len+])
         {
             ans.len++;
             if (ans.a[ans.len]>=) ans.a[ans.len+]+=ans.a[ans.len]/,ans.a[ans.len]%=;
         }
         while (ans.a[ans.len]== && ans.len>) ans.len--;
         return ans;
     }
     LLL operator + (const LLL &b)
     {
         LLL ans;
         for (int i=,to=max(len,b.len);i<=to;i++)
         {
             ans.a[i]+=a[i]+b.a[i];
             if (ans.a[i]>=)
             {
                 ans.a[i+]++;
                 ans.a[i]-=;
             }
         }
         ans.len=max(len,b.len);
         while (ans.a[ans.len+]) ans.len++;
         return ans;
     }
     void out()
     {
         printf("%d",a[len]);
         for (int i=len-;i>;i--)
         {
             for (int j=;j>;j/=) if (a[i]<j) putchar('');
             printf("%d",a[i]);
         }
     }
 }f[maxn][],g[maxn];

 int main()
 {
     scanf("%s%d",s+,&mod); len=strlen(s+);
     for (int i=;i<=len;i++) a[i]=s[len-i+]-'';
     lb=; while (len)
     {
         int tmp=;
         for (int i=len;i;i--) {int now=a[i]; a[i]=(tmp*+now)/mod; tmp=(tmp*+now)%mod;}
         for (;len && a[len]==;len--);
         b[++lb]=tmp;
     }
     f[][]=; f[][]=b[]+; g[]=mod--b[];
     for (int i=;i<=lb;i++)
     {
         f[i][]=f[i-][]*(b[i]+)+g[i-]*b[i];
         f[i][]=f[i-][]*(b[i]+);
         g[i]=(f[i-][]+f[i-][])*(mod--b[i])+g[i-]*(mod-b[i]);
     }
     f[lb][].out();
     return ;
 }