Label-free imaging flow cytometry is a powerful tool for biological and medical research as it overcomes technical challenges in conventional fluorescence-based imaging flow cytometry that predominantly relies on fluorescent labeling. To date, two distinct types of label-free imaging flow cytometry have been developed, namely optofluidic time-stretch quantitative phase imaging flow cytometry and stimulated Raman scattering (SRS) imaging flow cytometry. Unfortunately, these two methods are incapable of probing some important molecules such as starch and collagen. Here, we present another type of label-free imaging flow cytometry, namely multiphoton imaging flow cytometry, for visualizing starch and collagen in live cells with high throughput. Our multiphoton imaging flow cytometer is based on nonlinear optical imaging whose image contrast is provided by two optical nonlinear effects: four-wave mixing (FWM) and second-harmonic generation (SHG). It is composed of a microfluidic chip with an acoustic focuser, a lab-made laser scanning SHG-FWM microscope, and a high-speed image acquisition circuit to simultaneously acquire FWM and SHG images of flowing cells. As a result, it acquires FWM and SHG images (100 × 100 pixels) with a spatial resolution of 500 nm and a field of view of 50 μm × 50 μm at a high event rate of four to five events per second, corresponding to a high throughput of 560-700 kb/s, where the event is defined by the passage of a cell or a cell-like particle. To show the utility of our multiphoton imaging flow cytometer, we used it to characterize Chromochloris zofingiensis (NIES-2175), a unicellular green alga that has recently attracted attention from the industrial sector for its ability to efficiently produce valuable materials for bioplastics, food, and biofuel. Our statistical image analysis found that starch was distributed at the center of the cells at the early cell cycle stage and became delocalized at the later stage. Multiphoton imaging flow cytometry is expected to be an effective tool for statistical high-content studies of biological functions and optimizing the evolution of highly productive cell strains.