内存限制,赶紧去探望=v=

时间限制:10000ms

hihoCoder #1195 : 高斯消元·一,

岁月限制:10000ms 单点时限:1000ms 内存限制:256MB

日子范围:10000ms

单点时限:1000ms

描述

小Ho:喂不得了哇,这边便利店的薯片半价了!

小Hi:啥?!

小Ho:这边的便利店在促销让利啊。

小Hi:走走走,赶紧去看看=v=

于是小Hi和小Ho来到了便利店。

业主为了促销,推出了组合包的样式,将不同数额的各项商品打包成一个整合,顾客可以采取组合展开采购。比如2袋薯片,1听可乐的三结合只要5元,而1袋薯片,2听可乐的结合只要4元。

经过询问老总,小Hi和小Ho知道:一共有N种不同的货色和M种不同的货品组合;每一个结合的价格万分组合内货物售价之和,一个构成内一律件货物不会超过10件。

小Hi:这样算下来的话,一听可乐就是1元,而一包薯片是2元。小Ho,倘若你精晓所有的重组意况,你能分别算出每一件商品单独的标价么?

小Ho:当然可以了,这样的小问题怎么能难到自己吧?

指示:高斯消元

单点时限:1000ms

内存限制:256MB

输入

第1行:2个正整数,N,M。表示商品的数目N,组合的数量M。1≤N≤500, N≤M≤2*N

第2..M+1行:N+1个非负整数,第i+1行第j列表示在第i个组成中,商品j的数据a[i][j]。第i+1行第N+1个数表示该结合的售价c[i]。0≤a[i][j]≤10,
0≤c[i]≤10^9

内存限制:256MB

描述

小Ho:喂不得了呀,这边便利店的薯片半价了!

小Hi:啥?!

小Ho:这边的便利店在促销优惠啊。

小Hi:走走走,赶紧去探视=v=

于是乎小Hi和小Ho来到了便利店。

老总娘为了减价,推出了组合包的花样,将不同数量的各个货物打包成一个构成,顾客可以挑选组合展开购买。比如2袋薯片,1听可乐的三结合只要5元,而1袋薯片,2听可乐的结合只要4元。

透过通晓首席营业官,小Hi和小Ho知道:一共有N种不同的货物和M种不同的货物组合;每一个整合的价钱相当于组合内商品售价之和,一个重组内一样件货物不会超过10件。

小Hi:这样算下来的话,一听可乐就是1元,而一包薯片是2元。小Ho,如若您知道所有的整合情形,你能分别算出每一件商品单独的价位么?

小Ho:当然可以了,这样的小题目怎么能难到自家呢?

   

指示:高斯消元

 

输出

若没有章程总括出各样商品单独的价位,输出”No solutions”

若可能存在六个不等的结果,输出”Many solutions”

若存在唯一可能的结果,输出N行,每行一个非负整数,第i行表示第i个商品单独的售价。数据保证假诺存在唯一解,那么解一定恰好是非负整数解。

样例输入

2 2
2 1 5
1 2 4

样例输出

2

1

这坑爹oj没多少,害的自身拍了以深夜,

题相比简单,高斯消元的模版题,

注意eps要开double类型的

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<cmath>
 4 #include<algorithm>
 5 using namespace std;
 6 const int MAXN=1001;
 7 const double eps =1e-7;
 8 typedef double Matrix[MAXN][MAXN] ;
 9 inline void read(int &n)
10 {
11     char c=getchar();bool flag=0;n=0;
12     while(c<'0'||c>'9') c=='-'?flag==1,c=getchar():c=getchar();
13     while(c>='0'&&c<='9')    n=n*10+c-48,c=getchar();
14 }
15 int n,m;
16 Matrix a;
17 double t[MAXN];
18 void Gauss()
19 {
20     int r;
21     for(int i=1;i<=n;i++)// 枚举每一行 
22     {
23         r=i;
24         for(int j=m;j>=i+1;j--)// 枚举后面的行 
25             fabs(a[j][i])>fabs(a[r][i])?r=j:r=r;
26         if(r!=i)    swap(a[r],a[i]); 
27         if(r==i&&fabs(a[i][i])<eps)        printf("Many solutions"),exit(0);
28             
29         for(int j=i+1;j<=m;j++)// 行 
30         {
31             for(int k=n+1;k>i;k--)// 列 
32                 a[j][k]-=(a[j][i]/a[i][i])*a[i][k];
33             a[j][i]=0;
34         }
35     }
36     
37     for(int i=n,j;i<=m;i++)
38     {
39         for(j=1;j<=n;j++)        if(fabs(a[i][j])>eps)        break;
40         if(j==n+1&&fabs(a[i][n+1])>eps)        
41             printf("No solutions\n"),exit(0);
42     }// 是否有解 
43     
44     for(int i=n;i>=1;i--)// 枚举行 
45     {
46         for(int j=i+1;j<=n;j++)// 枚举列 
47             a[i][n+1]-=a[i][j]*a[j][n+1];
48         a[i][n+1]/=a[i][i];
49     }
50     for(int i=1;i<=n;i++)
51         printf("%d\n",(int)(a[i][n+1]+0.5));
52 }
53 int main()
54 {
55     freopen("a.in","r",stdin);
56     freopen("c.out","w",stdout);
57     read(n);read(m);
58     for(int i=1;i<=m;i++)
59         for(int j=1;j<=n+1;j++)
60             scanf("%lf",&a[i][j]);
61     Gauss();
62     
63     return 0;
64 }

 

    

