Gym102040 .Asia Dhaka Regional Contest(寒假自训第9场)

B .Counting Inversion

题意：给定L,R，求这个区间的逆序对数之和。(L,R<1e15)

思路：一看这个范围就知道是数位DP。 只是维护的东西稍微多一点，需要记录后面的各种数字的个数cnt，以及逆序对和sum，以及出现了多少种后缀num。

那么枚举到当前位时，假设为i ，那么sum+=cnt[i+1]+cnt[i+2]+....cnt[9];  cnt[i]+=num; 可以参考CF1073E。

#include<bits/stdc++.h>
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define ll long long
using namespace std;
struct in{
    ll num,cnt[],sum;
    in(){num=sum=; memset(cnt,,sizeof(cnt)); }
}dp[];
int q[],tot,vis[];
in dfs(int pos,int st,int lim)
{
    if(!lim&&vis[pos]) return dp[pos];
    if(pos==) {
        in res;  res.num=;
        return res;
    }
    int up=; in res,tmp; if(lim) up=q[pos-];
    rep(i,,up){
        tmp=dfs(pos-,i,lim&&i==up);
        res.sum+=tmp.sum;
        rep(j,i+,) res.sum+=tmp.cnt[j];
        rep(j,,) res.cnt[j]+=tmp.cnt[j];
        res.cnt[i]+=tmp.num;
        res.num+=tmp.num;
    }
    vis[pos]=;
    return dp[pos]=res;
}
ll cal(ll x)
{
    if(x<) return ;
    tot=; ll ans=;
    while(x) q[++tot]=x%,x/=;
    memset(dp,,sizeof(dp));
    memset(vis,,sizeof(vis));
    rep(i,,tot){
        ll up=; if(i==tot) up=q[tot];
        rep(j,,up){
             in tmp=dfs(i,j,(i==tot)&&(j==q[tot]));
             ans+=tmp.sum;
             rep(k,j+,) ans+=1LL*tmp.cnt[k];
        }
    }
    return ans;
}
int main()
{
    ll L,R;  int T,Ca=; scanf("%d",&T);
    while(T--){
       scanf("%lld%lld",&L,&R);
       printf("Case %d: %lld\n",++Ca,cal(R)-cal(L-));
    }
    return ;
}

C .Divisors of the Divisors of An Integer

题意：给出N，问N！的因子的因子个数和。

思路：唯一分解，对于一个素数p，假设它的幂次是x，那么因子的幂次有0，1，2，...x；那么因子的因子幂次就是(0)； (0，1)； ( 0，1，2)；   ... ; (0，1，2，...x)

所以就是一个累乘，对于每个素数p，ans*=(x+1)*(x+2)/2;    而阶乘的唯一分解只需要一直除即可。

