2D розділена сукупна сума


16

Виклик

Задавши матрицю M із r рядками та c стовпцями та два булеві списки V довжини r та H довжини c , обчисліть розділену кумулятивну вертикальну та горизонтальну суми.

Правила

  • r і c більше або дорівнюють одиниці

  • H і V починаються з справжнього значення

  • Значення в M знаходяться в розумному числовому домені вашої мови.

  • Розбиття та підбиття підсумків починається у верхньому лівому куті.

Пройдіться

Дано M :

┌──────────────┐
│ 1  2  3  4  5│
│ 6  7  8  9 10│
│11 12 13 14 15│
│16 17 18 19 20│
└──────────────┘

H :1 0 1 0 0

V :1 1 0 1

Розділіть M на групи стовпців, починаючи нову групу при кожному справжньому значенні H

┌─────┬────────┐
│ 1  2│ 3  4  5│
│ 6  7│ 8  9 10│
│11 12│13 14 15│
│16 17│18 19 20│
└─────┴────────┘

Розбийте кожну групу стовпців на групи рядків, починаючи нову групу при кожному справжньому значенні V :

┌─────┬────────┐
│ 1  2│ 3  4  5│
├─────┼────────┤
│ 6  7│ 8  9 10│
│11 12│13 14 15│
├─────┼────────┤
│16 17│18 19 20│
└─────┴────────┘

Сумарне підсумовування кожної комірки по горизонталі:

┌─────┬────────┐
│ 1  3│ 3  7 12│
├─────┼────────┤
│ 6 13│ 8 17 27│
│11 23│13 27 42│
├─────┼────────┤
│16 33│18 37 57│
└─────┴────────┘

Сукупно підсумовуйте кожну клітинку вертикально:

┌─────┬────────┐
│ 1  3│ 3  7 12│
├─────┼────────┤
│ 6 13│ 8 17 27│
│17 36│21 44 69│
├─────┼────────┤
│16 33│18 37 57│
└─────┴────────┘

Результат:

┌──────────────┐
│ 1  3  3  7 12│
│ 6 13  8 17 27│
│17 36 21 44 69│
│16 33 18 37 57│
└──────────────┘

Додаткові тестові випадки

M :

┌───────────┐
│15 11 11 17│
│13 20 18  8│
└───────────┘

H : 1 0 0 1V :1 0

Результат:

┌───────────┐
│15 26 37 17│
│28 59 88 25│
└───────────┘

M :

┌─┐
│7│
└─┘

Результат ( H і V повинен бути 1):

┌─┐
│7│
└─┘

M :

┌──┐
│ 3│
│-1│
│ 4│
└──┘

V : 1 1 0( H повинно бути 1)

Результат:

┌──┐
│ 3│
│-1│
│ 3│
└──┘

M :

┌───────────────────────────────────────────────────────┐
│10    7.7 1.9 1.5 5.4  1.2 7.8 0.6 4.3 1.2  4.5 5.4 0.3│
│ 2.3  3.8 4.1 4.5 1    7.7 3   3.4 6.9 5.8  9.5 1.3 7.5│
│ 9.1  3.7 7.2 9.8 3.9 10   7.6 9.6 7.3 6.2  3.3 9.2 9.4│
│ 4.3  4.9 7.6 2   1.4  5.8 8.1 2.4 1.1 2.3  7.3 3.6 6  │
│ 9.3 10   5.8 9.6 5.7  8.1 2.1 3.9 4   1.3  6.3 3.1 9  │
│ 6.6  1.4 0.5 6.5 4.6  2.1 7.5 4.3 9   7.2  2.8 3.6 4.6│
│ 1.7  9.9 2.4 4.5 1.3  2.6 6.4 7.8 6.2 3.2 10   5.2 8.9│
│ 9.9  5.3 4.5 6.3 1.4  3.1 2.3 7.9 7.8 7.9  9.6 4   5.8│
└───────────────────────────────────────────────────────┘

H :1 0 0 1 0 1 1 1 0 1 1 1 0

V :1 0 0 0 0 1 0 0

Результат:

┌────────────────────────────────────────────────────────────────┐
│10   17.7 19.6  1.5  6.9  1.2  7.8  0.6  4.9  1.2  4.5  5.4  5.7│
│12.3 23.8 29.8  6   12.4  8.9 10.8  4   15.2  7   14    6.7 14.5│
│21.4 36.6 49.8 15.8 26.1 18.9 18.4 13.6 32.1 13.2 17.3 15.9 33.1│
│25.7 45.8 66.6 17.8 29.5 24.7 26.5 16   35.6 15.5 24.6 19.5 42.7│
│35   65.1 91.7 27.4 44.8 32.8 28.6 19.9 43.5 16.8 30.9 22.6 54.8│
│ 6.6  8    8.5  6.5 11.1  2.1  7.5  4.3 13.3  7.2  2.8  3.6  8.2│
│ 8.3 19.6 22.5 11   16.9  4.7 13.9 12.1 27.3 10.4 12.8  8.8 22.3│
│18.2 34.8 42.2 17.3 24.6  7.8 16.2 20   43   18.3 22.4 12.8 32.1│
└────────────────────────────────────────────────────────────────┘