http://www.bkjia.com/cjjc/1222801.htmlwww.bkjia.comtruehttp://www.bkjia.com/cjjc/1222801.htmlTechArticlehihoCoder #1195 : 高斯消元一,
时间限制:10000ms单点时限:1000ms内存限制:256MB 描述
小Ho:喂不得了哇,这边便利店的薯片半价了! 小Hi:啥?! 小…

描述

小Ho:喂不得了啊,这边便利店的薯片半价了!

小Hi:啥?!

小Ho:这边的便利店在优惠促销啊。

小Hi:走走走,赶紧去探望=v=

于是小Hi和小Ho来到了便利店。

总经理娘为了让利,推出了组合包的形式,将不同数量的各项商品打包成一个结合,顾客可以采用组合展开选购。比如2袋薯片,1听可乐的组成只要5元,而1袋薯片,2听可乐的结缘只要4元。

经过摸底首席执行官,小Hi和小Ho知道:一共有N种不同的货品和M种不同的货品组合;每一个构成的价格万分组合内商品售价之和,一个组成内一律件货物不会超过10件。

小Hi:这样算下来的话,一听可乐就是1元,而一包薯片是2元。小Ho,假诺您了然所有的结合情状,你能分别算出每一件商品单独的价钱么?

小Ho:当然可以了,这样的小问题怎么能难到本人吗?

   

提示:高斯消元

 

输入

第1行:2个正整数,N,M。表示商品的多少N,组合的数量M。1≤N≤500, N≤M≤2*N

第2..M+1行:N+1个非负整数,第i+1行第j列表示在第i个结合中,商品j的数目a[i][j]。第i+1行第N+1个数表示该结合的售价c[i]。0≤a[i][j]≤10,
0≤c[i]≤10^9

输入

第1行:2个正整数,N,M。表示商品的数据N,组合的数量M。1≤N≤500, N≤M≤2*N

第2..M+1行:N+1个非负整数,第i+1行第j列表示在第i个组成中,商品j的数额a[i][j]。第i+1行第N+1个数表示该组合的售价c[i]。0≤a[i][j]≤10,
0≤c[i]≤10^9

输出

若没有主意总结出各样商品单独的标价,输出”No solutions”

若可能存在五个不等的结果,输出”Many solutions”

若存在唯一可能的结果,输出N行,每行一个非负整数,第i行代表第i个商品单独的售价。数据保证假如存在唯一解,那么解一定恰好是非负整数解。

样例输入

2 2
2 1 5
1 2 4

样例输出

2

1

这坑爹oj没多少,害的我拍了以早上,

题相比简单,高斯消元的沙盘题,

注意eps要开double类型的

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<cmath>
 4 #include<algorithm>
 5 using namespace std;
 6 const int MAXN=1001;
 7 const double eps =1e-7;
 8 typedef double Matrix[MAXN][MAXN] ;
 9 inline void read(int &n)
