Cody

# Problem 42503. Generating random matrix with given probability mass function

Solution 2555481

Submitted on 14 Jun 2020 by Rafael S.T. Vieira
This solution is locked. To view this solution, you need to provide a solution of the same size or smaller.

### Test Suite

Test Status Code Input and Output
1   Pass
rnd = sort(rand(randi([10,20]),1)); prob = vertcat(rnd(1,:),diff(rnd,1,1),1-rnd(end,:)); sz = [1 1e5;1e5 1;1e3 1e2;randi([100 200], 100, 2)]; sz = sz(randi(size(sz,1)),:); x = rndsampling(sz(1),sz(2),prob); prob_est = histcounts(x,1:numel(prob)+1,'Normalization','probability').'; err = mean(abs(prob_est - prob)) assert(err < 0.005 && isequal(size(x),sz) && all(~isnan(x(:))));

964 930 635 787 537 15 203 652 993 297 567 672 840 288 961 335 622 515 999 95 660 353 431 467 162 414 840 49 442 36 616 150 298 329 222 336 711 140 345 422 691 755 951 4 523 528 657 773 471 583 560 488 370 734 946 384 967 656 947 364 942 357 978 734 514 282 906 772 88 207 718 707 410 518 859 779 572 809 486 77 278 207 746 766 472 508 232 730 473 318 596 19 876 816 923 400 426 981 629 669 300 765 867 613 386 919 806 157 248 591 382 875 848 399 632 535 958 159 933 357 645 342 83 817 628 559 852 545 299 60 483 66 504 316 631 840 886 412 589 158 470 109 986 83 176 312 277 282 642 685 738 572 609 51 820 358 840 487 284 846 263 327 166 257 985 851 576 869 516 371 663 247 759 18 668 226 621 257 295 514 757 665 527 216 644 501 259 320 565 95 600 80 604 799 944 52 987 676 353 948 317 849 436 766 117 302 577 314 561 463 387 218 30 518 443 604 799 438 150 992 965 709 89 603 432 794 902 952 912 175 737 753 438 234 174 331 532 516 276 147 797 264 7 877 259 289 701 545 495 722 655 261 253 718 358 811 686 437 672 101 523 237 790 322 904 206 787 304 33 931 735 117 354 672 750 420 593 578 811 936 188 122 482 395 342 636 822 767 225 936 417 549 803 328 107 71 641 866 329 713 117 822 705 812 672 201 658 11 722 273 354 76 364 312 49 894 874 869 729 599 681 412 206 29 788 923 416 400 210 28 912 56 223 268 879 214 224 969 761 61 751 801 240 381 337 259 71 495 881 528 956 654 100 792 585 580 557 233 213 596 339 733 527 80 309 628 133 174 570 381 777 873 406 442 104 691 745 955 448 160 45 778 240 764 526 650 890 448 969 716 932 795 708 572 790 240 630 615 171 706 792 67 418 518 137 44 590 585 118 862 363 550 580 303 59 352 707 599 198 965 544 669 553 895 134 811 233 401 707 723 204 436 732 15 805 483 512 762 980 923 868 81 915 637 214 208 265 628 739 221 60 963 373 64 229 554 812 454 268 335 765 68 969 922 876 837 941 993 394 8 219 68 497 763 709 475 79 245 860 917 145 23 623 885 387 326 653 560 722 1000 16 378 8 539 124 304 333 958 117 753 938 998 881 98 856 550 36 713 19 764 983 217 696 812 399 755 332 595 589 650 714 542 518 681 937 752 228 406 977 299 282 590 147 225 945 846 32 271 677 953 358 797 695 730 805 355 394 328 256 291 650 723 916 462 605 858 890 167 175 174 669 973 171 156 944 6 924 146 405 547 860 739 436 75 173 599 793 617 290 443 198 822 20 958 19 250 599 145 559 741 284 444 665 549 83 34 951 273 136 999 949 498 200 594 124 744 716 335 423 421 486 84 440 93 606 260 131 694 454 639 987 65 838 838 343 426 884 710 684 529 275 253 945 647 844 969 227 179 103 715 852 446 742 60 535 38 960 574 491 472 373 949 780 671 243 774 67 771 723 377 93 981 693 731 540 196 841 348 62 701 976 237 760 350 558 399 129 605 958 701 59 2 248 632 596 133 573 619 47 523 234 958 409 327 368 881 470 785 387 138 80 705 955 144 436 27 781 430 544 562 675 421 376 951 614 211 88 784 673 214 511 426 853 881 655 795 395 298 25 926 687 331 434 76 655 571 31 38 119 836 791 974 347 758 813 850 902 915 463 232 139 530 974 733 684 450 521 400 548 506 792 652 770 661 540 422 729 139 591 927 180 498 471 995 671 806 365 257 533 700 227 664 559 431 602 727 697 138 871 620 623 469 649 913 151 86 767 862 109 520 230 479 525 802 624 280 536 308 512 969 770 390 769 173 270 511 712 887 457 560 115 774 894 554 18 926 367 79 715 655 69 281 735 391 730 524 937 30 178 178 763 809 226 904 373 768 509 669 511 87 594 246 489 270 442 300 975 947 23 818 712 110 477 977 486 93 452 227 166 97 684 442 209 830 100 923 586 572 411 292 669 648 414 501 247 258 472 203 913 320 725 866 368 439 36 341 848 359 584 942 938 766 828 951 99 21 487 830 822 349 63 564 489 306 809 554 235 321 406 824 682 659 643 866 958 275 306 103 972 896 3 645 321 571 848 814 726 202 482 763 797 147 761 651 653 704 139 876 509 331 197 181 752 280 144 843 748 732 888 284 909 283 870 637 178 459 943 706 582 172 581 40 755 387 350 173 227 134 664 282 631 181 446 85 139 39 461 100 37 58 974 434 200 691 905 942 471 231 776 725 534 818 849 321 255 782 480 101 590 2 373 366 971 194 650 831 197 151 385 869 749 235 180 28 702 105 240 345 636 378 873 689 603 833 877 620 587 695 634 696 999 799 711 850 626 261 872 940 615 608 267 561 417 197 246 734 52 882 352 987 442 660 299 512 974 630 784 296 505 826 709 621 168 672 711 814 51 244 609 725 202 40 848 755 820 989 852 174 956 862 9 277 948 89 628 794 32 368 681 517 162 529 198 661 338 52 390 896 583 87 336 68 790 841 13 726 505 260 8 34 843 97 555 874 44 819 260 689 152 386 840 493 403 469 73 971 700 645 492 839 724 669 833 876 365 764 816 159 744 520 826 537 746 235 787 877 911 834 625 518 491 236 150 55 858 296 90 483 435 590 149 145 620 984 147 858 802 957 204 969 130 145 808 191 175 774 126 835 421 418 459 18 538 706 756 382 325 789 505 156 319 643 374 890 861 966 794 735 298 414 173 158 412 268 474 984 156 918 159 625 477 973 238 417 498 807 946 983 176 222 650 110 591 664 179 756 138 769 460 221 974 182 121 211 492 381 851 862 492 69...