HUNNU-10307-最优分解问题




# include <queue>
# include <cstdio>
# include <cstring>
# include <iostream>
# include <algorithm>
using namespace std;

int v[4] = {1,10,100,1000},a[100000];

struct BigNumber{
    int d[1000];
    BigNumber() { memset(d,0,sizeof(d));    d[0] = 1; }
    BigNumber operator =(const int b){
        char s[40];
        sprintf(s,"%d",b);
        int len = strlen(s);
        memset(d,0,sizeof(d));
        d[0] = (len - 1)/4 + 1;
        for(int i = 0;i < len;i++)  d[(len-i-1)/4+1] += (s[i]-'0')*v[(len-i-1)%4];
        return *this;
    }
    void put(){
        printf("%d",d[d[0]]);
        for(int i = d[0]-1;i >= 1;i--)  printf("%04d",d[i]);
        printf("\n");
    }
};

BigNumber operator * (const BigNumber &a,const int &k)
{
    BigNumber c;
    c.d[0] = a.d[0];
    int i ,x = 0;
    for(int i = 1;i <=a.d[0];i++){
        x = a.d[i]*k + x;
        c.d[i] = x%10000;
        x /= 10000;
    }
    while(x > 0){
        c.d[++c.d[0]] = x%10000;
        x /= 10000;
    }
    while((c.d[0]>1) && (c.d[c.d[0]] == 0)) --c.d[0];
    return c;
}

int main()
{
    //freopen("a.txt","r",stdin);
    //freopen("b.txt","r",stdin);
    __int64 n,i;
    while(~scanf("%I64d",&n)){
        BigNumber ans;
        ans = 1;
        if(n <= 4){
            printf("%d\n",n);   continue;
        }
        //算出n的连续和;
        for(i = 2;i < n;i++){
            a[i-2] = i;
            n -= i;
        }
        
        int cnt = i - 3;
        while(n){
            a[cnt]++;
            n--;
            cnt -= 1;
            if(cnt == -1)   cnt = i - 3;
        }
        for(int j = 0;j <= i - 3;j++){
            //printf("%d ",a[j]);
            ans = ans * a[j];
        }
        ans.put();
    }
    return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
if __name__ == '__main__': # -------------Adjustable global parameters---------- n=512 # pixel number m=10 # number of time phases angle = 5 # #sample points = 360/angle on the boundary numOfAngles = int(180/angle) numOfContourPts = int(360/angle) labelID = 1 # 勾画的RS文件中第几个轮廓为GTV # path of the input data folder = 'E:\\MedData\\4DCT-202305\\' #patient = '0007921948' # 缺少时间信息 patient = '0000726380' # 病人的编号 # 呼吸曲线数据文件 vxpPath = folder+patient+'\\0000726380\\0000726380_20230420_143723.vxp' # Save the generated figures to the latex file path figPath = "D:\\HUNNU\\Research\\DMD\\4D-CT\\latex-DMD插值\\modify202305\\figure\\" # -------------Auto generated global parameters---------- # 每个dicom文件包含多少横截面 name = os.listdir(folder+patient+'\\0') cuts = [] for i in range(len(name)): if 'CT' in name[i][0:2]: cuts.append(i+1) cuts = np.array(cuts) # phase name times = np.linspace(0,90,10) # image pixel coordinate nums = np.linspace(0,n-1,n) x,y = np.meshgrid(nums,nums) # 输出dicom头文件信息 filename = folder+patient+'\\0\\CT.{}'.format(patient)+'.Image 1.dcm' print('CT dicom file information:') info = loadFileInformation(filename) # 像素之间的间距,包括列间距和行间距,单位mm SliceThickness = info['SliceThickness'] # Z轴的扫描分辨率,单位mm pixelSpace = info['pixelSpace'] # 一个像素所占的实际体积 pixelVol = float(pixelSpace[0])*float(pixelSpace[0])*float(SliceThickness) print('sliceThickness=',SliceThickness,' pixelSpace=',pixelSpace)
06-02

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值