10 {
11     char c=getchar();bool flag=0;n=0;
12     while(c<'0'||c>'9') c=='-'?flag==1,c=getchar():c=getchar();
13     while(c>='0'&&c<='9')    n=n*10+c-48,c=getchar();
14 }
15 int n,m;
16 Matrix a;
17 double t[MAXN];
18 void Gauss()
19 {
20     int r;
21     for(int i=1;i<=n;i++)// 枚举每一行 
22     {
23         r=i;
24         for(int j=m;j>=i+1;j--)// 枚举后面的行 
25             fabs(a[j][i])>fabs(a[r][i])?r=j:r=r;
26         if(r!=i)    swap(a[r],a[i]); 
27         if(r==i&&fabs(a[i][i])<eps)        printf("Many solutions"),exit(0);
28             
29         for(int j=i+1;j<=m;j++)// 行 
30         {
31             for(int k=n+1;k>i;k--)// 列 
32                 a[j][k]-=(a[j][i]/a[i][i])*a[i][k];
33             a[j][i]=0;
34         }
35     }
36     
37     for(int i=n,j;i<=m;i++)
38     {
39         for(j=1;j<=n;j++)        if(fabs(a[i][j])>eps)        break;
40         if(j==n+1&&fabs(a[i][n+1])>eps)        
41             printf("No solutions\n"),exit(0);
42     }// 是否有解 
43     
44     for(int i=n;i>=1;i--)// 枚举行 
45     {
46         for(int j=i+1;j<=n;j++)// 枚举列 
47             a[i][n+1]-=a[i][j]*a[j][n+1];
48         a[i][n+1]/=a[i][i];
49     }
50     for(int i=1;i<=n;i++)
51         printf("%d\n",(int)(a[i][n+1]+0.5));
52 }
53 int main()
54 {
55     freopen("a.in","r",stdin);
56     freopen("c.out","w",stdout);
57     read(n);read(m);
58     for(int i=1;i<=m;i++)
59         for(int j=1;j<=n+1;j++)
60             scanf("%lf",&a[i][j]);
61     Gauss();
62     
63     return 0;
64 }

 

    

输出

若没有章程总结出各类商品单独的标价,输出”No solutions”

若可能存在四个不同的结果,输出”Many solutions”

若存在唯一可能的结果,输出N行,每行一个非负整数,第i行代表第i个商品单独的售价。数据保证如果存在唯一解,那么解一定恰好是非负整数解。

样例输入

2 2
2 1 5
1 2 4

样例输出

2

1

那坑爹oj没多少,害的自身拍了以早上,

题相比较简单,高斯消元的模板题,

注意eps要开double类型的

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<cmath>
 4 #include<algorithm>
 5 using namespace std;
 6 const int MAXN=1001;
 7 const double eps =1e-7;
 8 typedef double Matrix[MAXN][MAXN] ;
 9 inline void read(int &n)
10 {
11     char c=getchar();bool flag=0;n=0;
12     while(c<'0'||c>'9') c=='-'?flag==1,c=getchar():c=getchar();
13     while(c>='0'&&c<='9')    n=n*10+c-48,c=getchar();
14 }
15 int n,m;
16 Matrix a;
17 double t[MAXN];
18 void Gauss()
19 {
20     int r;
21     for(int i=1;i<=n;i++)// 枚举每一行 
22     {
23         r=i;
24         for(int j=m;j>=i+1;j--)// 枚举后面的行 
25             fabs(a[j][i])>fabs(a[r][i])?r=j:r=r;
26         if(r!=i)    swap(a[r],a[i]); 
27         if(r==i&&fabs(a[i][i])<eps)        printf("Many solutions"),exit(0);
28             
29         for(int j=i+1;j<=m;j++)// 行 
30         {
31             for(int k=n+1;k>i;k--)// 列 
32                 a[j][k]-=(a[j][i]/a[i][i])*a[i][k];
33             a[j][i]=0;
34         }
35     }
36     
37     for(int i=n,j;i<=m;i++)
38     {
39         for(j=1;j<=n;j++)        if(fabs(a[i][j])>eps)        break;
40         if(j==n+1&&fabs(a[i][n+1])>eps)        
41             printf("No solutions\n"),exit(0);
42     }// 是否有解 
43     
44     for(int i=n;i>=1;i--)// 枚举行 
45     {
46         for(int j=i+1;j<=n;j++)// 枚举列 
47             a[i][n+1]-=a[i][j]*a[j][n+1];
48         a[i][n+1]/=a[i][i];
49     }
50     for(int i=1;i<=n;i++)
51         printf("%d\n",(int)(a[i][n+1]+0.5));
52 }
53 int main()
54 {
55     freopen("a.in","r",stdin);
56     freopen("c.out","w",stdout);
57     read(n);read(m);
58     for(int i=1;i<=m;i++)
59         for(int j=1;j<=n+1;j++)
60             scanf("%lf",&a[i][j]);
61     Gauss();
62     
63     return 0;
64 }

 

    

相关文章