HNOI-排队

ans=n!*A(n+1,2)*A(n+3,m)+n!*C(n+1,1)*2*m*A(n+2,m-1);
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<iostream>
using namespace std;

int n,m;

struct bign
{
    long long a[10010];
    int len;  
};

bign nn,aa;

void init()
{
    freopen("wiki.in","r",stdin);
    //freopen("wiki.out","w",stdout);
}

void readdata()
{
    scanf("%d%d",&n,&m);
}

void wocao()
{
    long long a[10010];memset(a,0,sizeof(a));
    a[1]=1;
    int len=1;
    for(int i=2;i<=n;i++)
    {
        long long c=0;
        for(int j=1;j<=len;j++)
        {
            long long s=a[j]*i+c;
            a[j]=s%100000;
            c=s/100000;
        }
        while(c > 0) 
        {
            len++;
            a[len]=c%100000;
            c/=100000;
        }
    }
    memcpy(&nn.a,&a,sizeof(a));
    nn.len=len;
    memcpy(&aa.a,&a,sizeof(a));
    aa.len=len;
}

void wogan(int u,int v)
{
    long long a[10010],c;memset(a,0,sizeof(a));
    a[1]=1;
    int len=1;
    for(int i=u;i>=u-v+1;i--)
    {
        c=0;
        for(int j=1;j<=len;j++)
        {
            long long s=a[j]*i+c;
            a[j]=s%100000;
            c=s/100000;
        }
        while(c > 0) 
        {
            len++;
            a[len]=c%100000;
            c/=100000;
        }
    }
    //*********************************
    long long b[10010],f[10010];
    int lenb,lenf=0;
    memset(f,0,sizeof(f));
    for(int i=1;i<=nn.len;i++)
    {
        memset(b,0,sizeof(b));
        long long x=nn.a[i],ss=0;int g=i;
        for(int j=1;j<=len;j++)
        {
            long long y=a[j]*x+ss;
            if(j==len)
            {
                b[g]=y%100000;
                b[g+1]=y/100000;
                if(y/100000!=0) lenb=g+1;
                else lenb=g;
            }
            else 
            {
                b[g++]=y%100000;
                ss=y/100000;
            }
        }
        lenf=max(lenf,lenb);
        for(int i=1;i<=lenf;i++)
        {
            f[i]=f[i]+b[i];
        }
        int ee;
        for(int i=1;i<=lenf;i++)
        {
            ee=f[i]/100000;
            f[i]=f[i]%100000;
            f[i+1]+=ee;
        }
        if(f[lenf+1]!=0) lenf++;
    }
    //*********************************
    memcpy(&nn.a,&f,sizeof(f));
    nn.len=lenf;
}

void wogan1(int u,int v)
{
    long long a[10010],c;memset(a,0,sizeof(a));
    a[1]=1;
    int len=1;
    for(int i=u;i>=u-v+1;i--)
    {
        c=0;
        for(int j=1;j<=len;j++)
        {
            long long s=a[j]*i+c;
            a[j]=s%100000;
            c=s/100000;
        }
        while(c > 0) 
        {
            len++;
            a[len]=c%100000;
            c/=100000;
        }
    }
    long long b[10010],f[10010];
    int lenb,lenf=0;
    memset(f,0,sizeof(f));
    for(int i=1;i<=aa.len;i++)
    {
        memset(b,0,sizeof(b));
        long long x=aa.a[i],ss=0;int g=i;
        for(int j=1;j<=len;j++)
        {
            long long y=a[j]*x+ss;
            if(j==len)
            {
                b[g]=y%100000;
                b[g+1]=y/100000;
                if(y/100000!=0) lenb=g+1;
                else lenb=g;
            }
            else 
            {
                b[g++]=y%100000;
                ss=y/100000;
            }
        }
        lenf=max(lenf,lenb);
        for(int i=1;i<=lenf;i++)
        {
            f[i]=f[i]+b[i];
        }
        int ee;
        for(int i=1;i<=lenf;i++)
        {
            ee=f[i]/100000;
            f[i]=f[i]%100000;
            f[i+1]+=ee;
        }
        if(f[lenf+1]!=0) lenf++;
    }
    //*********************************
    memcpy(&aa.a,&f,sizeof(f));
    aa.len=lenf;
}

void wori()
{
    nn.len=max(aa.len,nn.len);
    for(int i=1;i<=nn.len;i++)
    {
        nn.a[i]=nn.a[i]+aa.a[i];
    }
    int ee;
    for(int i=1;i<=nn.len;i++)
    {
        ee=nn.a[i]/100000;
        nn.a[i]=nn.a[i]%100000;
        nn.a[i+1]+=ee;
    }
    if(nn.a[nn.len+1]!=0) nn.len++;
}

void work()
{
    wocao();
    long long w=n*(n+1);
    long long c=0;
    for(int i=1;i<=nn.len;i++)
    {
        long long s=w*nn.a[i]+c;
        nn.a[i]=s%100000;
        c=s/100000;
    } 
    while(c > 0) 
    {
        nn.len++;
        nn.a[nn.len]=c%100000;
        c/=100000;
    }
    wogan(n+3,m);
    long long y=2*(n+1)*m;
    long long p=0;
    for(int i=1;i<=aa.len;i++)
    {
        long long h=y*aa.a[i]+p;
        aa.a[i]=h%100000;
        p=h/100000;
    } 
    while(p > 0) 
    {
        aa.len++;
        aa.a[aa.len]=p%100000;
        p/=100000;
    }
    //*******************
    wogan1(n+2,m-1);
    wori();
    for(int i=nn.len;i>=1;i--) 
    {
        if(i==nn.len) printf("%d",nn.a[i]);
        else printf("%05d",nn.a[i]);
    }
    printf("\n");
}

int main()
{
    init();
    readdata();
    work();
    return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
城市应急指挥系统是智慧城市建设的重要组成部分,旨在提高城市对突发事件的预防和处置能力。系统背景源于自然灾害和事故灾难频发,如汶川地震和日本大地震等,这些事件造成了巨大的人员伤亡和财产损失。随着城市化进程的加快,应急信息化建设面临信息资源分散、管理标准不统一等问题,需要通过统筹管理和技术创新来解决。 系统的设计思路是通过先进的技术手段,如物联网、射频识别、卫星定位等,构建一个具有强大信息感知和通信能力的网络和平台。这将促进不同部门和层次之间的信息共享、交流和整合,提高城市资源的利用效率,满足城市对各种信息的获取和使用需求。在“十二五”期间,应急信息化工作将依托这些技术,实现动态监控、风险管理、预警以及统一指挥调度。 应急指挥系统的建设目标是实现快速有效的应对各种突发事件,保障人民生命财产安全,减少社会危害和经济损失。系统将包括预测预警、模拟演练、辅助决策、态势分析等功能,以及应急值守、预案管理、GIS应用等基本应用。此外,还包括支撑平台的建设,如接警中心、视频会议、统一通信等基础设施。 系统的实施将涉及到应急网络建设、应急指挥、视频监控、卫星通信等多个方面。通过高度集成的系统,建立统一的信息接收和处理平台,实现多渠道接入和融合指挥调度。此外,还包括应急指挥中心基础平台建设、固定和移动应急指挥通信系统建设,以及应急队伍建设,确保能够迅速响应并有效处置各类突发事件。 项目的意义在于,它不仅是提升灾害监测预报水平和预警能力的重要科技支撑,也是实现预防和减轻重大灾害和事故损失的关键。通过实施城市应急指挥系统,可以加强社会管理和公共服务,构建和谐社会,为打造平安城市提供坚实的基础。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值