#include<bits/stdc++.h>
#define ll long long
#define rep(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
const int maxn=;
const int Mod=1e7+;
int a[maxn],N,p[maxn],vis[maxn],cnt;ll ans=;
int solve(int p)
{
    int tN=N,res=;
    while(tN) {
        res+=tN/p;
        if(res>Mod) res-=Mod;
        tN/=p;
    } return res;
}
int get(int p)
{
    if(p&) return 1LL*(p+)/*p%Mod;
    return 1LL*p/*(p+)%Mod;
}
int main()
{
    rep(i,,){
        if(!vis[i]) p[++cnt]=i;
        for(int j=;j<=cnt&&p[j]*i<=;j++){
            vis[p[j]*i]=;
            if(i%p[j]==) break;
        }
    }
    while(~scanf("%d",&N)&N){
       ans=;
       rep(i,,N) {
          if(!vis[i]) a[i]=solve(i);
       }
       rep(i,,N){
         if(a[i])
           ans=(ll)ans*get(a[i]+)%Mod;
       }
       printf("%lld\n",ans);
    }
    return ;
}

E.Helping the HR

题意：给定每个人的签到和离开时间，问每个人的..情况

思路：模拟； by许。

#include<bits/stdc++.h>
using namespace std;
char str[];
int main()
{
    int n,s=*,D=*,E=*;
    while(~scanf("%d",&n)&&n)
    {
        int cnt=;
        for(int cas=;cas<n;cas++)
        {
            int S=,T=,tot=,pre=,p=;
            scanf("%s",str);
            int len=strlen(str);
            for(int i=;i<len;i++)
            {
                if(str[i]!=':'&&i!=len-)pre=pre*+str[i]-'';
                else
                {
                    if(i==len-)pre=pre*+str[i]-'';
                    tot++;
                    if(tot<=)
                    {
                        S+=pre*p;
                        if(tot==)p=;
                        else p/=;
                    }
                    else T+=pre*p,p/=;
                    pre=;
                }
            }
            int flag=;
            if(str[]=='D'&&S>D)flag=;
            if(str[]=='E'&&S>E)flag=;
            int res=T-max(s,S);
            if(str[]=='D'&&res<*)flag=;
            if(str[]=='E'&&res<*)flag=;
            cnt+=flag;
        }
        if(!cnt)puts("All OK");
        else if(cnt<=)printf("%d Point(s) Deducted\n",cnt);
        else puts("Issue Show Cause Letter");
    }
}

F .Path Intersection

题意：给定一棵树, Q次询问，每次给定K条路经，求这K条路有多少个公共点.

思路：路剖,  即每条路经+1， 然后可以选一条路径看有多少个点被覆盖K次。

好久没写树剖了，开始写错的地方提醒下自己：  当top不同的时候，我们先操作dep[top[]]大的，然后把它变为fa[top[]];

而top相同的，正常的从小到大即可。

#include<bits/stdc++.h>
#define rep(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
const int maxn=;
int Laxt[maxn],Next[maxn],To[maxn],cnt,dep[maxn];
int sz[maxn],son[maxn],top[maxn],pos[maxn],N;
int Mx[maxn],num[maxn],Lazy[maxn],fa[maxn],tot;
void add(int u,int v){
    Next[++cnt]=Laxt[u]; Laxt[u]=cnt; To[cnt]=v;
}
void dfs1(int u,int f)
{
    sz[u]=;  fa[u]=f;
    dep[u]=dep[f]+; son[u]=;
    for(int i=Laxt[u];i;i=Next[i]){
        if(To[i]!=f){
            dfs1(To[i],u);
            if(sz[To[i]]>sz[son[u]]) son[u]=To[i];
        }
    }
}
void dfs2(int u,int Top)
{
    pos[u]=++tot; top[u]=Top;
    if(son[u]) dfs2(son[u],Top);
    for(int i=Laxt[u];i;i=Next[i]){
      if(To[i]!=fa[u]&&To[i]!=son[u])
         dfs2(To[i],To[i]);
    }
}
void build(int Now,int L,int R)
{
    Mx[Now]=Lazy[Now]=; num[Now]=R-L+;
    if(L==R) return; int Mid=(L+R)>>;
    build(Now<<,L,Mid); build(Now<<|,Mid+,R);
}
void pushdown(int Now)
{
    if(Lazy[Now]) {
        Mx[Now<<]+=Lazy[Now];Lazy[Now<<]+=Lazy[Now];
        Mx[Now<<|]+=Lazy[Now];Lazy[Now<<|]+=Lazy[Now];
        Lazy[Now]=;
    }
}
void pushup(int Now)
{
    Mx[Now]=Mx[Now<<]; num[Now]=num[Now<<];
    if(Mx[Now<<|]>Mx[Now])
      Mx[Now]=Mx[Now<<|],num[Now]=num[Now<<|];
    else if(Mx[Now<<|]==Mx[Now])
      num[Now]+=num[Now<<|];
}
void update(int Now,int L,int R,int l,int r,int val)
{
    if(l<=L&&r>=R){
        Mx[Now]+=val; Lazy[Now]+=val; return ;
    }
    pushdown(Now); int Mid=(L+R)>>;
    if(l<=Mid) update(Now<<,L,Mid,l,r,val);
    if(r>Mid) update(Now<<|,Mid+,R,l,r,val);
    pushup(Now);
}
int query(int Now,int L,int R,int l,int r,int K)
{
    if(Mx[Now]<K) return ;
    if(l<=L&&r>=R) return Mx[Now]==K?num[Now]:;
    pushdown(Now); int Mid=(L+R)>>,res=;
    if(l<=Mid) res+=query(Now<<,L,Mid,l,r,K);
    if(r>Mid) res+=query(Now<<|,Mid+,R,l,r,K);
    pushup(Now); return res;
}
void pathup(int u,int v,int val)
{
    while(top[u]!=top[v]){
        if(dep[top[u]]<dep[top[v]]) swap(u,v);
        update(,,N,pos[top[u]],pos[u],val);
        u=fa[top[u]];
    }
    if(dep[u]>dep[v]) swap(u,v);
    update(,,N,pos[u],pos[v],val);
}
int pathquery(int u,int v,int K)
{
    int res=;
    while(top[u]!=top[v]){
        if(dep[top[u]]<dep[top[v]]) swap(u,v);
        res+=query(,,N,pos[top[u]],pos[u],K);
        u=fa[top[u]];
    }
    if(dep[u]>dep[v]) swap(u,v);
    res+=query(,,N,pos[u],pos[v],K);
    return res;
}
int a[maxn],b[maxn];
int main()
{
    int T,Q,K,C=,u,v;
    scanf("%d",&T);
    while(T--){
        scanf("%d",&N);
        rep(i,,N) Laxt[i]=; cnt=; tot=;
        rep(i,,N-){
            scanf("%d%d",&u,&v);
            add(u,v); add(v,u);
        }
        dfs1(,); dfs2(,);
        build(,,N);
        scanf("%d",&Q);
        printf("Case %d:\n",++C);
        while(Q--){
            scanf("%d",&K);
            rep(i,,K) scanf("%d%d",&a[i],&b[i]);
            rep(i,,K) pathup(a[i],b[i],);
            printf("%d\n",pathquery(a[],b[],K));
            rep(i,,K) pathup(a[i],b[i],-);
        }
    }
    return ;
}

I .Triangles

题意：给定两个三维空间里的三角形，求最近距离。

思路：好像是不错的题，想补。

J. VAT Man

签到。 by许。

#include<bits/stdc++.h>
#define db double
using namespace std;
int main()
{
    int T;
    cin>>T;
    while(T--)
    {
        db x;
        cin>>x;
        printf("%.2lf\n",x*1.15);
    }
}