Відповіді:


9

Желе , 10 байт

Zœṗ@+\€Ẏð/

Спробуйте в Інтернеті! та Останній тестовий випадокGкінцем для читабельності).

Введення приймається як список [M, H, V].

Пояснення

Zœṗ@+\€Ẏð/  Input: [M, H, V]
        ð/  Insert the previous (f) as a dyadic link
            Forms f( f(M, H) , V)
            For f(x, y):
Z             Transpose x
 œṗ@          Partition the rows of x^T at each true in y
    +\€       Compute the cumulative sums in each partition
       Ẏ      Tighten (Joins all the lists at the next depth)

Ви можете використовувати колонтитул , як це так , що вам не доведеться возитися з фактичним кодом.
Ерік Аутгольфер

7

APL (Діалог) , 13 байт

Приймає IST з множинності в якості аргументу.

{⍉⊃,/+\¨⍺⊂⍵}/

Спробуйте в Інтернеті!

{}/ Вставити (зменшити на) наступну анонімну функцію, де термін зліва представлений ⍺, а термін праворуч представлений ⍵. Оскільки функції APL є правильним асоціативним, це, таким чином, V f ( H f M ).

⍺⊂⍵ розділ ⍵ відповідно до ⍺

+\¨ сукупна сума кожної частини

,/ зменшити шляхом конкатенації (це додає результат для зменшення рангу)

 розкрити

 переносити


6

Python 2 + numpy, 143 138 117 115 110 108 байт

-21 байт завдяки Adám !

lambda M,*L:reduce(lambda m,l:vstack(map(lambda p:cumsum(p,0),split(m,*where(l)))).T,L,M)
from numpy import*

Спробуйте в Інтернеті!


1
попросіть розділити, розділити та закінчити один раз, перенести, повторити.
Адам

@ Adám Спасибі, я чомусь не думав про це.
notjagan

Мені сподобався пошук у списку двох функцій :)
Джонатан Аллан

2
Будь ласка, зробіть заголовок "Python 3 + numpy"
Leaky Nun

5

Желе ,  15  14 байт

œṗ+\€Ẏ
ḢçЀZð⁺

Діадичне посилання H,Vзліва і Mсправа і повернення отриманої матриці.

Спробуйте в Інтернеті!

Крім того, як окремий рядок також для 14: Ḣœṗ+\€Ẏ$¥Ð€Zð⁺

Як?

œṗ+\€Ẏ - Link 1: partition and cumSum: list of partition bools, list of values
œṗ     - partition (the values) at truthy indexes (of the bools)
    €  - for €ach part:
  +\   -   cumulative reduce by addition
     Ẏ - tighten (flattens back into a list)

ḢçЀZð⁺ - Main link: list of lists, [H,V]; list of lists, M
      ⁺ - perform this twice:
     ð  - [it's a dyadic chain for the second pass, first pass is dyadic implicitly]
Ḣ       -   head, pop it & modify (so H the first time, V the second)
  Ѐ    -   map across right: (M the first time, the intermediate result the second)
 ç      -     the last link (1) as a dyad
    Z   -   transpose the result (do the rows first time, and the columns the second)

Попередній:

œṗ@+\€Ẏ
ç€Zç€⁵Z

Повна програма, що друкує представлення результату.


Whoa -50% від попередньої відповіді желе!
Adám

Хто що? Ого. Мені справді потрібно вивчити, як ти це зробив ... Неймовірно порівняно з моїм!
HyperNeutrino

О, це робить два напрямки окремо, правда? Розумний.
HyperNeutrino

Я думаю, це робить приблизно те ж саме ...
Джонатан Аллан

Хороший метод. Значить, я можу це перемогти APL. У мене 14 байт.
Adám

4

MATL , 19 байт

,!ix"0GYs@12XQ!g]v!

Вхідними даними є M(матриця), H(вектор стовпців), V(вектор стовпців). Роздільник рядків є ;.

Спробуйте в Інтернеті! Або перевірити всі тестові випадки: 1 , 2 , 3 , 4 , 5 .

Пояснення

Це робить сукупну суму горизонтально, потім вертикально.

,          % Do the following twice
  !        %   First time this inputs M implicitly. Transpose. Second time
           %   it transposes the result of the horizontal cumulative sum
  ix       %   Input H (first time) or V (second time). Delete it; but gets
           %   copied into clipboard G
  "        %   For each column of the matrix
    0G     %     Push most recent input: H (first time) or V (second)
    Ys     %     Cumulative sum. This produces a vector of integer values
           %     such that all columns (first time) or rows (second) of M 
           %     with the same value in this vector should be cumulatively
           %     summed
    @      %     Push current column of M transposed (first time) or M after
           %     horizontal cumulative sum (second time)
    12XQ   %     Cumulative sum. Gives a cell array of row vectors
    !g     %     Join those vectors into one row vector
  ]        %   End
  v        %   Concatenate the row vectors vertically into a matrix
  !        %   Transpose. This corrects for the fact that each column vector
           %   of the matrix was cumulatively summed into a row vector
           % Implicit end. Implicit display

