[bzoj1416]神奇的口袋

容易发现操作任意次并不会改变每一个点的概率（因为每一个点加d的概率相同，期望与原数成比例），然后直接输出即可（要用高精度）

#include<bits/stdc++.h>
using namespace std;
#define N 20005
int n,m,k,x,s,a[N],vis[N],p[N],mi[N],sum[N];
void calc(){
    memset(a,0,sizeof(a));
    a[0]=a[1]=1;
    for(int i=1;i<N-4;i++)
        for(int j=1;j<=sum[i];j++){
            a[1]=a[1]*i;
            for(int k=2;k<=a[0];k++){
                a[k]=a[k]*i+a[k-1]/10;
                a[k-1]%=10;
            }
            while(a[a[0]]>9){
                a[a[0]+1]=a[a[0]]/10;
                a[a[0]++]%=10;
            }
        }
    for(int i=a[0];i;i--)printf("%d",a[i]);
}
int main(){
    for(int i=2;i<N-4;i++){
        if (!vis[i])p[++p[0]]=i;
        for(int j=1;(j<=p[0])&&(i*p[j]<N-4);j++){
            vis[i*p[j]]=1;
            mi[i*p[j]]=p[j];
            if (i%p[j]==0)break;
        }
    }
    scanf("%d%d%d",&n,&m,&k);
    for(int i=1;i<=n;i++){
        scanf("%d",&a[i]);
        s+=a[i];
    }
    for(int i=1;i<=m;i++){
        scanf("%*d%d",&x);
        sum[a[x]]++;
        sum[s]--;
        a[x]+=k;
        s+=k;
    }
    for(int i=N-5;i;i--){
        if (vis[i]){
            sum[mi[i]]+=sum[i];
            sum[i/mi[i]]+=sum[i];
            sum[i]=0;
        }
    }
    calc();
    printf("/");
    for(int i=N-5;i;i--)sum[i]*=-1;
    calc();
}