1
Найбільш вражаюче. Я здогадуюсь, що Matlab був свого роду зроблений для подібних речей.
Adám

@ Adám Я впевнений, що довжина APL не буде дуже різною :-)
Луїс Мендо,

Моя контрольна реалізація, що використовується для створення тестових випадків, становить 26 байт.
Adám

@ Адам Дарн! APL побиття желе? Це неприпустимо! (повинен гольф моє рішення ... lol) xD
HyperNeutrino

@HyperNeutrino Ну, желе не має як рангів, так і глибини, як APL та J.
Adám

3

J , 20 байт

;@(<@(+/\);.1|:)&.>/

Спробуйте в Інтернеті!

Введення приймається як масив ящиків, що містять [V, H, M].

Пояснення

;@(<@(+/\);.1|:)&.>/  Input: [V H M]
  (     g      )   /  Insert g and reduce (right-to-left)
                      Forms V g H g M = V g (H g M)
                & >     Unbox each
             |:         Transpose the right arg
          ;.1           Partition
      +/\               Reduce each prefix using addition (cumulative sum)
   <@                   Box each partition
;@                      Raze (Concatenate the contents in each box)
                &.>     Box the result

2

Mathematica, 212 байт

(T=Transpose;A=AppendTo;J=Flatten;f[s_]:=Block[{},t=2;r=1;w={};While[t<=Length@s,If[s[[t]]==0,r++,w~A~r;r=1];t++];w~A~r];K[x_,y_]:=Accumulate/@#&/@(FoldPairList[TakeDrop,#,f@y]&/@x);d=J/@K[#,#2];T[J/@K[T@d,#3]])&


вхід
[M, H, V]

[{{15, 11, 11, 17}, {13, 20, 18, 8}}, {1, 0, 0, 1}, {1, 0}]


2

C # (.NET Core) , 164 байти

(M,H,V)=>{int a=M.Length,b=M[0].Length,i,j;for(i=0;i<a;i++)for(j=0;j<b;j++)if(!H[j])M[i][j]+=M[i][j-1];for(i=0;i<a;i++)for(j=0;j<b;j++)if(!V[i])M[i][j]+=M[i-1][j];}

Спробуйте в Інтернеті!

В основному це робить саме так, як зазначено в ОП. Спочатку повторюється підсумовувати по горизонталі, а потім повторювати по вертикалі.


2

Haskell , 129 байт 119 байт

s m v=tail$scanl(\a(x,s)->if s then x else zipWith(+)a x)[](zip m v)
t=foldr(zipWith(:))$repeat[]
f m h v=t$s(t$s m v)h

Спробуйте в Інтернеті!

Збережено 10 байт завдяки @ceasedtoturncounterclockwis

t(для транспонування) перемикає рядки та стовпці. Швидке пояснення:

foldr(zipWith(:))(repeat[])(r1,...,rn) =
zipWith(:) r1 (zipWith(:) r2 (... zipWith(:) rn (repeat [])))

Читайте справа наліво: ми переглядаємо рядки знизу вгору та натискаємо кожне значення у стовпчику призначення.

s в основному є сукупною сумою векторів, але скидається, коли виникає значення True v

fпідсумовує рядки з sнаступними vі робимо те ж саме, що й стовпці, що слідуютьh


t=foldr(zipWith(:))(repeat[]). Не тільки коротше, але й значно менш неефективно.
перестали повертати проти годинника,

@ceasedtoturncounter clockwis Дякую за пораду.
jferard

1

JavaScript (ES6), 88 байт

(a,h,v)=>a.map(b=>b.map((e,i)=>t=h[i]?e:t+e)).map((b,j)=>t=v[j]?b:t.map((e,i)=>e+b[i]))

0

Желе , 31 байт

+\€€
œṗḊZ€⁵œṗ$€Ḋ€Ç€ÇZ€€Z€;/€€;/

Спробуйте в Інтернеті!

Gah, це занадто довго для Jelly xD

BTW, 11/31 байт у цій програмі складається з символів євро. Це більше третини програми!


Занадто багато євро.
Adám

@ Adám Мої думки саме: P Робота з матрицями з подвійною розподіленою формою не така весела, як я думав, що це буде, тому що я займаюся
картою

Чому ти витрачаєш свої гроші так € - €
Куртуа
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.